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’s) and Computational Fluid Dynamics (CFD) 


have in recent years become equally important in accurate speed prediction for sailing yachts. 


ABSTRACT 


In this paper an overview is presented of the state-of-the-art 
of predicting the speed of sailing yachts. After explaining the 
importance of this topic in the design and racing of sailing 
yachts, the basic hydrodynamic and aerodynamic forces and 
their inter-dependence, which ultimately determine the speed 
potential at the various points of sail, are described. The dif- 
ferent possibilities of determining these forces, together with 
their centers of effort, is the principal subject of this paper. 


The use of model testing techniques and associated data 
reduction methods for the determination of the various hydro- 
dynamic forces are first of all dealt with. Details of how these 
model tests are best carried out are given. A method is defined 
with which it is possible to divide the total measured hydrodyn- 
amic resistance into its six basic components. Extrapolation of 
model test results to arrive at full-scale values, and the use 
thereof for speed prediction, are described. 


A major part of the paper is associated with the state-of-the- 
art of calculating sailing yacht performance as this is per- 
formed in semi-empirical models or so-called velocity predic- 
tion programs (VPP’s). Each of the elements involved are con- 
sidered, e.g. the side forces on hull and appendages, induced 
resistance, viscous and wave resistance, the resistance due to 
heel, the added resistance due to waves, etc. 


The recent availability of computational fluid dynamics 
programs and their present and future impact on sailing yacht 
performance prediction is addressed in the final part of the 


paper. 
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NOMENCLATURE 
A; = Area of jib or genoa; 
An = Area of mainsail; 


An = Normalized sail area; 

Ag, = Area of spinnaker; 

A, = Lateral area of body, normally expressed as 
maximum span times average chord length; 

Aix = Lateral area of keel; 

Aur = Total lateral area of canoe body and 
appendages; 

AP = Projected area in wind resistance calculation; 

AR = Geometric aspect ratio of lifting surface, equal 
to the ratio of the geometric span of the quarter- 
chord line divided by the average chord length; 

AR, = Effective aspect ratio (in many cases considered 
equal to two times the geometric aspect ratio); 

AR. = Effective aspect ratio of rudder, including the 
effect of the free surface; 

AR,s = Effective aspect ratio of sail plan; 

Awp = Area of waterplane of canoe body; 

BKT;, = Blanketing factor for jib or genoa; 

BKT,, = Blanketing factor for mainsail; 

BM = Metacentric radius of total underwater con- 
figuration in the zero heel condition; 

Binax = Maximum beam of deck; 

BMcs = Metacentric radius of canoe body in the zero 
heel condition; 

Bw. = Maximum beam of canoe body on waterline; 

Cz = Block coefficient of canoe body; 

Cp = Drag coefficient; 

Coro = Parasitic drag coefficient of hull and crew ex- 
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posed to the wind; 

Parasitic drag coefficient of jib or genoa; 
Parasitic drag coefficient of fully-battened main- 
sail; 

Parasitic drag coefficient of mainsail; 

Parasitic drag coefficient of spinnaker; 
Coefficient of frictional resistance of the subject 
body, usually considered equivalent to that of a 
flat plate with the same Reynolds number; 
Coefficient of frictional resistance of canoe 
body; 

Coefficient of frictional resistance of keel; 
Coefficient of resultant force; 

Lift coefficient; 

Coefficient of reduction in lift due to attitude of 
hull, usually associated with flow separation, 
otherwise not accounted for; 

Lift coefficient of jib or genoa; 

Lift coefficient of fully-battened mainsail; 

Lift coefficient of mainsail; 

Lift coefficient of spinnaker; 

Coefficient of linear component of lift; 
Coefficient of non-linear component of lift; 
Coefficient of lift at zero angle of leeway; 

Area coefficient of transverse section of canoe 
body with greatest area; 

Prismatic coefficient of canoe body; 

Resistance coefficient; 

Coefficient of resistance associated with the lift- 
attitude of the hull in causing flow separation or 
other behaviour, otherwise not accounted for; 
Coefficient of induced resistance; 

Coefficient of induced resistance associated with 
non-linear lift; 

Coefficient of residual stability; 

Coefficient of total resistance; 

Coefficient of total hydrodynamic resistance of 
canoe body; 

Coefficient of total resistance at zero heel and 
leeway; 

Coefficient of resistance associated with heel 
angle at zero side force; 

Coefficient of viscous resistance associated with 
heel angle at zero side force; 

Coefficient of resistance associated with flow 
separation; 

Coefficients in formula for wetted surface area 
of canoe body; 

Side force coefficient; 

Coefficient of linear component of side force; 
Coefficient of non-linear component of side 
force; 

Coefficient of viscous resistance; 

Coefficient of viscous pressure resistance; 
Coefficient of wave resistance; 

Area coefficient of waterplane of canoe body; 
Ratio of flap area to total area of lifting surface; 
Force in general; 

Froude number (Fy= V/(g.Ly)!” ); 

Froude number based on average chord length 
of the rudder; 

Resultant force; 

Lever of the hydrostatic righting moment of 
total configuration; 

Height of wind-sensing instruments at the mast- 
head above the water surface; 

Height of fore triangle of sail plan; 

Base of fore triangle of sail plan; 
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Vertical distance of the center of buoyancy of 
the complete underwater configuration from the 
bottom of the canoe body or the bottom of the 
deepest appendage, for the zero heel condition; 
Vertical distance of the center of buoyancy of 
appendage i from the bottom of the canoe body 
or the bottom of the deepest appendage, for the 
zero heel condition; 

Vertical distance of the center of buoyancy of 
the canoe body from the bottom thereof or the 
bottom of the deepest appendage, for the zero 
heel condition; 

Vertical distance of the center of gravity from 
the bottom of the canoe body or the bottom of 
the deepest appendage, for the zero heel con- 
dition; 

Separation drag coefficient of jib or genoa; 
Separation drag coefficient of mainsail; 
Separation drag coefficient of spinnaker; 

Lift force; 

Loss in lift due to lift-attitude of hull, usually 
associated with flow separation, otherwise not 
accounted for; 

Factor accounting for the effect of the bulb at 
the tip of a lifting surface on the lift thereof; 
Longitudinal position of center of buoyancy, in 
percent of Le, aft of the mid-Lc, location; 
Length of canoe body on waterline; 

Factor accounting for the effect of the canoe 
body on the lift of an appendages (keel, rudder, 
etc); 

Factor accounting for the effect of the canoe 
body on the lift of the keel; 

Total lift on keel, including that induced by 
canoe body and bulb and the lift associated with 
a trim tab, but excluding the carry-over lift on 
the canoe body and bulb (if fitted); 

Length of hull on waterline; 

Linear component of lift; 

Non-linear component of lift; 

Lift at zero angle of leeway; 

Total lift on lifting surface, including that in- 
duced by canoe body and bulb and the lift as- 
sociated with a flap or a trim tab, but excluding 
the carry-over lift on the canoe body and bulb 
(if fitted); 

Righting moment of crew sitting on weather 
rail; 

External heel moment exerted on model to 
obtain specific heel angle during model test; 
Residual stability lever divided by sing; 
Resistance associated with the attitude of the 
hull in generating lift, otherwise not accounted 
for; 

Frictional resistance; 

Induced resistance; 

Resistance of boundary layer interaction at the 
junction of appendage and canoe body; 

Induced resistance associated with non-linear 
component of lift; 

Reynolds number of canoe body; 

Residual resistance; 

Total resistance; 

Total resistance of canoe body; 

Total hydrodynamic resistance at zero heel; 
Viscous resistance; 

Viscous resistance of canoe body; 

Viscous resistance of keel; 
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Viscous resistance of rudder; 

Wave resistance; 

Total resistance associated with heel angle at 
zero side force; 

Viscous resistance associated with heel angle at 
zero side force; 

Wetted surface area of body; 

Wetted surface area of canoe body; 

Side force; 

Total loss in hydrodynamic side force due to lift 
attitude of hull, usually associated with flow 
separation, otherwise not accounted for; 

Side force of keel; 

Linear component of side force; 

Non-linear component of side force; 

Wetted surface area of keel; 

Factor accounting for additional drag for ex- 
posed mast of fractional rig; 

Factor accounting for the reduced lift of main- 
sail for fractional rigs; 

Thrust; 

Maximum draft of canoe body at zero heel; 
Draft of canoe body at the location of the 
quarter-chord line of keel; 

Effective draft of the canoe body in calculation 
Of Leg fact 5 

Effective draft values of the canoe body in cal- 
culation of Tea, ; 

Maximum draft of hull to the bottom of deepest 
appendage at zero heel; 

Taper ratio of lifting appendage, defined as the 
ratio of the chord length at the tip of the lifting 
surface to that at the root; 

Flow velocity; 

Apparent wind velocity in horizontal plane; 
Apparent wind velocity in plane normal to the 
mast; 

Boat velocity; 

Component of boat velocity opposite to true 
wind velocity (called the speed-made-good to 
windward); 

Minimum speed; 

True wind velocity in horizontal plane; 

True wind velocity in the plane normal to the 
mast at the masthead; 

True wind velocity in plane normal to the mast; 
Geometric span to the tip of the lifting surface 
where the quarter-chord line intersects the tip; 
Effective span of lifting surface; 

Effective span of keel, including effect of wing- 
lets or centerboard; 

Geometric span of keel; 

Average chord length of lifting surface; 

Average chord length of keel; 

Constant in definition of coefficient of lift reduc- 
tion, usually associated with flow separation, 
otherwise not accounted for, due to lift-attitude 
of hull; 

Constant in the formula for the linear lift coef- 
ficient depending on the shape of the tip; 
Constant in definition of coefficient of non-lin- 
ear lift; 

Constant in the formula for the non-linear lift 
coefficient depending on the shape of the tip; 
Average chord length of rudder; 

Constant in definition of coefficient of induced 
resistance due to non-linear lift; 

Length of root chord of lifting surface, exclu- 
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Constant in formula for separation resistance 
coefficient according to Prohaska; 

Length of tip chord of lifting surface; 

Length of tip chord of keel; 

Factor accounting for reduced lift of sail due to 
sail flattening; 

Free-surface influence factor indicating whether 
or not the free-surface will effect the lift on a 
lifting surface; 

Free-surface influence factor for the rudder; 
Acceleration due to gravity (9.81 m/sec? ); 
Maximum (vertical) height of bulb; 

Depth of top of appendage below water surface 
at the speed of interest; 

Depth of top of rudder below water surface at 
the speed of interest; 

Height of roughness used for tripping boundary 
layer; 

Factor accounting for effect of viscosity on the 
slope of the two-dimensional, linear lift coef- 
ficient curve; 

Factor accounting for effect of viscosity on flap- 
effectiveness factor; 

Mass of crew sitting on weather rail; 

Factor accounting for the reduced hoist of sail 
when reefing; 

Ratio of maximum radius of fillet to chord 
length, at root of appendage; 

Average thickness-chord ratio of lifting surface; 
Thickness-chord ratio of lifting surface at the 
root without the thickness of fillets; 

Average thickness-chord ratio of keel; 

Maximum thickness of appendage at the root, 
excluding the thickness of fillets; 

Taylor wake fraction; 

Horizontal distance of quarter-chord line of keel 
from the origin of the x,y,z axis system adopted 
(usually the forward or aft ending of the water- 
line length of the canoe body); 

Horizontal distance of quarter-chord line of 
rudder from the origin of the x,y,z axis system 
adopted {usually the forward or aft ending of the 
waterline length of the canoe body); 

Height of center of resultant aerodynamic side 
force above bottom of canoe body or deepest 
appendage; 

Height of center of resultant hydrodynamic side 
force above bottom of canoe body or deepest 
appendage; 

Height of center of hydrodynamic side force on 
appendage i above bottom of canoe body or 
deepest appendage; 

Height of tow force exerted on model above 
bottom of canoe body or deepest appendage; 
Angle-of-attack of a body to the incident flow; 
Effective deadrise angle of the canoe body in 
way of the quarter-chord line of the lifting 
surface, relative to a point a distance of 0.5c up 
from where the top of the quarter chord line 
intersects the canoe body, measured along the 
girth thereof; 

Induced flow angle at the lifting surface; 
Induced flow angle at the quarter-chord location 
of the rudder due to the downwash of the keel; 
Additional angle-of-attack of average section of 
lifting surface relative to reference plane (usu- 
ally a plane parallel to plane of symmetry of 
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yacht); 

Average zero-lift angle-of-attack of (cambered) 
lifting surface; 

Leeway angle; 

Apparent wind angle relative to centerline of 
yacht in horizontal plane; 

Apparent wind angle relative to centerline of 
yacht in plane normal to the mast; 

Effective leeway angle of average section of 
lifting surface in the case this has camber or 
additional angle-of-attack, or both, such as in 
the case of a leeboard; 

True wind angle relative to centerline of yacht 
in horizontal plane; 

True wind angle relative to centerline of yacht 
in plane normal to the mast; 

Ratio of the effective to the geometric span of a 
lifting surface, accounting for the effect of the 
canoe body on the lift; 

Ratio of the effective to the geometric span of a 
lifting surface, accounting for the effect of a 
bulb fitted at the tip of the lifting surface on the 
lift; 

Ratio of the effective to the geometric span of 
keel, accounting for the effect of the canoe body 
on the lift; 

Total weight of displacement of canoe body and 
appendages; 

Total weight of displacement of canoe body; 
Additional side force due to flap angle; 

Rotation angle of flap of lifting surface relative 
to centerline of yacht; 

Rotation angle of keel relative to centerline of 
yacht; 

Hydrodynamic drag angle; 

Aerodynamic drag angle; 

Sweep angle of quarter-chord line of lifting sur- 
face; 

Sweep angle of hinge line of flap of lifting sur- 
face; 

Scale factor of model; 

Kinematic viscosity of water at prevailing temp- 
erature; 

3.141592653; 

Density of water, as a function of temperature; 
Density of air, as a function of temperature; 
Half of average trailing-edge angle of lifting 
surface; 

Heel angle; 

Three-dimensional form factor accounting for 
the additional resistance of the subject body over 
that of the equivalent flat plate; 
Three-dimensional form factor of canoe body; 
Three-dimensional form factor of canoe body at 
heel angle ¢; 

Three-dimensional form factor of keel; 
Correction factor in formula for induced resis- 
tance for non-elliptical spanwise lift distribution; 
Flap-effectiveness factor, equal to slope of the 
angle-of-attack curve with respect to the flap 
angle 6; 

Effective slope of linear lift curve with respect 
to angle-of-attack of canoe body and appen- 
dages, when considered as a lifting body; 
Lift-curve slope with respect to angle-of-attack 
of lifting surface, assumed to be constant over 
the range in angle-of-attack of interest; 
Lift-curve slope with respect to angle-of-attack 


of keel, assumed to be constant over the range 
in angle-of-attack of interest; 


Vv = Volume of total underwater configuration; 
V api = Volume of appendage i; 
Ven = Volume of displacement of canoe body. 


Notes on subscripts: 

- Values for model scale are denoted by the subscript (M); 

- The subscript TH refers to the total hydrodynamic value for 
canoe body and appendages; 

- The subscript TA refers to the total aerodynamic value for 
the yacht; 

- The subscripts CB, K, R and W denote canoe body, keel, 
rudder and winglets, respectively; 

- The subscript "math" refers to mathematical model; 

- The subscripts u and d in the formulae for aerodynamic sail 
forces refers to upwind and downwind, respectively. 


1 INTRODUCTION 


Sailing yacht technology has improved over the last 10 years 
at a significant rate. In 1983 “Australia II" won the America’s 
Cup using advanced model testing and mathematical modelling 
techniques. For the conditions prevailing off Newport, Rhode 
Island, a near-perfect hull and (winged) keel concept had been 
developed within the confines of the applicable design rules, 
using superior technology. Ever since, the development of 
mathematical models for the prediction of the performance of 
sailing yachts and the further development of (model) testing 
techniques and associated data reduction methods have been 
pursued intensively, particularly in the context of America’s 
Cup competition. 


It is the ability of accurately predicting the speed of a sailing 
yacht in particular that has lead to the development of better 
hull forms and more efficient appendages. This ability will 
eventually also lead to equitable racing between yachts of 
different size and proportions. The recent world-wide introduc- 
tion of the so-called International Measurement System (IMS), 
developed in the United States, is a start in this direction. It is 
based on a semi-empirical model to predict the speed of dif- 
ferent types of sailing yachts and although far from perfect, 
represents a completely new approach to equitable yacht racing 
and handicapping. 


In this paper an overview is presented of the state-of-the-art 
of predicting the speed potential of sailing yachts. The suc- 
cessful execution thereof is closely associated with the ability 
to determine the various forces acting on a specific configu- 
ration, together with their centets of effort. Different ex- 
perimental and numerical methods are presently available for 
this purpose, each with its specific problems and deficiencies. 
It should be emphasized that the determination of these forces, 
their centers of effort and their inter-dependence, when the 
yacht makes heel and leeway, constitutes one of the most 
difficult problems in naval architecture today, particularly 
when the effect of waves on the performance has to be con- 
sidered as well. 


A full account is given of the model testing methods that can 
be adopted in determining the speed of a sailing yacht. Since 
Davidson’s memorable. towing tank investigations (Davidson, 
1936), the development of hydrodynamic testing techniques 
and data reduction methods have been considerably refined. 
Particular attention is focussed on important test data obtained. 


A major part of the paper is associated with the state-of-the- 
art of calculating sailing yacht performance as this is executed 
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in so-called Velocity Prediction Programs (VPP). Each of the 
parts involved herein are considered, viz: the model for the 
calculation of the side forces on canoe body and appendages, 
the induced resistance, viscous and wave resistance, the resis- 
tance due to heel, the added resistance due to waves, the 
aerodynamic forces on the various sails and the above-water 
hull and spars, etc. The final calculation of the various entities 
from force and moment equilibrium considerations, such as 
boat speed, heel angle, leeway angle and rudder angle, are 
dealt with. 


The recent availability of Computational Fluid Dynamics 
(CFD) programs and their impact on sailing yacht performance 
prediction is addressed in the final part of the paper. Although 
it is not yet possible to cover the total subject of speed 
prediction through CFD, programs are now available to ad- 
dress many parts of this problem. A description is given of 
some of the more suitable programs for this purpose. An 
outline is presented of how these CFD techniques will even- 
tually fulfil the quest for ever-better yacht performance predic- 
tion techniques. 


This paper is in part an attempt to show that when more 
exactness and physical reasoning is introduced into the analysis 
of sailing yacht performance, particularly with respect to 
mathematical modelling and interpretation of model towing 
tank data, the overall result is an increase in prediction ac- 
curacy. In various areas, such as in the analyses of model test 
data for the purpose of VPP development (such as for equitable 
racing and handicapping purposes) and for important projects 
requiring the highest possible accuracy (America’s Cup 
projects for example), approximations are often made that lead 
to unnecessary inaccuracies. The assumptions, for example, 
that the total hydrodynamic side force is a linear function of 
leeway angle, that the total resistance due to side force is a 
function of the square of that side force, that at large angles of 
heel the side force is zero for zero leeway angle, etc, all need 
to be studied from case to case. When this is done it will often 
be found that such assumptions are not valid and that a more 
meticulous approach, allowing for greater detail in specific 
situations, will lead to more accurate speed predictions. 


2. BASIC CONSIDERATIONS 


2.1 Equilibrium of Aerodynamic and Hydrodynamic For- 
ces 


When running before the wind the resultant force developed 
by the action of the wind on the sails of a sailing yacht is 
approximately in the direction of the heading of the yacht. In 
that situation a sailing yacht does not distinguish itself from 
other craft with respect to the attitude of the hull in the water 
and the hydrodynamic forces caused by the flow of water along 
the hull (usually called the canoe body) and the appendages. 
When beating to windward or when on a reaching course 
however the resultant force of the wind on the sails has a 
component at right angles to the yacht’s heading, as seen in 
plan view. This transverse force can be up to about five times 
greater than the forward thrust component of the sail force. 
Since in any steady sailing condition equilibrium must exist 
between the aerodynamic and hydrodynamic forces on the 
yacht, it follows that the transverse component of the wind 
force must be in equilibrium with the transverse component of 
the resultant hydrodynamic force on the canoe body, keel and 
rudder. This equilibrium of aerodynamic and hydrodynamic 
forces is shown in Figure 1. 


The only way for a yacht to develop the required hydro- 


dynamic transverse force, normally referred to as side force, is 
to adopt an angle of yaw 6 (called leeway angle in sailing 
yacht terminology) relative to the track of the yacht through 
the water. The keel and rudder of a yacht do not normally 
have camber, requiring that the total side force be generated by 
angle-of-attack only. Camber in the lifting appendages cannot 
be adopted because a sailing yacht must perform equally well 
on port and starboard tack (an aircraft does not have to per- 
form equally well when flying upside down). Rotatable keels 
and forward rudders as fitted to the latest generation of 
America’s Cup yachts, and leeboards as fitted on traditional 
Dutch sailing craft, allow for the development of side force 
without the necessity of leeway angle. The side force re- 
quirement in moderate to high wind speeds is such that ad- 
ditional side force usually needs to be generated by leeway 
angle in these cases as well however. 


To design a sailing yacht with good ability to sail to wind- 
ward, it is necessary that the required side force be generated 
at a small leeway angle. A large leeway angle seriously im- 
pairs the windward ability of a yacht. In part this is due to the 
increase in hydrodynamic resistance of the hull with increasing 
leeway angle. Mainly, however, a large leeway angle results in 
a lower value for the component of the boat speed in the 
direction opposite to the true wind speed V,,, than in the case 
of a small leeway angle, as can be seen in Figure 1. This 
component of the boat speed, called the speed-made-good to 
windward Vy, is equal to: 


Vac = Vp-Cos(Brw+8) 


The most important single problem in the design of a yacht, 
when windward performance is pre-eminent, can be formulated 
as how to design an underwater configuration that generates 
the required side force for a minimum of resistance at a small 
leeway angle or, how to maximize the side force-resistance, or 






Baw: Ea +€H-B 
Ys ~Ey - Ey 


Figure 1. Definition of forces and angles in the horizontal 
plane as used in mathematical relations for a yacht sailing at an 
angle of By to the true wind direction with velocity V, and 
leeway angle 6. 
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lift-drag ratio, for a given value of the side force, at a smail 
leeway angle. 


2.2 Hydrodynamic Side Force 


As the hull of a sailing yacht adopts a leeway angle to coun- 
teract the aerodynamic side force, it also adopts a heel angle. 
After equilibrium between all forces acting on the yacht is 
attained and a steady state value of the leeway angle 6 and heel 
angle ¢ is reached, the angle of incidence of the flow to the 
keel, which is normally the main part of the yacht with respect 
to the production of side force, is @cos@ (in fact 
arctan(tan@.cos@) to be precise), where @ is the angle of heel. 
If the canoe body is considered as a lifting body then its angle- 
of-attack to the flow is also Bcos@. In the case of the rudder 
however, when in the normal position behind the keel, the 
angle-of-attack is usually reduced because of the downwash of 
the flow behind the keel. A first estimate of this downwash 
effect is that it reduces the effective value of Bcos@ by half 
(Scott, 1974). As is described in Paragraph 4.3.10, the effec- 
tive angle of incidence of the flow to the rudder appears to be 
significantly more complex to determine, being not only depen- 
dent on keel geometry and keel lift but heel angle and Froude 
number as well. 


The horizontal side force component of the lift produced by 
the keel can be expressed as: 


SF = 4pVz7.(8C,x/da).Bcosd.Ayx-COS.Lep xéac 


where p = fluid density; 

Ve = boat speed; 

OC\x/da = lift-curve slope with respect to angle-of- 
attack, often assumed constant over the 
range of angle-of-attack of interest; 

Atx = lateral area of keel (equal to span times 


average chord length); 
= factor accounting for the effect of the 
canoe body on the lift of the keel. 


and Loca Ktact 


Similar expressions can be defined for the rudder, based on 
the appropriate angle-of-attack, and any other lifting appen- 
dages the yacht may possess. For the small leeway angles 
normally experienced (generally less than about 8 degrees), the 
lift-curve slope factor 0C,,/da can usually be considered as a 
constant for a given planform, if the aspect ratio is not too 
small. It follows that the side force on keel and rudder can be 
considered linear with respect to leeway angle. This cannot be 
assumed for the case of the canoe body as it normally con- 
tributes an important non-linear side force component, approxi- 
mately proportional to the square of the angle of incidence. 
This is described more fully in Paragraph 2.3.6.2. Because of 
this non-linear side force component, the graph of total side 
force against angle-of-attack Bcos¢ is usually not a straight line 
but a curve with increasing gradient, with increasing value of 
Bcos¢. 


The influence of the canoe body on the lift of the keel can be 
appreciable. Depending on the length-beam ratio, the beam- 
draft ratio, the deadrise angle of the canoe body, the heel 
angle, and the degree of wave-making exhibited by the yacht, 
an important increase in the side force on the keel is often 
found. By similar interaction, the keel enhances the side force 
on the canoe body. At relatively large angles of heel the inter- 
action phenomena involved result in less enhancement of the 
side force. This is discussed more fully in Paragraph 4.3.7. 


A typical example of the total side force, as measured on a 


model in a towing tank for different speeds, leeway angles and 
heel angles, is given in Figure 2. The model concerned is a 1:7 
scale model of an International America’s Cup Class (ACC) 
yacht, the results of which will be used to illustrate various 
matters discussed in this paper (see also Van Oossanen, 1992). 


2.3. Hydrodynamic Resistance 
2.3.1 Introduction 


The hydrodynamic resistance experienced by a sailing yacht, 
for the purpose of speed prediction, can be divided into six 
main components. Three of these are common to all types of 
vessels, viz: viscous resistance, wave resistance and added 
resistance in waves. Because of the significant hydrodynamic 
side force on the underwater body of a sailing yacht however, 
three additional resistance components are incurred. The first 
of these is a consequence of the significant angle of heel 
adopted by the hull, which results in variations to the viscous 





ZERO 


IDE FORCE 
ANGLE FOR HEEL 
ANGLE |( @ ) = 30° 














° 25 5 75 10 125 
6 COS 9 IN DEGREES 


Figure 2. Measured side force on a 1:7 scale model of an 
International America’s Cup Class (ACC) Yacht. Note the 
more than linear increase in side force with increasing angle- 
of-attack at heel angles of 10 and 20 degrees and, at 30 
degrees of heel, the decrease in side force at leeway angles in 
excess of about 7 degrees. The effect of Froude number on 
side force is here less significant than is usually the case, 
allowing the values for 4 different speeds to be plotted on one 
curve for constant heel angle. Also note that zero side force at 
30 degrees of heel occurs at Bcos# = 0.7 degrees, or B = 
0.81 degrees. 
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resistance, the wave resistance, and the added-resistance-in- 
waves component, as found for the zero-heel condition. This 
component is referred to as the resistance due to heel. The 
second component, called induced resistance, is associated with 
the induced flow (the upwash in front of a lifting surface and 
the downwash behind it) in that it causes the lift to possess a 
vector component acting opposite to the direction of motion. 
The third resistance component is the direct result of the 
attitude of the hull in terms of the leeway angle, keel rotation 
angle, or whatever variable utilized to produce the side force, 
in causing variations to the basic resistance components des- 
cribed above. In the usual case when side force is produced by 
the yacht in adopting leeway, this resistance component is 
mainly associated with the separation of flow from the canoe 
body. This component can be termed the resistance due to hull 
attitude. 





The main properties and features of these six resistance 
components will now be described. 


2.3.2 Viscous Resistance 


The viscous resistance is considered to be equivalent to the 
sum of the frictional resistance and a pressure resistance com- 
ponent of viscous origin. The frictional resistance is associated 
with the force required to overcome the tangential stresses 
developed between the canoe body, keel and rudder and the 
viscous fluid. The viscous pressure resistance is due to the 
pressure difference between the forebody and the aftbody as a 
consequence of the growth of the viscous boundary layer. This 
causes the pressure on the aft-body to be lower then on the 
fore-body, leading to a resultant pressure force on the hull. 
Since this pressure force has a viscous origin, as distinct from 
the forces associated with the pressure distribution over a body 
in potential flow, it is accordingly termed the viscous pressure 
resistance. 


Since the pressure resistance component is usually small, the 
viscous resistance is often expressed as: 


Ry = “4pV2.Cy.S 


where Cy = Cp(1+k) 

in which C, = specific frictional resistance coefficient of 
the body, here assumed equivalent to that 
of a flat plate with the same Reynolds num- 
ber as that of the subject body; 

1+k = three-dimensional form factor accounting 

for the additional resistance of the subject 
body over that of the equivalent flat plate; 

and S = wetted surface area of the body. 


It should be noted that in this expression the form factor not 
only accounts for the viscous pressure resistance but also for 
the difference between the frictional coefficient of the equiva- 
lent flat plate and that of the subject body. This difference is 
due to the fact that the wall shear stress in the boundary layer 
along a three-dimensional body is different to that along a flat 
plate. 


2.3.3. Wave Resistance 


At normal sailing speeds, for Froude number values around 
F,=0.35, the flow along the canoe body leads to an increase 
in pressure in the region of the bow, a decrease in pressure 
amidships and an increase in pressure at the stern. This leads 
to a wave peak at the bow, a wave trough amidships and a 


wave peak at the stern. The energy expended by the hull in 
causing these waves is significant and greatly dependent on the 
length-displacement ratio L¢,/Vc,'3. To illustrate this we can 
write, for the canoe body, at zero heel and leeway, in still 
water: 


Rice = Rvcp + Rw = 4eVe*.Cerce-Scp 


in which Cetce = Crcg(] +Keg) + Cy 


where Rycg = viscous resistance of canoe body; 
Ry = wave resistance of canoe body; 
Cercp = total hydrodynamic resistance coefficient of 
canoe body: 
Scp = wetted surface area of canoe body; 
and Cy = wave resistance coefficient. 


A useful formulation of the wave resistance coefficient is: 
Cy = a.F," 


in which a and n are constants. From Prohaska-type plots of 
the upright resistance (Prohaska, 1966), it can be seen that for 
many Sailing yacht models, the value of n is approximately 
equal to 4. On using n=4, when the coefficient a usually only 
varies to a small extent over the speed range up to F,=0.35, it 
follows that: 


Ry = 4apV—2Sep-Fn" 
i.e. Ry = YapVp° Scp/(g-Lepy” 


From studies of the wetted surface area Sc, of similar types 
of canoe bodies it is known that: 


Cs\(Vep-Lcp)” 


It 


Sce 
ae 2: 
or Sep = Cgo-Ven3, 


where C,, and C,, are coefficients. Note that the latter for- 
mulation can always be expressed as: 


i} 


Scp Cs(Vep-Lce)” 


since (Veg-L cp)!” Von? (Lee! Von!'3)"2 
It follows that: 


AapVg.Csi(V cp-Leg)'(g-Lcp)? 


iH 


Rw 


It 


i.e. Ry A(apCg,/g?)Vp°/(Lcp/V cp)? 


while, on using Scg = Cg.Vep"?, we obtain: 
Rw = A(apCg)/g”)Ve'Lce/Ven?)? 


It is therefore clear that when the length-displacement factor 
Lep/Vcp'? increases by the factor f, the wave resistance de- 
creases by at least the factor 1/f 3”, if not 1/f?. 


It is this component of the total resistance which, in a practi- 
cal sense, limits the attainable speed of mono-hull, ballasted 
sailing yachts with a low to moderate value of the length- 
displacement ratio (ower than about 5.5), to a value in knots 
between 2.35(L¢,)'” and 2.70(Lc,)'”, where Leg is the length 
of the canoe body on the waterline in meters. This "maximum" 
speed is often referred to as the “hull speed". It should be 
noted that such a speed does not physically exist since the 
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speed keeps on increasing with increasing thrust. 
2.3.4 Added Resistance in Waves 


The combined heave, pitch and roll motions experienced in 
head and bow-quartering seas and the roll motion experienced 
in beam seas, sometimes in combination with smaller ampli- 
tude surge, yaw and sway motions, are considered to be the 
main cause of the so-called added resistance in waves, some- 
times simply called the added resistance. In actual fact this 
added resistance is the sum of a number of different com- 
ponents, the first and most important of which is the resistance 
associated with unsteady hull motions in producing damping 
waves. Another component caused by unsteady hull motions is 
that due to the unsteady flow and flow separation along the 
hull, particularly if this unsteady flow is caused by roll 
motions. A third component of the added resistance is associa- 
ted with the phase shift between the waves encountered and the 
damping waves the hull causes. A fourth component is as- 
sociated with the reflection and refraction of waves from the 
hull. A fifth, usually very small component is caused by the 
unsteady rudder movement necessary to sail a straight course 
in waves. 
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Figure 3. Comparison of the resistance at zero heel and lee- 
way in calm water and in waves for the International 12 Meter 
Class yacht “Australia I". The curve shown for calm water is 
based on model test results. The curve for the added resistance 
in waves was calculated for a specific set of wave conditions 
measured at the America’s Cup buoy off Newport, Rhode 
Island, representative of the average conditions in the early 
afternoon during the month of September. 


The severity of pitching is particularly found to influence the 
level of the added resistance component in waves. On reducing 
pitch motions, the level of the added resistance is reduced. 
Another important feature of this resistance component is that 
it is independent of the still water resistance and that the added 
resistance in waves is proportional to the square of the wave 
height. 


The maximum added resistance usually occurs at or near 
pitch synchronism. If pitch synchronism occurs at low boat 
speeds the total resistance curve actually displays a local 
minimum such as shown in Figure 3. No steady-state sailing 
condition can be sustained in this speed region however, 
mainly due to the variation in total sail thrust with the non- 
steady motion of the yacht. 


In following seas the overall effect of the waves often results 
in a speed increase of the yacht. Particularly when the length 
of the waves is more than 5 times the length of the yacht, the 
length-displacement ratio of the yacht is more than about 6, 
and the following waves have a low encounter frequency, the 
yacht will surf down the face of the overtaking wave, momen- 
tarily reaching very high speeds. At high encounter frequencies 
and in less regular seas the effect of the following waves is 
generally still found to increase the resistance. 


2.3.5 Resistance due to Heel 


The resistance due to heel is defined as the sum of the 
changes to the viscous, wave and added resistance components, 
as caused by heel only. It is normally only determined for the 
still water condition and the effect of the change in added 
resistance with heel is in most studies (conveniently) assumed 
to be negligible, although measurements have revealed that at 
least for 12 Meter yachts this is significant (Klaka and 
Penrose, 1987). Since the resistance due to heel does not 
include the induced resistance, its only means of experimental 
determination is to extrapolate the curve representing the 
resistance increase found with increasing side force (due to 
leeway angle for example), for constant heel angle and forward 
speed, to the condition at which no side force is exerted on the 
model. This entails carrying out model tests at varying values 
of the transverse stability of the model since it is necessary to 
maintain the same heel angle while the leeway, i.e. the side 
force, varies. The alternative method of carrying out straight- 
forward towing tests at zero yaw angle while the heel angle is 
increased at constant speed, by re-positioning weights in the 
model, is not to be preferred above the method of varying the 
side force at constant heel angle because in general it cannot be 
assumed that the side force on the canoe body remains zero 
while the heel angle is varied in that case. 


The resistance due to heel is sometimes less than zero for 
heel angles less than about 20 degrees. This is usually due to a 
reduction in wetted surface area as the yacht heels, particularly 
when the canoe body has very little topsides flare. In that case 
the decrease in viscous resistance is often greater than the 
resistance increase associated with the lesser quality of the flow 
over the canoe body. 


2.3.6 Induced resistance 


2.3.6.1 Induced Resistance Associated with Linear Lift 


Lift on a body such as a keel or a rudder is generated by 
deflecting the incoming flow over an angle sideways. In terms 
of the better-known terminology relating to aircraft wings, this 
change in the flow relative to the lifting body is referred to as 
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an "upwash” in front and a "downwash" behind the body. The 
resultant force on the body acts approximately at right angles 
to the resultant flow velocity at the body. The lift or side force 
is defined to be the useful component of this resultant force 
acting at right angles to the undisturbed flow direction as 
shown in Figure 4. It follows that in terms of this lift force L 
the resultant force F,,, is equal to L/cosa; and the component in 
the direction of the undisturbed flow is R,=F,,,.sina, or 
R,=L.tana;. This last component is experienced as a resis- 
tance. It is termed the induced resistance because it is as- 
sociated with the deflected or induced flow angle a, which 
causes the resultant force to exhibit a component in the direc- 
tion of motion in the first place. 


It is a well known fact (see e.g. Kuethe and Chow, 1986, 
page 148 to 152) that the induced flow angle at the body a, is 


approximately equal to C,(1+0)/7AR,, where 
C.=L/(‘4pV".A,), so that for small angles a: 
R, = L.tana, = L.a; 
i.e. R, = 4pV?.A,.C,2(1 +.0)/(@AR,) 
in which R, = induced resistance; 
A, | = lateral area of body, normally expressed as 


span times average chord length; 

lift coefficient; 

correction factor for non-elliptical lift dis- 

tribution on lifting surface; 

and AR, effective aspect ratio (in many cases equal 
to two times the geometric aspect ratio). 


+ 
Q 
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It follows that the induced resistance is a linear function of 
the square of the total lift or the square of the effective angle- 
of-attack, when lift is a linear function of the angle-of-attack. It 
should be noted that the effect of heel only influences the 
induced resistance through the effect thereof on the aspect ratio 
AR, and, generally, on the value of 1+o. On utilizing the 
measured side force for the calculation of the induced resis- 
tance, this must be divided by the cosine of the heel angle first 
to find the appropriate value of the lift force or lift coefficient. 












* Cos.ai 
RESULTANT 
INFLOW 
Ry= F SIN. Gi 
UNDISTURBED FLOW DIRECTION i tans 


Figure 4. Induced resistance component of the resultant force 
F,,, on a lifting body as a consequence of the induced flow 
angle a;, as constituted by the upwash in front of the lifting 
surface and the downwash behind it. 


2.3.6.2 Induced Resistance Associated with Non-Linear Lift 


The resistance due to lift of bodies that possess an important 
non-linear lift term in the angle-of-attack (notably in the case 
of bodies with geometric aspect ratio values less than about 
0.5) is only in part calculated with the expression given above 
for the induced resistance associated with linear lift. The 
theory indicates that such bodies have a linear lift term in the 
angle-of-attack and a lift term approximately dependent on the 
square of the angle-of-attack (Hoerner, 1965, page 7-16 to 7- 
20). The linear term is considered to be associated with the in- 
duced flow as in planforms with greater aspect ratio, while the 
non-linear lift term is considered to be associated with the 
separated flow from the lateral edges of the body. 


The induced resistance associated with the non-linear lift 
component L,, is defined similarly to the induced resistance 
associated with linear lift: 


Ry = Lyp.tana; or Cpiago = Cyg2-tane; 


Here however C,,. is proportional to a? and not to a. Together 
with the fact that in small aspect ratios the total induced angle 
is approximately equal to the total angle of incidence (the flow 
becomes aligned with the body before the trailing edge is 
reached, with all of the lift being produced close to the leading 
edge), it follows that for small a, a, = a@/2, i.e.: 


R,o = “’eV?.Ay.Cy,9-0 
where the coefficient of non-linear lift is: 
Cree = Cig2-07 


in which c,,.) iS a constant for a given geometry, speed, heel 
and leeway angle. 


It follows that whenever the total lift on a body possesses a 
significant non-linear component, the total resistance due to lift 
can no longer be considered to be strictly proportional to the 
square of that lift or to the square of the angle-of-attack. For 
the leeway angles adopted by the hull of a well-designed 
sailing yacht when sailing to windward, the resistance due to 
this non-linear lift is small and often not noticeable in the data 
derived from towing tank tests. In that case this resistance 
component will only come into play during tacking and gybing 
manoeuvres, when the leeway angle can momentarily attain 
large values. When the side force produced by the appendages 
is too small or when the appendages are particularly inef- 
ficient, requiring large angles-of-attack to generate the required 
side force for equilibrium, this component of resistance will 
clearly play an important part in the overall performance of the 
yacht. 


The extent to which non-linear lift is present in model tank 
test data can be ascertained from the measured values of the 
total side force for a given speed and heel angle as is described 
in Paragraph 3.5.4. The determination of the corresponding 
Ria. resistance component, dependent on the cube of the angle 
of incidence of the flow, is described in Paragraph 3.5.6.2. 


From these considerations it follows that to reduce the magni- 
tude of the induced resistance it is extremely important to 
maximize the value of the effective aspect ratio AR,. Apart 
from design considerations involving novel means to increase 
the effective aspect ratio through, for example, the use of 
winglets, the aspect ratio of the keel is often limited by restric- 
tions to the maximum draft. For International 5.5, 6, 8 and 12 
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Meter class yachts for example the geometric aspect ratio of 
the keel, equal to the span divided by the average chord 
length, is typically no more than about 0.3 to 0.5 and the 
effective aspect ratio approximately double thereof. For such 
yachts the induced resistance component is appreciable. 


2.3.7 Remaining Resistance Caused by the Generation of 
Side Force: The Resistance Due to the Lift Attitude of 
the Hull 


It is important to observe that the resistance associated with 
the generation of side force, as can be found from model test 
results by subtracting the value of the resistance at zero heel 
(Ryy=0) from the total measured resistance at non-zero values 
of heel, for the same speed (R;), is generally not only com- 
posed of resistance due to heel and the induced resistance. At 
values of the leeway angle in excess of about 8 degrees or 
more, and with increasing beam-draft ratio, the flow along the 
canoe body will generally have a greater viscous resistance 
(i.e. a greater value of the form factor), while the wave resis- 
tance and the added resistance in waves will also be subject to 
change. Often, significant flow separation from the canoe body 
will start to occur as well. It is necessary therefore to also con- 
sider the effect on resistance of the attitude or orientation of 
the hull and of the (rotatable) appendages, required to generate 
lift, as distinct from the resistance that can be directly related 
to the level of the side force, viz: R, and R,,). 


In the same way that heel angle is generally associated with a 
change in resistance, so is leeway angle or the rotation angle of 
one or more of the appendages generally associated with a 
change in resistance. This change in resistance is often found 
to possess a pronounced coupling with heel angle. This com- 
ponent of resistance can be termed the resistance due to the 
“lift-generation attitude” of the hull or simply the "resistance 
due to attitude” of the hull, R,,,, not accounted for in the sum 
of the resistance due to heel and the induced resistance. 


2.4 Examples of the Resistance Caused Through the Gene- 
ration of Side Force, with Particular Reference to the 
Resistance Due to Attitude R,,, 


The additional resistance forces that come into play when the 
hull creates side force are defined above to be the resistance 
due to heel, the induced resistance, and a component associated 
with the lift attitude of the hull. As stated above, this last com- 
ponent is usually not related to the level of the side force 
generated but is directly associated with the orientation or 
attitude of the canoe body and appendages in producing the 
side force. Very often this resistance is associated with flow 
separation. The relative importance of these 3 resistance com- 
ponents will now be illustrated by 2 examples. 


2.4.1 R,,, Associated with Flow Separation 


Figure 5 shows the results of resistance measurements on the 
same model for which the side force measurements are given 
in Figure 2, for 30 degrees of heel and for 4 forward speeds. 
These resistance values are set out against the square of the 
effective angle-of-attack Gcos@, which representation results in 
a straight line for the induced resistance associated with linear 
side force. The resistance due to heel for the lowest and 
highest of the 4 speeds is indicated. For the Froude number 
value of 0.273, for which relatively large leeway angles were 
considered during the towing tank tests, a sharp increase in 
total resistance is observed. For the same leeway angle values 
a small, yet conspicuous decrease in side force was found to 
occur (see Figure 2). This fact points to the probability that 


considerable flow separation occurs along the canoe body at 
these leeway angles. Here the resistance due to heel (at zero 
side force) is relatively small compared to the induced resis- 
tance at nominal values of the leeway angle, such as are 
adopted by the hull when sailing to windward. The resistance 
associated with the attitude of the hull at larger leeway angles, 
such as could occur during a tacking manoeuvre, is here 
significant. 


2.4.2 R,,, Associated with Wave-Making 


Figure 6 shows the results of towing tank tests with another 
1:7 scale model of an IACC yacht for which the side force was 
generated by rotating the entire keel fin (not the bulb attached 
thereto). The model itself (including the large keel bulb) was at 
zero leeway, while the rudder was fixed at zero helm. In this 
arrangement, the keel is the only lift-producing part of the 
yacht, aside from some lift carry-over onto the canoe body. 
The angle-of-attack is 6,.cosd, where d, is the rotation angle 
of the keel fin. The values of all of the resistance components 
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Figure 5. Measured resistance on the same model for which 
the side force measurements are given in Figure 2, at 4 dif- 
ferent speeds and a constant heel angle of 30 degrees. Note the 
increase in resistance with increasing value of the square of the 
effective angle-of-attack for the lowest speed, due to separation 
of the flow from the canoe body. This resistance-due-to-at- 
titude component is not proportional to the value of the side 
force and is a consequence of the attitude of the hull only. The 
relative values of the other resistance components is indicated 
for the same speed (Fy = 0.273). 
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for the Froude number value of 0.270 are indicated. The value 
of the frictional resistance R; herein was determined by cal- 
culation but all other components follow from the test data 
analysis method discussed in Paragraph 3.5. 


The induced resistance varies linearly with the square of the 
angle-of-attack as could be expected. No appreciable canoe 
body lift (other than that associated with canoe body - keel 
interaction) can develop with the hull at zero leeway angle and, 
since the keel fin has a high effective aspect ratio, no non- 
linear lift exists. Accordingly, the straight line shown, deter- 
mined by defining the tangent to the curve through the 
measurement points at the location of zero side force, repre- 
sents the entire resistance related to the value of the side force. 


From Figure 6 it can be seen that at relatively large values of 
the rotation angle of the keel fin a significant increase in 
resistance occurs. This resistance increase is obviously as- 
sociated with some feature of the rotatable keel. Observations 
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Figure 6. Results of resistance measurements on a 1:7 scale 
model of an International America’s Cup Class (IACC) yacht, 
for a heel angle of 30 degrees, fitted with a rotatable keel fin. 
The model itself (including the large keel bulb) was at zero 
leeway, while the rudder was fixed at zero helm. The resis- 
tance due to attitude R,,, at large values of the side force, is 
here due to the gap between keel fin and canoe body which 
occurs as the fin is turned through too large a rotation angle 
6x. This was seen to give rise to a pronounced wave at the 
location of the keel. Here the form factor 1+kcg, to derive the 
value of Ry, was determined from the so-called Prohaska plot 
(Prohaska, 1966), as described in Paragraph 3.5.2. 


revealed a significant keel wave being formed as the gap 
between the top of the keel at the leading and trailing edges 
and the canoe body opened up, at relatively large values of the 
rotation angle. Accordingly, this resistance is at least in part 
associated with greater wave-making. It is possible that the 
occurrence of the gap between keel fin and canoe body also 
reduces the effective aspect ratio of the keel, in which case 
part of the resistance increase will be due to an increase in the 
induced resistance. The results of the side force measurements 
do not indicate a significant decrease in lift however, which 
would have to be the case if the effective aspect ratio were to 
decrease. 


It is important to note that when the canoe body plays a role 
in the development of side force, such as is normally the case 
when the hull adopts leeway, the slope of the curve depicting 
the induced resistance will generally vary with speed. This can 
be clearly seen in Figure 5, in which the slope of the lines for 
Rymy increase with increasing Froude number. The change in 
the interaction between keel and canoe body and the differen- 
ces in the wave pattern between the windward and leeward 
sides of the yacht, at varying speeds, leeway angles and heel 
angles, causes the effectiveness of the canoe body in deve- 
loping side force to differ from case to case. At zero leeway 
the canoe body plays a lesser role in this respect, causing the 
slope of the lines depicting the induced resistance to be almost 
independent of speed, as follows from Figure 6. 


2.5 Aerodynamic Forces 
2.5.1 Components of the Total Aerodynamic Force 


The aerodynamic forces experienced by a sailing yacht are 
those associated with the wind acting on the sails, hull, super- 
structure, mast, standing and running rigging, crew, etc. In 
general, these are all friction, viscous pressure and flow sepa- 
ration resistance forces except for the lift forces acting on the 
sails. It should be noted that in the field of aerodynamics, it is 
common to refer to resistance forces as drag forces and to 
refer to viscous resistance as parasitic drag. 


On windward and reaching points of sail, the induced resis- 
tance (induced drag) of the sails is by far the most significant 
of the drag forces. This is a consequence of the relatively very 
high lift coefficient the sails then have. Sailing downwind, the 
sails of common displacement-type yachts are in what is 
termed the “drag mode" and act similarly to that of a para- 
chute. In this situation the sails usually experience a relatively 
small lift force but a near-maximum attainable drag force. 


2.5.2 Lift and Drag Forces on the Sails 


The determination of the aerodynamic forces on the sails is 
usually much more difficult to achieve then the determination 
of the hydrodynamic forces on, for example, the keel. This is a 
consequence of the fact that no accurate experimental technique 
has yet been devised for the experimental determination of the 
lift and drag components of the force on the sails other than 
carrying out very expensive large scale tests on the actual 
yacht or on special test rigs. There are three prime reasons for 
this. The first of these is associated with the fact that it is 
extremely difficult to model the stretch and other properties of 
the sail material on small scale. The second is associated with 
the fact that when the sails are set to their optimum trim, yiel- 
ding the highest boat speed, the sails are often found to operate 
close to their maximum lift coefficient on reaching points of 
sail and at their maximum drag coefficient on downwind points 
of sail. Such sail settings invariably invoke a significant 


Predicting the Speed of Sailing Yachts 347 


amount of ‘separation of flow from the sails. When, as is 
normally the case, the full-scale Reynolds number cannot be 
adhered to in the test facility (usually a wind tunnel), this fact 
leads to significant inaccuracies since flow separation phe- 
nomena are prone to large scale effects when the Reynolds 
number is not adhered to. 


Thirdly, the significant interaction between the various sails 
is dependent on many variables, such as the amount of sail 
overlap (between genoa and mainsail and between mainsail and 
spinnaker), sheeting angles, sail shape (as associated with 
aspect ratio and the amount of roach or leech hollow), sail 
trim, etc, which makes it almost impossible to accurately 
predict the sail forces at full scale from previously-attained 
model or full-scale test results, irrespective of any other 
problem. This alone is significant enough to warrant, in 
America’s Cup competition for example, special video camera 
arrangements on deck and in the mast to ensure that previously 
recorded sail shapes and sail trim situations at which good boat 
speeds were logged can be reproduced. 


Up to now the problem of accurately predicting sail lift and 
drag forces has still not been dealt with adequately. Many 
noteworthy experimental studies have been carried out (see 
e.g. Marchai, 1988) which have helped to provide insight into 
this difficult subject and which have yielded important quan- 
titative information. A notable project in this respect, utilizing 
a 40% model of an IACC yacht, was recently reported on 
(Milgram et al, 1993). Also, many numerical studies have been 
carried out to derive design data for sail lift and induced drag 
(see e.g. Milgram, 1968, and Euerle and Greeley, 1993). None 
of these studies have been conclusive, however, with respect to 
the lift and drag data obtained. 


Notwithstanding this, a situation has presently been reached 
in which reasonably accurate sail lift and drag data have been 
defined for the purpose of determining maximum boat speed 
potential. This data represents the upper bound of attainable 
sail performance, assuming optimum sail trim and sheeting 
angles and, hence, cannot be used to determine the actual 
performance on the water at any one time, at non-optimum sail 
settings. This data assumes that the yacht is at all times sailing 
to its full potential. 


A representative set of lift and drag coefficients in this 
respect is that defined in the International Measurement System 
for use in predicting speed potential for rating and handicap- 
ping purposes (Poor, 1986). These values are given and dis- 
cussed more fully later in Paragraph 4.9. Such values, when 
multiplied by the appropriate sail areas and the dynamic pres- 
sure of the apparent wind, yield boat speed values which have 
been found to correlate reasonably well to full-scale values 
whenever accurate full-scale measurements are taken such as 
during research and development programs carried out for the 
America’s Cup. The IMS system defines correction factors for 
sail blanketing and overlap effects, for the effect of the ex- 
posed portion of the mast of fractional rigs, and the effect of 
using full-length battens. No effect of the aspect-ratio of the 
individual sails are included which, in the case of the lift 
coefficient, is a particular problem that needs to be addressed 
when implementing these sail coefficients for studies involving 
different rig heights. 


2.5.3, Aerodynamic Parasitic Drag Forces 
In most speed prediction methods, whether these be based on 


model tests or on numerical procedures, the calculation of the 
aerodynamic resistance of the hull, superstructure, mast, 


rigging and crew is usually not carried out in a very refined 
way. At best, a constant drag coefficient is adopted and the 
representative frontal area calculated. The reason for this is 
that the forces involved are normally relatively small compared 
to the forces exerted by the wind on the sails. When sailing 
high on-the-wind, however, close to the maximum attainable 
Vac, this drag force has an appreciable effect on performance. 
Accordingly, there is a reward in considering this force in 
some detail as well. This too will be discussed in more detail, 
in Paragraph 4.9. 


2.6 Hydrodynamic and Aerodynamic Force Inter-Depen- 
dence 


It is important to note that the resultant hydrodynamic and 
aerodynamic forces acting on a sailing yacht are inter-depen- 
dent. When the force of the wind on the sails is increased and 
a straight course is maintained (normally requiring a helm 
angle adjustment), the yacht will adopt a greater leeway angle 
thereby producing the necessary increase in the hydrodynamic 
lift force to maintain equilibrium. In actual fact this process is 
much more complex because the heel angle will increase at the 
same time when the sail force increases, causing initial loss in 
side force from the keel, in itself requiring a helm angle cor- 
rection and a greater leeway angle to maintain equilibrium. 


A study of Figure 1 reveals that for equilibrium it is neces- 
sary that the following equation is satisfied: 


Baw + B = ey + €q 


where Gay = angle of the apparent wind relative to the 
centerline of the yacht in the horizontal 
plane; : 
B = leeway angle; 
€y = hydrodynamic drag angle in horizontal plane; 
and Eq = aerodynamic drag angle in horizontal plane. 


The hydrodynamic and aerodynamic drag angles are defined as 
follows: 


€y = arctan(Ryy/SFy,) 
and EA = arctan(R;,/SF,,) 
where Ry, = total hydrodynamic resistance; 
SFr, = total hydrodynamic side force; 
Rr, = component of total aerodynamic resistance in 
horizontal plane; 
SF;, = component of total aerodynamic side force in 


horizontal plane. 


This equation is attributed to F. W. Lanchester (Marchaj, 
1988, page 12). It is utilized in most mathematical models in 
one way or another for expressing the inter-dependence of the 
hydrodynamic and aerodynamic side forces and resistance. 


3. HYDRODYNAMIC MODEL TESTING AS A BASIS 
FOR SPEED PREDICTION 


3.1 Introduction 


It is interesting to note that until 1983 the use of the towing 
tank for the speed prediction of sailing yachts was considered 
by some as a liability. In many America’s Cup projects the 
towing tank had failed in producing the winning design. Prior 
to 1936 this can be attributed to the lack of understanding of 
the nature of lateral resistance, for which tests at heel and 
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leeway need to be carried out. Although Hiley (Allan et al, 
1957) and Kempf (Barkla, 1960) have been attributed the 
honour of being the first to test yacht models at combined 
values of heel and Jeeway, yacht models were basically only 
tested at zero leeway until 1932. It was Davidson (1936) who 
showed for the first time that the lateral resistance is in fact as- 
sociated with the lift of the hull and that the increase in resis- 
tance caused thereby is associated with this same lift force in 
the form of induced resistance. 


When Davidson started his investigations in 1932 he found 
that for the study of on-the-wind performance it was necessary 
to tow a model at specific combinations of heel and leeway. 
Under his direction the first lateral dynamometer test set-up 
was designed and built. At the same time he found that rela- 
tively small models were very prone to maintain a laminar 
boundary layer and that the boundary layer needed to be 
tripped into turbulence to reduce the effects of model scale in 
predicting full-scale behaviour to acceptable proportions. 
Davidson’s 1936 SNAME paper, in which these two elements 
are presented and discussed for the first time, still makes com- 
pelling reading. 


Between 1936 and 1983, the failure of the towing tank to pin- 
point the winning or best design in a number of cases is more 
difficult to explain. Disbelief in the towing tank was mainly 
found in the United States. In part, this can be traced back to 
projects in which erroneous interpretations were supposedly 
made of towing tank test data such as for "Mariner" (Vaughan, 
1975). 


The development of "Australia II", first and foremost by use 
of the towing tank but also by velocity prediction programs 
(VPP) and computational fluid dynamics (CFD) programs (Van 
Oossanen, 1985), has shown that there is as yet no substitute 
for the towing tank in accurate speed prediction work. 


3.2 Testing Techniques 


Over the years many different test techniques have been 
investigated. Today, only two practical methods are in every- 
day use that yield accurate results for speed prediction at all 
points of sail. The first of these is the method devised by 
Davidson at the Stevens Institute of Technology, already 
referred to. In this method the model is restrained in the 
athwartships position at a pre-determined leeway angle by 
means of arms fixed to the model fore and aft. These arms are 
constructed in such a way that no vertical force is exerted on 
the model. This allows the model to freely rise, sink and trim 
with speed. The horizontal component of the lift exerted on the 
model at speed is measured at each of these arms where they 
are fixed to the model, at equal heights above the water surface 
at rest, by means of block strain-gage dynamometers. The 
arms are fixed to the model through these block gages in a way 
allowing the model to adopt a heel angle while the arms re- 
main vertical. The required angle of heel at speed is usually 
realized by moving a weight transversely along the deck before 
the run commences, which position is then often fine-tuned to 
compensate for the change in transverse stability of the model 
with speed. The resistance of the model, equal to the towing 
force, is measured in the towing connection between carriage 
and model, which connection is also able to accommodate any 
vertical movement of the model. Ballast weights need to be 
moved longitudinally between each run to compensate for the 
fact that the model is not towed from the center of the resultant 
aerodynamic force. The compensating trimming moment is 
directly proportional to the resistance of the model. Accor- 
ding-ly, an estimate of the resistance during each run needs to 


be made beforehand. A photograph of this set-up, as utilized at 
the Delft University of Technology towing tank, is shown in 
Figure 7. Various modifications of this relatively simple testing 
method are in use at different laboratories (see for example 
Herreshoff, 1964, and Kirkman, 1979). 


The second method adopted in various towing tanks is that 
originally defined (but never tested) at the National Physical 
Laboratory in England. in 1953 (Allan et al, 1957). This 
method, aimed at reducing the number of test runs required to 
analyze the performance of one configuration, is based on 
towing the model from the (estimated) effective center of effort 
of the resultant aerodynamic force on the yacht by means of a 
special tow-mast set-up. Here too, the test rig requires the yaw 
angle to be set at the beginning of each run. While the model 
is accelerated, side force is generated leading to a force and a 
torsional moment in the tow mast. The corresponding strain 
gage dynamometer signals are used to drive an electric motor 
to lower the vertical sub-carriage mounted on the main tow- 
tank carriage attached to the tow mast, thereby allowing the 
model to adopt that heel angle at which the axial force in the 
tow mast is zero. At the same time another electric motor 
drives the tow mast to that longitudinal position in the model at 
which the torque in the tow mast is zero. The equilibrium 
Situation is normally attained between 10 to 20 seconds. 


In the Davidson test set-up a matrix of heel and leeway 
values is run at each speed of interest. Here the choice of heel 
angle for every leeway angle is arbitrary and the transverse 
stability of the model is constantly adjusted. In the so-called 
free-running test set-up the transverse stability of the model is 
according to that of the full-scale yacht. In towing from the 
center of the aerodynamic force, with no other attachment to 
the carriage whatsoever, the full-scale heel angle will be 
obtained for every leeway angle adopted. 


Many initial problems had to be overcome before this free- 
running method yielded reliable and accurate results. This 
system was initially used by the Fluid Dynamics Laboratory of 
the Saunders Roe Division of Westland Aircraft Ltd, on the 
Isle of Wight (Crago, 1962). The National Research Council of 
Canada (Murdey, 1978) and the Netherlands Ship Model 
Basin, now called the Maritime Research Institute Netherlands 





Figure 7. Photograph of test set-up at the Hydrodynamics 
Laboratory of the Delft University of Technology. This test 
set-up is representative of the original Davidson technique. 
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(Gommers and Van Oossanen, 1983) perfected it. At the Mari- 
time Research Institute Netherlands (MARIN) the remaining 
problems with this method, mainly involving lack of equi- 
librium for some rudder settings, trim tab settings and tow- 
mast positions, were resolved by using the servo motor system 
described above to drive the tow-mast set-up in the model to 
the position at which the moment of all forces in the horizontal 
plane is zero. A photograph of this test set-up is given in 
Figure 8. 


3.3 Model Size and Turbulence Stimulation Requirements 


While significant studies have been carried out over the years 
concerning the effects of model length and turbulence stimu- 
lation techniques on the accuracy of test results, no definitive 
conclusions have been reached concerning the minimum re- 
quirements. The laboratories with Jarge towing tanks employ- 
ing large models maintain that the models need to be large, 
while the smaller towing tanks have repeatedly found that 
when paying particular attention to turbulence stimulation of 
the boundary layer on canoe body and appendages, the results 
obtained are accurate enough for most purposes, as initially 
demonstrated by Davidson himself. 


Noteworthy studies on the accuracy of yacht model tests have 
been carried out. By comparing the results of a series of 5 dif- 
ferent scale (geosim) models of the 5.5 meter yacht "Antiope" 
(with scale ratio’s of 1:6, 1:5, 1:3, 1:1.85 and 1:1 for the full- 
scale yacht) Kirkman and Pedrick (1974) showed that unaccep- 
table differences between the results for the different models, 
when extrapolated to full scale, start to occur when the water- 
line length of the smallest model becomes less than about 2.5 
meters. Similar comparisons for other yachts, for which 3 or 
more models of different size were tested, showed similar 
results. As the model size becomes smaller, the differences 
between the extrapolated test results for the smallest and the 
largest of the models become significantly greater. 


ne 





Research Institute Netherlands. With this test set-up, originally 
defined by the National Physical Laboratory, every run yields 
a heel-leeway combination in accordance with the transverse 
stability of the model. 


This is not surprising on considering the fact that a 
measurement error of 0.1 Newton in model resistance or side 
force, for example, represents a full-scale error of about 173 
Newton when the model scale ratio is 1:12 (0.1 x 12°) and 
only 2.7 Newton for a scale ratio of 1:3 (0.1 x 3%). This 
demonstrates that the demand to measure accurately increases 
dramatically as the model becomes smaller. 


One of the most serious problems incurred when adopting a 
small model is that associated with the error in extrapolated 
resistance when, by necessity, Froude scaling has to be adop- 
ted because of lack of knowledge about the appropriate 
Reynolds law, or when it is not known what part of the meas- 
ured model resistance is Reynolds dependent. This situation 
arises when the (common) practice is adhered to of scaling the 
resistance due to heel or the resistance associated with the at- 
titude of the hull by Froude’s law of similitude, when in actual 
fact Reynolds scaling needs to be adopted. The following 
example will illustrate this. 


Assume that in a particular case the resistance associated with 
heel is entirely related to a change in the viscous resistance of 
the canoe body, defined by: 


Ry = %pVp2-Sca-Crg-Crca(l +kep) 


with the coefficient of resistance due to heel determined from 
the results of model tests according to: 


Crom = Rooy/('Ze ey Vaemy Scaemy-Cecaenyl +kepomy)) 


If the scale ratio is \, then with Cp, = Cpga, Ve’ = A-Veuny> 
Scp = Scam, and Keg = Keg), We obtain: 


Ry = “pd. Vgmy’)-?.Scne)-Cromy-Cres(] + Keay) 


On adopting the 1957 ITTC flat plate friction coefficient 
correlation formulation for turbulent flow (see Van Manen and 
Van Oossanen, 1988, page 10 to 15, for an overview of the 
different friction formulations) we obtain: 


Crcam = 0.075/(Log(Ryceap) - 2)? 
and Crcp = 0.075/(Log(Rycp) - 2) 


Since Rycp = Vy-Lep/v, it follows that when differences be- 
tween the viscosity of sea water and fresh water is disregarded: 


Ryce = (A)? Vimy) O-Lepay)/? 
and Ry = Ryay-A2(Log(Rycam)) - 2)?/(Log(Rncaay-*”) - 2)? 


If Froude scaling is adopted, again disregarding differences 
associated with fluid density between model and full scale, we 
obtain: 


Ry = Re-*" 
The value of the expression: 
(Log(Rycam) ~ 2)/ (Log(Rycaa.*”) - 2), 


which constitutes the difference between the two scaling laws, 
is given in Table 1 for various values of the Reynolds number 
and the scale ratio. From this table it follows that as the scale 
ratio becomes larger, i.e. the model becomes smaller, Froude 
scaling of viscous resistance Jeads to too great a full-scale 
value of the resistance. For Rycpa) = 500,000 and A=12 for 
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example, the full-scale resistance is over-predicted by a factor 
equal to 1/0.484 = 2.07, while for a Reynolds number of 
5,000,000 and A=3 this factor is 1/0.753 = 1.33. If the full- 
scale value of the resistance due to heel is, for example, 5% of 
the total resistance, the over-prediction of the total resistance is 
then 5.35% against 1.65%, respectively. 


Table 1. Ratio of the value of the ITTC flat plate friction 
resistance coefficient on the full scale to that on model scale 
for different values of the Reynolds number and the scale ratio. 





Rycam) (Log(Rycam)) - 2)°/(Log(Rncamy-\*”) - 29 
for \=1 for \=3 for \=7 for \=12 

100,000 1.000 0.650 0.494 0.422 

$00,000 1.000 0.702 0.555 0.484 

1,000,000 1.000 0.720 0.577 0.507 

5,000,000 1.000 0.753 0.620 0.553 


It is the author’s experience that large differences between 
test results for different geosim models such as reported by 
Kirkman and Pedrick are not necessarily all a consequence of 
physical scale effects. Human "error" can also be the cause of 
as much difference, even between models of significant size, 
due to lack of duplication of the test set-up between one model 
and the next, lack of model manufacturing accuracy, errors in 
the adjustment of ballast weights and transverse stability char- 
acteristics, errors in dynamometer alignment and dynamometer 
calibration, etc. This can even occur when repeating ex- 
periments with one and the same model in the same laboratory, 
when all relevant model preparations are repeated. 


Present experience indicates that, on ruling out human error, 
minimum model size - when basing this on the requirement of 
avoiding significant physical scale effects - is dependent on the 
chord length of appendages such as keels and rudders rather 
than on the waterline length of the model. A minimum 
Reynolds number must be adhered to at which a trip wire, a 
row of studs, or distributed roughness such as sand or car- 
borundum particles, can still trip the boundary layer into 
turbulence. Various studies carried out reveal that boundary 
layer tripping on streamlined lifting surfaces at Reynolds 
numbers less than about 100,000 to 150,000 does not always 
lead to the desired result (Van den Bosch and Pinkster, 1971). 
Laminar flow and subsequent laminar boundary layer sepa- 
ration will then sometimes still occur, particularly in a favour- 
able pressure gradient. In those cases both the lift and the drag 
of the lifting surfaces are affected. 


The minimum requirement of a Reynolds number value of 
100,000 based on the average chord length of the smailest 
appendage (the rudder for example) leads to: 

Vaan (Ce/A)/y = 108 


in which Vain minimum speed of model during testing; 


CR = average chord length of full-scale rudder; 
mn = required scale factor for smallest model; 
and v = kinematic viscosity of water at prevailing 


temperature. 


The minimum speed of interest during model testing depends 


on the nature of the project being carried out. For current 
projects involving America’s Cup races off San Diego, where 
light wind conditions generally prevail, this minimum speed of 


- interest corresponds to a Froude number value of around 0.2. 


In this case it follows that: 
Vuan/(gZ-Lep/d)'? = 0.2 
or Vain = 0.626(L¢p/A)!” 
The model scale factor is then calculated from: 
h = (0.626Leg!”.cp/(105))*4 


For fresh water at 15°C (for which y=1.142 x 10° m*/sec), 
with Lep=19 m and cp=0.75 m for example, a value for 
=6.85 is calculated, corresponding with a length of the canoe 
body model, on the waterline, equal to 2.77 m. 


These criteria serve as guidance only. For studies not invol- 
ving racing yachts, it is sometimes acceptable to apply the 
above minimum Reynolds number rule to the keel and not to 
the rudder, depending on whether or not the tests to be carried 
out are to address helm balance and associated problems invol- 
ving rudder angle as a variable. 


Model test accuracy will always be served by adopting the 
largest possible model that can be handled in the towing tank 
without incurring problems associated with tank blockage. It is 
always advisable to discuss the whole matter of model scale at 
length with the experts at the laboratory involved. Herreshoff 
(1964), Kirkman (1974), Campbell and Claughton (1987), and 
DeBord and Teeters (1990) provide useful information on test 
accuracy. 
turbulence 


Comparable tests carried out with different 


‘stimulation techniques reveal that when heel angle is a test 


variable it is best to adopt at least 2 rows of studs. or 2 strips 
of distributed roughness (sand or carborundum) in a transverse 
plane along the girth on the canoe body. This is to ensure that 
all of the flow along the model is tripped into turbulence and 
not only the flow over the forebody. The advantage of using 
studs is that the parasitic resistance thereof can be approx- 
imately determined by calculation, requiring no further testing 
to correct resistance measurement results. When using strips of 
distributed roughness it is necessary to carry out resistance 
measurements with at least 2 different strip widths in order to 
be able to determine the viscous resistance (including the 
separation resistance) thereof. This is normally carried out by 
extrapolating the results of these measurements to the case of 
zero strip width. 


The position of the forward row of studs or strip of rough- 
ness on the canoe body, for it to be effective, should be such 
that the local flow Reynolds number is at least 6 x 10*. At a 
model speed of 1 m/sec for example, this amounts to a dis- 
tance from the stem (or the origin of the boundary layer) of 
about 7 cm. On appendages, where this requirement would 
lead to. the positioning of the studs or roughness aft of the 
location of maximum thickness (minimum pressure), it is satis- 
factory to position the studs or roughness at that maximum 
thickness location. At zero angle-of-attack, natural transition of 
the boundary layer on most types of (smooth) keels and rud- 
ders will occur at that location on the full scale between 
Reynolds number values of about 5 x 10° and 5 x 10°. 


The use of distributed roughness with a height k, such that 
the value of the Reynolds number based on k is at least 120 


Predicting the Speed of Sailing Yachts 351 


(i.e. Vaw-k/y > 120) has been found to give satisfactory 
results in both favourable and adverse pressure gradients (see 
Figure 17.45 of Schlichting, 1979). A disadvantage of using 
sand or carborundum roughness however is the need, as stated 
above, to carry out additional resistance tests to determine the 
parasitic resistance thereof. It has been found that on perfor- 
ming such tests for the canoe body alone and for the yacht with 
keel and rudder fitted, different results are obtained, which 
indicate that the resistance of the roughness on keel and rudder 
is appreciably greater then that on the canoe body. 


3.4 Matrix of Variable Values to be Investigated 


A model test program for the determination of the full-scale 
performance at all possible points of sail needs to address 
adequate ranges in model speed, heel and leeway angle. In 
addition, it is advisable to carry out additional runs with the 
model at 2 non-zero values of the rudder angle to obtain 
insight into how the rudder will augment total side force and 
how the location of the hydrodynamic center of side force is 
affected, particularly if the yacht, on windward points of sail, 
requires an optimum helm angle from a performance point of 
view. 


The model test program should preferably commence by 
carrying out the upright resistance test with the model of the 
canoe body without any appendages fitted and then with the 
model with all appendages fitted. Runs should be carried out 
for every 0.025 increment in Froude number, starting with the 
value Fy=0.1. This low value of the speed is required to 
discern. at what speed the flow can be considered to become 
fully turbulent (see Paragraph 3.5). The maximum speed to be 
selected is dependent on the length-displacement ratio of the 
model. For L¢,/Vc,'? less than about 5 (for heavy displace- 
ment yachts) the maximum speed rarely exceeds Fy = 0.5. 
When the length-displacement ratio is around 7, it is advisable 
to adopt a maximum Froude number value of about 0.75. In- 
between these 2 values of the length-displacement ratio a 
corresponding maximum Froude number value in-between 0.5 
and 0.75 or greater should be adopted. 


The bow-down trimming moment caused by the high position 
of the center of effort of the sail plan should be accounted for 
during these tests. If the model is to be towed from a lower 
height it is necessary to shift ballast weights between each run 
to compensate. To do this, it is necessary that the upright 
resistance curve be estimated beforehand. Rise or sinkage of 
the model at the center of gravity, and the trim fore and aft, 
should be measured in addition to the total resistance. It is 
essential that accurate upright resistance values are obtained for 
all speeds at which runs at non-zero values of heel and leeway 
are to be carried out. 


The tests at heel and leeway need to be performed for speeds 
between the Froude number values of 0.25 and 0.40. This 
speed range needs to be accounted for with a minimum of 4 
speeds. Since the Froude number value of about 0.35 repre- 
sents the speed at which the yacht usually attains its best 
performance (highest speed) to windward, it is essential that 
one speed is chosen close to this Froude number. At each 
speed a matrix of 4 heel angles and 4 leeway angles need to be 
investigated. When the transverse stability of the model is a 
variable, such as when using the Davidson test method (see 
Paragraph 3.2), it is recommended to adopt standard heel angle 
values of 0, 10, 20 and 30 degrees and to perform at least 4 
runs at different leeway angles for each heel angle, covering as 
large a variation in leeway angle as the limits in transverse 
Stability will allow. When the transverse stability of the model 


cannot be varied, it is advisable to carry out at least 5 runs per 
speed, each at a greater leeway angle then before, until a heel 
angle of 30 degrees or more is attained. 


The minimum number of runs necessary when employing the 
Davidson type test method is about 64 (4 speeds x 4 leeway 
angles x 4 heel angles) while when using the free-running test 
method, with modelled full-scale transverse stability, this will 
amount to about 20 runs. In both cases however it is usually 
beneficial to adopt more than 4 speeds, e.g. speeds correspon- 
ding to Froude number values of 0.25, 0.3, 0.325, 0.35, 0.375 
and 0.4. To analyze many of the effects of forward speed on 
side force and resistance (which effects are strongly non- 
linear), such speed data is invaluable (see Paragraphs 3.5 
through 3.6). 


3.5 Analysis and Extrapolation of Test Results to Full Scale 


The extrapolation of the test results to full scale should com- 
mence with the analysis of the upright resistance tests. In the 
first instance these results should be subjected to a study of the 
possibility that unwanted laminar flow has influenced the 
results, followed by a determination of the form factor of the 
canoe body (the value of 1+kc,). This is best carried out by 
the resistance measurement analysis technique proposed by 
Prohaska (1966). One of two methods can be utilized, depen- 
ding on whether or not a separate resistance test was carried 
out with the canoe body only, without appendages fitted. In 
extensive model testing programs, such as is usually carried 
out for America’s Cup projects, this should. be viewed as an 
essential requirement. ; 


3.5.1 Upright Resistance Test Without Appendages 


The specific resistance of the canoe body model is first deter- 
mined, equal to: - *¢ 


Crrcaim) = Rycaewy!/(AP cm) Vac Scam) 


where Rycg) = measured (total) upright resistance of canoe 
body model. 


Following Prohaska’s principle, the upright resistance of the 
model is divided into a viscous part dependent on Reynolds 
number and a part dependent on Froude number, viz: 


Crrcamy = Crcaa(] +Kcp) + Cw 


where Cy = a.F," 

The value of Cerca divided by the frictional resistance coef- 
ficient of the canoe body Crcpyy, here considered equal to the 
1957 ITTC formulation, viz: 


Ceca = 0.075/(log(Rncemy) - 2) 


in which Rycam = Vau)-Lcaay/¥amy 

is then plotted against Fy*/Cgcpy. In calculating Crop the 
actual length of the waterline of the canoe body Le, should be 
taken and not 0.7L¢, as originally proposed by Davidson 
(1936). No account need here be taken for parts of the hull 
with an appreciably shorter waterline length. An example of 
this analysis is given in Appendix 1, for the 12 Meter yacht 
“Australia II" without keel and rudder fitted. 


Laminar flow is present on the model if the value for 
Crrcao/Crcam for a specific speed is below the (straight) line 
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that fits the data points best for Fy*/Crcgay) values less then 
about 3 (see also Paragraph 3.5.3). The form factor of the 
canoe body 1+kcg is the value for Cyrcgim/Crcaqay On this line 
for Fy, = 0. 


With the value of the form factor thus found the viscous 


resistance of the canoe body in the upright condition for the 
model follows from: 


Rycam) = pV Cerca +kep)Scam) 


The wave resistance of the canoe body model is then found 
from: 


Rw = Rresm) - Rvcaa 


from which the wave resistance coefficient is determined from: 


Cy = Rwem/(A mV pany’ Scam) 
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Figure 9. Prohaska plots to detect the presence of laminar flow 
on the model and to determine the form factor of the canoe 
body, for the case of a 1:3 model of “Australia II" fitted with 
keel, winglets and rudder. The plot indicated by Ceyny/Crcam 
is obtained when the viscous resistance of keel, winglets and 
rudder is subtracted from the measured value to obtain a quasi- 
measured resistance for the canoe body at low speeds. The plot 
indicated by Cycaan/Crcam, derived from the test with the 
canoe body only, is included for comparison. The inset figure, 
at a smaller scale, shows the difference in resistance between 
the two plots, indicating the wave resistance associated with 
the appendages. 


which value, for a given speed, is considered constant irrespec- 
tive of the scale. 


The full-scale resistance is then determined from: 
Races = Rucn + Rw = 20Ve"-Sco(Crce(] +kep) + Cw) 
3.5.2 Upright Resistance Test With Appendages 


When the results of an upright resistance test with the canoe 
body only are not available, the analysis and extrapolation 
method described above needs to be performed for the upright 
resistance test for the model with all appendages fitted. This is 
best executed by calculating approximate values for the viscous 
resistance of keel and rudder models (and any other ap- 
pendages) and subtracting these from the total measured resis- 
tance. In that way quasi-measured upright resistance values are 
obtained which at low speed can be considered to be associated 
with the canoe body alone. At higher speeds this is not the 
case since the keel, rudder and other appendages, wherever 
positioned, normally have a wave resistance component. The 
procedure to be adopted in this case will now be described. 


The viscous resistance of a planar, keel-like appendage is ap- 
proximately as follows: 


Rvx = YapVp_?.Cey(1 +ky)Sx 


where 1+k, = 1 + 1.2(/c)y + 70(t/c),4 


which relation is adopted from Hoerner (1965). Here: 
(Ue), = average thickness-chord ratio of keel. 


When the flow over the keel is assumed turbulent, the coef- 
ficient of friction can be determined from the 1957 ITTC 
friction line. A similar formula holds for the viscous resistance 
of the rudder. In the case winglets or a bulb, or both, are fitted 
to the bottom of the keel, equivalent approximate viscous 
resistance formulations for these appendages have to be used. 


Non-dimensional values for the quasi-measured resistance of 
the canoe body (including the wave resistance component of 
the appendages) are then determined from: 


Cer) — (Ram = Rye : Rypm)/C AP ay seay Scp) 


The value of Cpr) divided by the frictional resistance coef- 
ficient of the canoe body Cycgm) is then again plotted against 
Fy‘/Cecamy- Note that in calculating Cyc, the full length of the 
waterline of the canoe body should again be taken. In subtrac- 
ting out the influence of the viscous resistance of keel and 
rudder and any other appendage on the measured upright resis- 
tance there is no need to utilize a shorter length of flow. 


The results of an analysis of this type is given in Appendix 2 
(again for the model of the 12 Meter yacht "Australia II", but 
now fitted with keel and rudder), analogous to that of Appen- 
dix 1. Figure 9 shows the Prohaska plots for the calculations 
carried out in Appendix 1 and 2. 


As before, the wave resistance of the model is then found 
from: 


Rw a (Rim) = Rye Ryko ba Ryra) 


where = Rycaaay = 20cm) Vaumy-Crcnmy(] + kce)Scaem 
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From Cw = Rwan/(2eanV em Scan) 
the full-scale resistance of the canoe body and appendages is 
then determined using: 


Ry = “’pVs? (Cecea(l+kep)Scp + Ce (1 +kS, 
+ Cep(1 +kp)Sp + Cy-Scp) 


3.5.3 Experience Obtained with the Prohaska Analysis 
Method 


It should be realized that the basic resistance formulation as 
used by Prohaska (1966), based on work by Hughes (1966), 
Viz: 


Cy = Ce(1+k) + Cy 


requires that no Separation is present on the model since the 
viscous resistance, as defined here, is simply a constant incre- 
ment of C, over the speed range considered. This increment 
only accounts for the change in C,; due to the elongated and 
curved streamlines over that of the equivalent flat plate and the 
viscous pressure resistance. No account can be made in this 
way for the effects of flow separation. Although most sailing 
yachts have comparatively very Jow block coefficients (usually 
in the range between 0.35 and 0.45), some canoe bodies will 
display important flow separation phenomena, often associated 
with the significant fullness of the canoe body in the aftbody, 
when a bustle or large skeg is present. 


Experience reveals that it is a requirement of the Prohaska 
principle that for sailing yacht forms the C,/C, plot on Fy‘/C, 
essentially constitutes a straight or slightly concave curve for 
F,‘/C, values between-0 and 4. If flow separation is present 
the plot will be of a convex nature such as shown in Figure 10. 
In that case, the analysis needs to be carried out utilizing a 
separation coefficient, e.g. such that: 


Cy = b + cy.Cp 
viz: Cy = Cp(1+k+c,) + a.Fy* + b 


Unless the model is very large, i.e. the model scale factor 
very small, there is little point in extrapolating the form factor 
and the components of resistance resulting from an analysis of 
this type. This is due to the fact that Reynolds number has a 
strong effect on separation and due to the fact that such flow 
behaviour is erratic, changing abruptly with small modifica- 
tions in trim for example. In sailing yacht design and research 
it is better practice in such a case to attempt to eliminate the 
flow separation, by decreasing the volume of the aftbody for 
example. 


Interesting results of various analyses of resistance tests using 
the Prohaska principle have been published by DeBord and 
Teeters (1990) and by Teeters (1993). The first of these publi- 
cations refers to experience gained in the USA with this 
method when applied to 12 Meter class yachts. The second 
paper gives results of various studies for the new America’s 
Cup Class. The method of subtracting the viscous resistance of 
appendages from the total measured resistance to obtain quasi- 
measured values for the canoe body at low speed is here 
referred to as “stripping”. Additional formulas are given to 
‘carry out this process. 


On the basis of many analyses of this type it would appear 
that when no separation is present on the model the value of 


the constant a in Cy = a.F,‘ varies between 0.07 to 0.15, 
valid up to about Fy=0.3, depending on the form of the canoe 
body. The value of the form factor itself would seem to vary 
between 1.07 for very slender, low-block canoe bodies, to 
about 1.15 for relatively high-block forms. Most well-designed 
yachts have values for 1+k_, between 1.07 and 1.11. 


3.5.4 Analysis of Side Force Measurements 


The analysis of side force measurements is necessary to 
determine whether or not non-linear lift terms are present {see 
Paragraphs 2.2 and 2.3.6.2) and to ascertain the value of the 
effective aspect ratio of the yacht as a whole when considered 
as a lifting body. To this end the total lift can be plotted 
against the angle-of-attack Bcosd for constant values of speed 
and heel angle. The tangent drawn to these curves, at zero or 
small side force values, can then be associated with the linear 
lift of the yacht as a whole. It should be noted that the lift is 
not always obtained from the measured side force by dividing 
the latter by cosé. When the model has appendages not orien- 
ted vertically at zero heel, such as winglets, this is not the case 
and care should be taken in calculating the appropriate value of 
the lift on which to base the value of the effective aspect ratio 
and the induced resistance. 
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Figure 10. Example of a Prohaska plot from which no ac- 
curate value of the form factor can be derived. The "convex" 
nature of the plot indicates flow separation is present. Visual- 
ization studies carried out later for this model (a 1:3.3 scale 
model of the canoe body of the IOR maxi yacht "Kialoa V"), 
revealed flow separation to occur in way of the skeg. 
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When, as is usual at non-zero heel angles, it is found that the 
lift is not zero at zero angle-of-attack, the zero lift angle 
(Bcos@), should be noted. For this purpose the analysis should 
preferably commence by plotting side force or lift as a function 
of angle-of-attack, similarly to Figure 2. 


The linear component of the side force on the model can be 
found by inspection using the relation: 


SFornimy = eV pny (OC iry/9Q),. (Bcosd-(Beosd)o)Ay ryCOS@ 


= total linear side force at a particular 


where SF ota) 
Vay 8 and ¢; 


(dC, +y/da), = linear lift-curve slope of model as a 
whole at a particular Vay) and ¢; 
(BCOS) zero lift angle-of-attack; 


and Avro) total lateral area of model. 


The linear lift-curve slope factor can be calculated from: 
(8C\ry/4a).= SFarnan/(AP cay V any (BCos-(Bc0s))A tymyCOS$) 


which relationship can be used for any value of SF yyy) and 
Bcos@ if the side force curve is sufficiently linear. 


When the linear lift-curve slope for the various heel and 
model speed combinations is known, the effective aspect ratio 
of the total model configuration can be obtained using the well- 
known Whicker and Fehlner formula (1958), viz: 


(8Cyqy/90), = 27ARM(2 + (ARZ + 4)!) 
that is: AR, = 82(8C,q,/4a),/(42° - (@C,7y/00),”) 


This effective aspect ratio is later used for the analysis of the 
induced resistance (see Paragraph 3.5.6). 


The coefficient of the non-linear component of the side force, 
dependent on the square of the angle-of-attack, can be deter- 
mined by inspection for every combination of speed and heel 
angle, using: 


SFoomum) = YapeyV amy -Croztn-ArtimyCos 
where CLaty = CLyzty-(Bcosé-(Bcos@),)? 
from which follows that: 
Create = SFartnan!/('42eanV ae -Arta-(Bcos@-(Bcosd),)?.cos¢) 
Likewise, the sometimes significant loss in side force that 
occurs due to flow separation, which can usually be considered 


to be a function of the third power of the angle-of-attack, can 
be written as: 


SFyctuomy = 20ayV ae) -Crarn-Avtim-COS@ 
where Cunth = Ctanty-(Bcos¢-(Bcos¢)y)* 
from which follows that: 
Conta = SFaeraan/CA Pay aay -Arraw-(Bcosd-(Bcosd))’.cos¢) 
An alternative to this tedious inspection method for the 
analysis of the side force consists of utilizing a multiple regres- 
sion analysis technique or a least-squares method with which 


this analysis is performed quickly and accurately. In the latter 
Case it is necessary that the dependencies on heel angle, Froude 


number and angle-of-attack are known. This is generally not 
the case so that a multiple regression analysis method is always 
to be preferred. The following general model is best used for 
this purpose: 


(Curnen)man = a0 + a6? + a3Fy + aFy? + ash.Fy 
+ (by + b,@ + byd? + b,Fy + DF? + bob.Fy).Bcosd 
+ (Cy + C6 + Ch? + Fy + CyFy? + €56-Fy).(Bcos¢)? 
+ (do + did + dd? + Fy + dF? + do@-Fx).(Bcosd)? 


Here, Cypeoyraamy = a:b + a.¢? + asFy + agFy? + ac6.Fy 


(OC, 74/40), = bo + bd + bg? + bsFy + DFy? + Bsh-Fy 


Craaty = Co + CO + Ob? + C3Fy + CyFy? + C56.Fy 


and Cute = dy + do + Gop? + G,Fy + dyFy? + dob.Fy 
The data file used to generate this mathematical model must 
list Fy, @. 8 or Bcos@, and SFr myC2e (my Vey -Arri)-COS¢). 


The introduction of the terms independent of Bcos@ are re- 
quired to model the asymmetry of the model associated with 
heel angle, i.e. to model the fact that the side force is general- 
ly not equal to zero at zero leeway, for non-zero heel angles. 
The value of Cy,¢~oyrumy thus determined is the value of the lift 
coefficient of the complete model at non zero heel angles, 
when the angle-of-attack is zero. 


An example of the calculations involved in this analysis is 
given in Appendix 3. The case considered as an example is 
that of the model of the IACC yacht already referred to in 
Figures 2 and 5. 


A note of caution is here necessary with respect to the use, in 
multiple regression, of a mathematical model with a large 
number of independent terms (variables) relative to the number 
of data points, such as employed in Appendix 3 (and Appendix 
4). This can only result in meaningful results if a so-called 
forward selection method is used together with a high value for 
the F statistic. In forward selection, the first variable con- 
sidered for entry into the equation is the one with the largest 
positive or negative correlation with the dependent variable. 
The F statistic for the hypothesis that the coefficient of the 
entered variable is 0 needs to be calculated. To determine 
whether this variable (and each succeeding variable) is entered, 
the F value is compared to a prescribed criterion. In the mul- 
tiple regression routine included in the well-known SPSS suite 
of statistical programs (Statistical Package for the Social Scien- 
ces, available from SPSS Inc. in The United States), one crite- 
rion is the minimum value of the F statistic that a variable 
must achieve in order to enter, called "F-to-enter". Another 
criterion associated with the F statistic, is called "probability of 
F-to-enter". Minimum values for both these F levels need to be 
prescribed. On choosing sufficiently high values thereof only 
those independent variables will be selected that have a very 
high degree of correlation with all of the test data, yielding a 
sufficiently accurate mathematical model. 


From the results of this multiple regression analysis, cross- 
curves for the effective aspect ratio and the four lift coef- 
ficient components can be drawn as a function of angle-of- 
attack, model speed and heel angle, as shown in Figure 11 for 
the example considered in Appendix 3. Again it should be 
noted that the characteristics of the JACC yacht here con- 
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sidered as an example, which has a length-beam ratio on the 
waterline of 4.70, display less dependence on Froude number 
than do many other yachts, particularly with respect to the 
coupling between Froude number and heel angle. 


Although both the non-linear lift and, at larger values of the 
angle-of-attack, the subsequent loss in lift are strongly as- 
sociated with viscosity, and therefore strongly dependent on 
Reynolds number, the total side force can only be extrapolated 
to full scale by assuming a Froude dependence only. Accor- 
dingly, the full-scale value of the side force for a given boat 
speed Vz, heel angle @ and leeway angle 6 follows from: 


SFu, = “’pV~e". Arcos Cy scorn + (8C.7/dx),.Bcos¢ 

+ Comm (Bcosd)? + Cram (Bcosd)*), or 
SFy, = “’pV,'-Apcoso((8C, 7,/da),(Bcosp-(Bcos¢),) 

+ Carta Bcosd-(Bcosd)))? + Craemn(Bcosd-(Bcos¢)y)*) 
if the inspection method is used. 

As the angle of heel and leeway increase, the separation 

phenomena associated with the non-linear component of the 
side force will change. Usually a decrease in this side force 


component is found due to the more rounded shape of the bilge 
of the canoe body presented to the incoming flow, from which 
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Figure 11. Results of the analysis of the side force on the 1:7 
scale model of the IACC yacht carried out in Appendix 3. All 
four components, and their dependency on angle-of-attack, heel 
angle and Froude number, are indicated. The inset figure gives 
the value of the effective aspect ratio when the model, as a 
whole, is considered as a lifting body. 





the flow is seen to separate more gradually then from the 
sharper forefoot or centerline crease, when upright. This can 
be seen in Figure 2, for 30 degrees of heel, for which the non- 
linear lift component is less important than for 10 degrees of 
heel. 


3.5.5 Resistance Due to Heel 


To determine the resistance due to heel it is necessary to 
extrapolate the curve of measured resistance as a function of 
angles-of-attack, for a given model speed and heel angle, to 
zero side force. This is carried out by plotting the total meas- 
ured resistance as a function of the square of the side force. 


As discussed in Paragraph 2.3.5, the resulting curve will 
generally not always be a straight line (as is often assumed) if 
in the earlier analysis of the side force a significant non-linear 
lift component was found to exist. Accordingly, the plotted 
curves for each of the speed-heel combinations will generally 
reveal an increasing slope as the value of the square of the side 
force increases. At low values of the side force however the 
non-linear lift term and the associated drag due to (non-linear) 
lift will generally be negligible and a straight line can usually 
be defined to pass through this part of the curve which can be 
extrapolated to obtain a value for zero side force. 


On subtracting from the thus derived value of the resistance 
the known upright resistance of the model for the same speed, 
the resistance associated with heel only can be found. 


Extrapolation to the full scale of this component of resistance 
is usually carried out by assuming dependence on Froude 
number only (the velocity-squared law). This assumes that this 
component is primarily of wave-making origin. It is known 
that this is often not the case and that the change in wetted sur- 
face area with heel is the primary cause of this resistance. 


No simple solution exists to the problem of how to scale the 
resistance due to heel or indeed the resistance due to the at- 
titude of the hull. When the resistance due to heel is "negative" 
it can safely be assumed that the net effect of heel is associated 
with reduced wetted surface area, in which case this com- 
ponent of resistance can be scaled by: 


Ryy = ’pV,g".Cecp(] +kcp)Crgv-Scp 


requiring full-scale values for Crcg and Scg. Here the coef- 
ficient C,,y is interpreted to represent the ratio of the wetted 
surface area of the hull at the heel angle under consideration to 
that at zero heel. This coefficient is determined from the 
measurement data using: 


Crov = Revow/(4 aa V een -Crca(] +Kep)Scecmy) 


When it is probable that the resistance due to heel is involved 
with a change in wave-making, the resistance due to heel is 
scaled according to: 


R, = ’pVe -Cro-Scp 


where Cry = Rocay(20 em Veo Scam) 

Similarly to the analysis of the side force, an efficient method 
to carry out this resistance-due-to-heel analysis is to adopt a 
multiple regression analysis technique. This can be carried out 
best by determining the total resistance coefficient associated 
with the development of side force and to find the best fit to 
the data, adopting the following general mathematical model: 
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(Cornea -Cernenye=o)mah = &16 + e.g" + @sFy + egFy? 

+ es.Fy + €gh?.Fy + €76-Fy? + €967.Fy? + egFy' 

+ €:90-Fy* + €),67.Fy* + (fo + £16 + fhe? + Fy + fF? 
+ f6b.Fy + f56?.Fy + f56-Fy? + f967.Fy? + foFx* + fiod-Fy! 
+ f6?-Fy®). (Cusp commen + Cy ry/80),-BC08$) yy) @AR,) 

+ (Bo + £10 + 226° + BaF + Bay? + B56-Fy + B6¢7.Fy 
+ 856.Fy? + B6?.Fy’ + BoF N+ 2100-Fy'+ £1197 -Fy*).(Bcos)’ 
where (Cig-oyrnam + (@Cr/9e),.BCOSd) = a6 + ang? + a3Fy 
+ ajFy? + ash.Fy + (bo + b)d + by¢? + DFy + DAF? 

+ bo.Fx).Bcosd 

Here the resistance associated with heel is: 

Crm) = C16 + C26? + O3Fy + egFy? + €56.Fy + €6¢?.Fy 

+ e.Fy2 + egh?.Fy? + Fy + €196-Fy* + €1,¢7-Fry* 
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Figure 12. Results of analysis of the resistance at heel and 
leeway, measured on 1:7 scale model of LACC yacht (see also 
Figure 5). This analysis, carried out in Appendix 4, reveals the 
significance of each of the three resistance components and 
their dependence on angle-of-attack, heel angle and Froude 
number. 


The remaining part of this general formula for the resistance 
associated with the development of side force is explained in 
the following paragraphs. An example of this analysis is given 
in Appendix 4. 


3.5.6 Induced Resistance 


3.5.6.1 Induced Resistance Associated with Linear Lift 


The induced resistance is associated with the induced flow 
angle. It is represented by the straight line drawn through the 
data points at low values of the side force in the plot of the 
resistance versus the square of the side force, discussed .in 
Paragraph 3.5.5. 


The general formula for the induced resistance associated 
with linear lift is that given in Paragraph 2.3.6.1 which, for 
the model as a whole, becomes: 


Rum) = Yo pay gem -Atomy-Cratamy +0)/(7AR,) 


where Cuotn = (6Cyry/0),.BCosd 

With the value of the effective aspect ratio determined from 
the side force measurements, it is possible to determine the 
value of the 1+o0 factor, accounting for the deviation of the 
model as a lifting body from a body with an elliptical spanwise 
loading distribution, for each combination of Vgq4) and heel 
angle @, according to: 


1+o = Ryw-TARS (A payVamy-Artiowy-Cremmy) 


The extrapolation of this resistance component to the full 
scale is carried out according to: 


R, = 4pVep"-Apr-Crara“(1 + 0)/(@AR,) 


The part of the mathematical model given in Paragraph 3.5.5. 
associated with induced resistance due to linear lift is: 


Rymy = (fo + fo + fp? + fFy + f)Fy? + £56-Fy + f.¢? Fy 
+ f)6.Fy? + f96?-Fy? + foFy* + fiod-Fy* + f11¢?-Fy‘). 
A(Cygeornany + (OC, /8e2),-BCOSH) man) (AR) 


An example of this analysis is given in Appendix 4 for the 
1:7 scale LACC model referred to before, the results of which 
are also shown in Figure 12. The value of 1+0 as follows 
from the analysis, shown in Figure 13, is seen to be a function 
of Froude number. Usually this efficiency factor is a function 
of both heel angle and speed. 


3.5.6.2 Induced Resistance Associated with Non-Linear Lift 


It will not always be possible to distinguish between the 
induced resistance associated with non-linear lift and the resis- 
tance increase associated with the attitude of the hull, i.e. the 
increase in wave or viscous resistance of the canoe body (due 
to excessive leeway angle for example), otherwise not ac- 
counted for. This is usually not a serious problem since the 
scaling of these components of resistance must by necessity be 
according to Froude’s.law of similarity. Accordingly, there is 
sometimes little to be gained by differentiating between these 
resistance components. 


When the analysis of the induced resistance due to non-linear 
lift is of interest however, this can be undertaken by assuming 
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that this component is proportional to the cube of the angle-of- 
attack, viz: 


‘ast 2 
Ria) oe 2eayVam -Craamy-Artim) 


where Cram) = Crioa(Bcose)’ 
The value of cp,,. is determined for every combination of speed 
and heel angle from: 


Crier = Ryaray/CZ 0m Veuy-Arran(Bcosd)*) 


As already stated, this resistance component is strongly 
dependent on Reynolds number because of the viscous nature 
of the flow separation involved. Scaling rules are not available 
and there is no choice but to extrapolate this component accor- 
ding to Froude’s Law of similarity. Again, however, as long as 
the model is not too small the errors made hereby should not 
be significant unless the canoe body and not the appendages 
produces most of the side force. The full-scale value is deter- 
mined from: 


Rian = ’pVp"-Art-Crio2(Bcoso)? 


in which the full-scale value for A,,; is used. The resistance 
coefficient Cp),. must be considered as independent of scale. 


The part of the general formula given in Paragraph 3.5.5 for 
the analysis of the total resistance of the model associated with 
this component, and with R,,, is given by: 

Crtoam) + Cram) = (Bo + 816 + 26? + B3Fy + 2.F,? 
+ g50.Fy + 26?.Fy + 876-Fy? + g9¢?.Fx? + BoFy* 
+ £100-Fy* + 21:67.Fx‘).(Bcose)? 
In appendix 4 the induced resistance associated with non- 


linear lift-is found not to be significant enough for the used 
regression analysis method to select any of these terms. Only 
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Figure 13. Value of the induced resistance efficiency factor 
1+o for the 1:7 IACC model considered in Figure 12. Values 
greater than 1.0 indicate a greater induced resistance than that 
associated with elliptical spanwise loading (see Appendix 4) 


the component of resistance associated with the attitude of the 
hull is shown to exist. 


3.5.7 Analysis of Resistance Due to the Attitude of the Hull 


As explained in Paragraph 3.5.6.2 the resistance associated 
with the generation of side force that remains after subtraction 
of the resistance due to heel at zero leeway and the induced 
resistance due to the linear component of lift, can generally not 
be isolated when the R,,, component also exists. This is be- 
cause its analysis must normally be undertaken by adopting an 
expression also possessing terms involving the cube of the 
angle-of-attack, similarly to that for C,),9. 


By first principles, however, the analysis of this resistance 
component is not identical to that described in Paragraph 
3.5.6.2 for R,.o. In general no specific relation with side force 
or angle-of-attack may be presumed, and a resistance coef- 
ficient for every data point needs to be defined using: 


a 2 
Cran) ans Rermy-Rvemy-RweyRoomyRiany Rica)! (Ap. om) ¥ Bow) Aur) 


The resulting values for Cram) are then cross-faired on the 
basis of “4payVau-Arrmy-Bcosd, setting out different curves 
for Froude number and heel angle in particular. The result 
thereof can then be investigated with respect to determining 
whether a quadratic, cubic or other relation in Bcos¢ exists, 
and the dependency thereof on heel angle and Froude number. 
In the case the resulting relation is a quadratic in the angle-of- 
attack it is possible that a mistake will have been made in 
determining the induced resistance associated with linear lift, 
as this is a quadratic function with respect to angle-of-attack. 
In the case the resulting curve is a cubic in the angle-of-attack 
the associated resistance can be linked with the induced resis- 
tance due to non-linear lift if this was found to exist using the 
analysis procedure of Paragraph 3.5.4. When this is not the 
case and a third or higher degree function in the angle-of- 
attack presents itself, it is likely that this resistance component 
is associated with changes in wave resistance or in viscous 
resistance (particularly as a result of flow separation) or in 
both, caused by the leeway angle or other hull attitude varia- 
bles generating the side force. 


Extrapolation to the full scale is carried out by considering 
the resistance coefficient independent of scale, viz: 


Ra = YpVs-Cran-Art 


using, for every speed, heel and leeway angle of interest, the 
cross-faired curves for Cpaymy a8 a basis. 


In appendix 4 an important R,, component is found to exist 
for the 1:7 IACC model there analyzed. Here the angle-of- 
attack Bcos@ at which this component presents itself is higher 
then about 7.5 degrees, as can be seen from Figures 5 and 12. 
Insufficient data points are here available to detect any re- 
lationship with heel or leeway angle. Even its relationship with 
the angle-of-attack could not properly be defined using regres- 
sion analysis. A cubic function in angle-of-attack was found by 
inspection. 


3.6 Determination of Transverse Stability from Model 
Tests 


Irrespective of the test technique adopted, the transverse 
stability characteristics of the model can be ascertained from 
the following moment equation, for the equilibrium condition: 
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Awy:-GZiay = Mueexmy + rr - Zerumy)SFi/cose 


where Ay 


(M) 
heel(M) 


total weight of displacement of model; 

righting moment lever of model; 

external hee! moment exerted on model to 

obtain the required angle of heel; 

height of tow force exerted on model above 

the bottom of the canoe body or some other 

reference point, measured in the plane of 
symmetry; 

Zcrum) = height of center of resultant side force 
acting on canoe body and appendages 
above the bottom of the canoe body or 
some other reference point, measured in 
the plane of symmetry; 

= hydrodynamic side force acting on model. 


tou 


" 


ZtF 


and SF) 

In the case the external heel moment exerted on the model is 
caused by moving a weight transversely across the deck, this 
moment is equal to: 


Mreevm) = mg-h,.cos 


where m = mass of weight displaced from the center- 
line of the model to bring about the re- 
quired heel angle; 

and h, = transverse distance the mass is moved from 


the centerline of the model. 


In the free-running test set-up no external heeling moment 
acts on the model when the height of the tow force is at the 
estimated center of effort of the resultant sail force and the 
location of the center of gravity of the model is properly 
modelled. In that case Mreeymy iS zero. In the case of a very 
low position of the center of gravity of the full-scale yacht it is 
not always possible to model the location of the center of 
gravity, in which case an external heeling moment needs to be 
applied in this case as well. Here, however, this heeling 
moment needs to provide the model with a greater stability. 


As is explained in Paragraph 4.8.1.2 the transverse righting 
moment at speed is usually different to the (hydrostatic) 
righting moment at rest. Accordingly, it is necessary to ex- 
perimentally determine the height of the center of the resultant 
side force on canoe body and appendages and to consider the 
righting moment of the model as the unknown for which to 
solve the heeling moment equation. The measurement of the 
location of the side force is normally carried out by utilizing 
multi-component Strain gages, allowing the forces and 
moments in 2 or 3 orthogonal planes to be determined. If the 
vertical location of the center of side force cannot be experi- 
mentally determined, both the righting moment and this side 
force location can be approximately determined from test data 
analysis by using a least-squares method. 


Rather than solve the heeling moment equation given above 
for the total righting moment lever GZ), it is better to solve 
directly for that part of the righting lever that is influenced by 
speed. For this purpose this equation needs to be re-written as 
follows: 


4an-MNiwy-Sind + SFoyy-Zcrimy-COSS = Moccia 
~ Awy-(KBay) + BM - KGoyy)-sing + SFy)-Z¢-cosd 
Here, the lever of residual stability MN,).Sing (and sometimes 


also Zcrimy) is unknown. The value of all other parameters is 
either determined from the measurement data (Aw), SF), 


Mreeumy> KGa, Zre and ) or is obtained from hydrostatic 
calculations (KB, and BM). 


Comparison of the experimental value of MN.) with the 
corresponding hydrostatic value, obtained from calculations, 
will yield the required information with respect to the loss of 
transverse stability with speed for the various heel angles 
considered. 


The experimentally-determined values for MN) and Zcrum) 
need to be available as functions of Froude number and heel 
angle to facilitate interpolation when calculating the full-scale 
performance of the yacht. This is again best carried out by 
using the multiple regression procedure described in Paragraph 
3.5.4, viz: 


MNowyman = hie + hod? + hyFy + hyFy? + hoo.Fy 

+ (ip + i + ip? + Fy + igFy? + ip6.Fy).Bcosd 

+ (io + JO + joe? + jsFy + jaFy? + i5¢-Fy)-(Bcosd) 

+ (ky + ky + k,@? + kyFy + kyFy? + ke.Fy).(Bcoso)? 
Zcrum)mah = 1) + 16? + Fy + Fy? + 15¢-Fy 

+ (mp) + m@ + m¢? + mFy + mF,y? +m.¢.Fx).Bcosd 
+ (ny + nd + no@? + nFy + n,Fy? + nsb.Fy).(Bcosd)* 
+ (09 + 0,6 + 067 + 0,Fy + O,Fy? + 056.Fy).(Bcoso)? 


Extrapolation of these two entities is simply carried out by 
multiplying the model values by the scale factor ). 


3.7 Final Performance Prediction 


The final performance prediction for a complete range of true’ 
wind speeds and true wind angles is carried out using a similar 
procedure to that described in Paragraphs 4.9 and 4.10 for the 
VPP assessment of the resultant sail force and the subsequent 
calculation of the equilibrium true wind speed and true wind 
angle. Here, however, the experimentally-determined side 
force, resistance and stability values are used instead of the 
VPP values described in Paragraphs 4.3 through 4.8. The 
procedure to be followed will now be described. 


Starting with the lowest of about 10 or more boat speeds 
(depending on whether or not a computer program is used for 
this procedure), covering the range of interest, the model side 
force value is first calculated from the expressions (or graphs) 
obtained by means of the procedure described in Paragraph 
3.5.4. This has to be executed for a large number of leeway 
angles (20 or more) to cover all the points of sail, starting with 
zero leeway angle and zero rotation angle of the rotatable 
appendages, corresponding with the straight-downwind point. of 
sail. An initial, zero heel angle value is to be adopted for each 
leeway angle case. The found side force value is then scaled up 
to full scale. The corresponding model values for MN and Zcr 
are then calculated from the expressions (or graphs) discussed 
in Paragraph 3.6 and also scaled up. A new value of the 
equilibrium heel angle can then be found from: 


tan@ = SF(Zcp - Zep) KA(KB + BM + MN - KG)) 
in which SF, MN and Z.z, follow from the model test results 


as described above, the height of the center of the resultant:sail 
force Zc_ from VPP calculations (or otherwise), and KB, BM 
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and KG from the hydrostatic and weight calculations for the 
full-scale yacht. 


Since the heel angle found will, in general, not be equal to 
the initial value assumed, new values for SF, MN and Zcp, 
need to be determined for this new heel angle value, thereby 
initiating an iteration procedure. After 3 or 4 iterations conver- 
gence is usually reached yielding final values for side force and 
heel angle, for the considered Froude number and leeway 
angle. With these values the total model resistance is calculated 
from the expressions (or graphs) discussed in paragraphs 3.5.5 
through 3.5.7 and also scaled up for the full-sized yacht. 


Once the full-scale side force and total resistance for a par- 
ticular set of values for forward speed, heel and leeway angles 
are known, the resultant sail force is then determined utilizing 
the IMS sail coefficients or other sail force formulations. The 
procedure to be followed to obtain the resulting true and 
apparent wind speeds and angles is described in detail in 
Paragraphs 4.9 and 4.10. 


When the free-running test set-up is utilized the procedure 
described above can also be adhered to. Since this test method 
yields the same (unique) leeway - heel angle relationship as 
that for the full-scale yacht (discussed in Paragraph 3.2), the 
calculation of full-scale performance is simplified to some 
extent in that the specification of a set of values for speed and 
leeway angle immediately leads to the corresponding heel angle 
without having to solve the heeling moment equation. For this 
test technique alternative procedures have been used by Van 
Oossanen (1985) and by Henriksen and Jensen (1989). 


3.8 Noteworthy Model Test Results in the Public Domain 


Results of tests with models of sailing yachts are not often 
published or become generally available. The reason for this is 
that most experimental studies are performed for racing yachts 
and are carried out on behalf of yacht designers and naval 
architects specializing in this field. On publishing the results 
obtained, the insight gained will become common knowledge 
and any advantage that a particular designer will have obtained 
through model testing over others, will be lost. 


While this is an understandable position for designers to 
adopt, the recent trend for some workers in the field to make 
sweeping Statements and important conclusions in papers on 
results of sailing yacht research without providing sufficient 
results for others to corroborate or otherwise, is a most deplor- 
able practice. The custom of presenting graphs and figures 
without providing scales of reference, the total lack of: suppor- 
ting results altogether, or the practice of not providing infor- 
mation about the underlying methods and assumptions made, is 
unacceptable when addressing professionals or the scientific 
community. 


Nevertheless, some very useful papers have been published 
over the years on towing tank work for sailing yachts. First 
and foremost are those published since 1980 by the Delft 
University of Technology on the results of the Delft yacht 
model series. This series now forms a very important data 
base, in particular with respect to resistance and side force. 
Important parts of the IMS VPP are in fact based on these 
results. The main papers on this work are those of Gerritsma et 
al (1981, 1988, 1991 and 1992/1993). 


Another noteworthy series of publications giving useful 
results of tank tests are those by Claughton and co-workers of 
the University of Southampton. Examples are the papers by 


Claughton et al (1987) and Campbell and Claughton (1987). 


The work carried out at different laboratories on various 
models of the 5.5 Meter yacht "Antiope" also constitutes an 
important set of results. The main publications on this study 
are those by Herreshoff and Newman (1967), Kirkman (1974), 
Kirkman and Pedrick (1974), and Letcher (1975). 


Recent results of experimental research on modern 12 Meter 
Class yachts have been published by Van Oossanen (1985) and 
Van Oossanen and Joubert (1986). 


Results of important experimental studies covering ‘steering 
properties and the behaviour in waves have been published by 
Spens et al (1967), Gerritsma (1971), Pedrick (1974), Nomoto 
(1979), and Masuyama et al (1993), amongst others. 


Although important experimental studies on the new class of 
IACC yachts have been carried out, much of the papers on this 
work suffers from the problem indicated above of not actually 
providing any useful data to workers in the field. Noteworthy 
exceptions are the recent papers by this author (1992) and 
Teeters (1993). 


A careful search of the technical literature will reveal that in 
addition to the above many other useful results of experimental 
studies are available. Most of these can be found in the fol- 
lowing series of technical symposia: 


- The bi-annual "Chesapeake Sailing Yacht Symposium” (first 
held in 1974), the main organizer of which is The Society of 
Naval Architects and Marine Engineers; 

- The "Ancient Interface" Symposia (first held in 1969), the 
main organizer of which is the American Institute of Aero- 
nautics and Astronautics; 

- The "New England Sailing Yacht Symposia" (first held in 
1976), the main organizer of which is The Society of Naval 
Architects and Marine Engineers; 

- The bi-annual “Symposium on Yacht Architecture” (first held 
in 1971), the main organizer of which is the Netherlands 
National Association of Watersport Industries (HISWA); 

- The bi-annual "Fachseminar Jachtentwurf und Jachtbau" (first 
held in 1980), the main organizer of which is the Deutscher 
Boots- und Schiffsbauer-Verband. 


4. SPEED PREDICTION BASED ON 
EMPIRICAL VELOCITY PREDICTION 
MODELS 


(SEMI-) 
(VPP) 


4.1 Introduction 


For some years now semi-empirical speed prediction models 
have been extremely popular with all those involved in the 
design and racing of sailing yachts. Many design and race 
handicap problems are addressed through these computer 
programs. Their usefulness will continue for some time yet but 
it has already become clear that during the years to come this 
type of performance model will make way for the com- 
putational fluid dynamics programs that are now making their 
presence felt. 


One of two basically different methods can be adopted in 
setting up a mathematical model for the calculation of the 
speed of a sailing craft in a VPP. The first and most frequently 
adopted method is to determine the thrust and side force 
developed by the wind on the sails and then to calculate at 
what boat speed, heel and leeway angles these aerodynamic 
forces make equilibrium with the hydrodynamic forces. Here, 
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3 equations are available to find these 3 unknowns, viz: force 
equilibrium in the horizontal plane in the direction of the boat 
speed vector and in the direction at 90 degrees thereto (the 
hydrodynamic side force direction), while the third equation 
involves the heel angle and moment equilibrium of all forces in 
the vertical, transverse plane. The principle of this model is 
outlined in Figure 14. 


The second method, less often adhered to, is to assume a 
specific value for the boat speed and the leeway angle, and to 
calculate the hydrodynamic side force and resistance in a 
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Figure 14. Typical calculation scheme employed in velocity 
prediction programs starting with the calculation of the aero- 
dynamic forces. 


straight-forward manner. Knowing the boat speed and the 
resultant hydrodynamic force, the calculation of the resultant 
aerodynamic force on the sails, for equilibrium, is a relatively 
simple matter. This, in turn, permits the calculation of the 
apparent and true wind speeds and angles. On performing these 
calculations for a series of leeway angles from 0 degrees to, 
say, 12 degrees, all the points of sail ranging from running 
square before the wind to beating into the wind, are covered. 
An interpolation method then yields boat speed, heel angle, 
leeway angle, etc, at any required value of the true or apparent 
wind speed and angle. The basic principle involved in this 
model is outlined in Figure 15. 


This second method has been adopted by the author. Apart 
from the mathematically more simple process for finding the 
conditions for aerodynamic-hydrodynamic equilibrium, it 
places the hydrodynamic properties of the canoe body and its 
appendages in a more central role then does the first method. . 
This results in a favourable situation with respect to mathe- 
matical modelling since much more is known about hull hydro- 
dynamics then is known about rig and sail aerodynamics, by 
virtue of the fact that towing tank testing has yielded a wealth 
of information while, by comparison, little is still known about 
the forces on the individual sails. 


A description will now be given of the more important 
elements in the method that has been developed by the author 
since 1976 (see Van Oossanen, 1979 and 1981). 


4.2 Identification of Main Elements in VPP Modelling 


In the model developed by the author the overall step-by-step 
process involved is given below. Many of these elements are 
common to other VPP models. 


4.2.1 Input Part 


1. Input of values for all 150 parameters describing the 
geometry of canoe body, appendages, spars and sails; 

2. Specification of values for calculation of density and 
viscosity of air and water; 

3. Specification of values for selection of matrix of boat 
speed and leeway angles for which calculations have to be 
carried out; 

4. Input values for the rotation angle of rotatable appendages 
(such as trim tab, keel fin and rudder) and specification of 
whether or not the rudder angle has to be calculated to 
obtain moment equilibrium of all forces in the horizontal 
plane; 


4.2.2 Hydrodynamic Part 


Calculate first estimate of heel angle; 

Calculate side force of canoe body and all appendages; 

Calculate induced drag of canoe body and all appendages; 

Calculate viscous resistance of canoe body and all ap- 

pendages; 

9. Calculate wave resistance of hull and of any appendages 
piercing the water surface; 

10. Calculate resistance due to heel; 

11. Calculate righting moment; 

12. Calculate centers of all hydrodynamic forces on canoe 
body and appendages and calculate resultant force and its 
point of application (in x, y and z coordinates); 

13. Calculate or input the height of the vertical center of the 
resultant force on the sails; 

14. Calculate heeling moment; 

15. Calculate the effect of crew weight and position on 
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Figure 15. Calculation scheme employed in velocity prediction 
program starting with the calculation of the hydrodynamic 
forces on the canoe body and appendages. 
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18. 


righting moment; 

Calculate new heel angle by equating heeling moment to 
righting moment; 

If the determined heel angle is greater then a previously 
specified maximum value then the heel angle is limited to 
this value and sails are flattened (reducing the lift coef- 
ficient), followed by either a reefing procedure or a pro- 
cedure to reduce the foot length of the sails. The height of 
the vertical center of the sail force is reduced when con- 
ventional reefing is specified; 

Repeat steps 6 through 17 until heel angle value has 
converged; 


4.2.3, Aerodynamic Part 


19. 


20. 
2s 


22. 


23. 


24. 


25. 


26. 


27. 


Calculate a first value for the apparent wind angle and 
make a first estimate of what sails are being used (genoa 
or spinnaker, mizzen spinnaker or mizzen staysail); 
Calculate the lift and drag coefficients of each sail, based 
on the assumed apparent wind angle; 

Calculate the resultant aerodynamic force coefficient and 
associated drag angle equal to arctan(Cpy,/Cy74); 

Calculate a new value for the apparent wind angle by 
equating the resultant aerodynamic force direction op- 
posite to the resultant hydrodynamic force direction (the 
sum of the hydrodynamic and aerodynamic drag angles is 
equal to the apparent wind angle plus the leeway angle, 
see Paragraph 2.6); 

Determine which sails are being used, based on criteria 
involving the calculated apparent wind angle and the 
forward thrust provided by each sail combination, when 
more than one type of sail can be used (spinnaker versus 
jib or genoa, mizzen spinnaker versus mizzen staysail); 
Repeat steps 20 through 23 until the value of the apparent 
wind angle has converged; 

Calculate the apparent wind speed by equating the resul- 
tant aerodynamic force. equal to the resultant hydro- 
dynamic force; 

Calculate true wind speed and angle from the apparent 
wind speed and angle, and the boat speed; 

Correct the calculated true wind speed and apparent wind 
speed and angle for a standard wind gradient, to derive 
the values valid for the top of the mast where the wind 
sensing gear is fitted; 


4.2.4 Added Resistance due to Waves Part 


28. 


29. 


If the true wind is forward of abeam, calculate the added 
resistance due to waves on providing either inputs for 
wave height and period, or fetch and duration of wind 
over fetch; 

Repeat steps 12 through 28 until wind direction and speed 
no longer change and the calculated added resistance due 
to waves has converged; 


4.2.5 Interpolation Part 


30. 


31. 


Repeat steps 5 through 29 for all boat speeds and leeway 
angles prescribed (usually about 30 boat speeds and about 
100 leeway angles per boat speed), filling a number of 
two-dimensional arrays for such important entities as boat 
speed, leeway angle, heel angle, total side force, total 
resistance, the so-called "flat", "foot" and "reef" values 
for the sails, apparent wind speed and angle, true wind 
speed and angle, speed-made-good to windward, etc; 
Interpolate to find the required performance for previously 
input matrix of true wind speed and angle values; 
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4.2.6 Output Part 
32. Output of performance in tables and graphs. 


The VPP here described, as do other programs. has more 
then one mode of operation. The first mode concerns a detailed 
output mode which provides the results of all of the hydro- 
dynamic and aerodynamic calculations for each speed and 
leeway angle combinations considered. This mode is utilized 
when detailed design problems need to be addressed and not 
the overall performance of the yacht. Because of the vast 
amount of output generated in this mode, only a few boat 
speed-leeway angle combinations can be dealt with in one run 
of the program. 


The second mode of operation gives one line of output for 
each boat speed-leeway angle combination, printing the cal- 
culated values of the main parameters, such as heel angle, total 
resistance, hydrodynamic side force, apparent and true wind 
speeds and angles, the speed-made-good to windward, etc. 


The third mode of operation only outputs the results of the 
interpolated values for boat speed, speed-made-good to 
windward and true and apparent wind angles relative to the 
track of the yacht through the water, for the prescribed true 
wind speeds and angles, on which results the polar perfor- 
mance diagram is based. This program mode writes the polar 
performance data to an output file, which is read by the pro- 
gram that draws the polar diagram. 


To facilitate the editing of input files, special file editing 
programs need to be available. 


4.3 Calculation of Hydrodynamic Side Force 


4.3.1 The General Expression for the Side Force of Planar 
Lifting Surfaces 


The hydrodynamic side force of every lift-producing element 
of the yacht needs to be considered. As already discussed in 
Paragraph 2.2, the side force is the horizontal component of 
the lift force developed as a consequence of leeway or some 
other means such as camber, keel rotation, trim tab angle, 
rudder angle, etc. The basic expression used for all lifting 
surfaces, such as keels, winglets and rudders, for the side force 
in the horizontal plane, is: 


SF= pV eo (Craiip(IC_/4a)Bcosd Leg tart Craa(BCos)’)A, .cosd 


where p = density of water; 
Va = boat velocity; 
CLatip) = coefficient depending on the shape of the 
up; 


0C,/da = slope of (linear) lift coefficient curve with 
respect to angle-of-attack; 

8 = leeway angle in radians (note that Bcos¢ is 
the effective angle-of-attack); 


ty = heel angle in radians; 
Lop ta = factor accounting for the lift interaction 
: between canoe body and the appendage; 
Cia? = constant in definition of coefficient of non- 
linear lift; 
and AL = lateral area of lifting surface, equal to the 


maximum span times the average chord 
length. 


The term dependent on the square of the angle-of-attack ac- 
counts for the non-linear component of lift, described in Para- 


graph 2.2. This component is only significant when the effec- 
tive aspect ratio assumes small values. As shown by Whicker 
and Fehliner (1958), the constant in the definition of the coef- 
ficient of non-linear lift c,,, (defined in Paragraph 2.3.6.2) is 
dependent on the effective aspect ratio, the taper ratio and the 
shape of the tp. A suitable formula is: 


Creo = (CraaiyTR + 0.1)/AR, 


= coefficient depending on the shape of the tip; 

= taper ratio of the appendage, defined as the 
ratio of the chord length at the tip of the 
lifting surface to that at the root of the lifting 
surface; 

and AR, = effective aspect ratio of lifting surface. 


Here Cyaritipy 
TR 


4.3.2 The Effect of the Shape of the Tip 


From the data given by Whicker and Fehlner (1958) the 
value Of Cyacip) ANd Cyo2qip) CAN be found to be 1.0 and 1.60 
respectively for square tips, while for well-rounded tips, faired 
with circular arcs, the values are approximately: 


Crattipy = 1 - 0.135/AR2° (for AR, > 0.75) 
and Craaip) = 0.70 
4.3.3 The Slope of the Linear Lift Curve 


The slope of the linear lift coefficient curve with angle-of- 
attack is based on the expression developed by Whicker and 
Fehlner (1958), viz: 


AC, /aa = 27k, AR/(2k, + cosA(ARZ/(cosA)* + 4)") 


k, = factor accounting for the. influence of viscosity 
on the slope of the two-dimensional lift coef- 
ficient curve with respect to the angle-of-attack; 

and A = sweep angle of quarter-chord line of lifting 

surface. 


Here 


The representative sweep angle of the quarter-chord line A in 
the case that the quarter-chord line is not a straight line, is best 
determined by calculating the average value along the span. 
4.3.4 Effect of Viscosity on Lift-Curve Slope 


The effect of viscosity on the two-dimensional linear lift 
curve slope is approximately (Van Oossanen, 1981): 


ks = 1 + 0.82(t/c) - tanzy.(0.117/(t/c) + 3.2(t/c) + 3.9(t/c)’) 


where t/c = average thickness-chord length ratio of lifting 
surface; 

and 7; = half of average trailing-edge angle of lifting 
surface; 


4.3.5 The Value of the Effective Aspect Ratio 


The effective aspect ratio is considered to be two times the 
geometric aspect ratio. This is the value applicable when the 
surface from which the lifting surface projects is a flat plane, 
orthogonal to the plane of the lifting surface, that is: 


AR, = 2b/c 


where b = geometric span to the tip of the lifting surface 


where the quarter-chord line intersects the tip; 
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and c = average chord length of the lifting surface; 


The placement of a bulb, winglets or a centerboard at the tip 
of the lifting surface (or winglets protruding from a bulb 
attached to the bottom of the lifting surface) augments the lift 
considerably. Computational fluid dynamics programs, dis- 
cussed in Paragraph 5, should be used to accurately determine 
the effective aspect ratio, lift curve slope and induced resis- 
tance properties for such a complex geometry configuration, 
which values can then be introduced into the VPP in a suitable 
form to account for the performance thereof. In situations 
when this is not possible or when estimates for the perfor- 
mance thereof are acceptable, the effect of winglets, a center- 
board and a bulb on the lift of the vertical, planar lifting 
surface to which these bodies are attached can be accounted for 
by increasing the effective aspect ratio (not the lateral area) by 
using corrected values for the span thereof. For winglets or a 
centerboard attached to the bottom of the keel, the effective 
span is approximately (here written for the case of winglets): 


ber = bq + Bw(Croow/Ctipk) 
where by = span of the keel to where the top surface of 


the winglets fairs into the keel surface, ex- 
cluding the fillet at that location; 


by | = span of (one) winglet measured along the top 
surface thereof, along the quarter-chord line; 
Craw = length of root chord of winglet; 
and Cipk = length of tip chord of keel. 


For a bulb attached to the bottom of the lifting surface the 
effective span of the keel is approximately: 


Derek = bx + 0.2h un 


= maximum vertical height of the bulb, as- 
sumed located in such a way that this 
maximum height is situated close to the 
location of the quarter-chord line of the keel 
at the tip. 


where hyn, 
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Figure 16. Comparison of calculated and measured lift coef- 
ficient on the yawed, surface-piercing hydrofoil studied by 
Kuhn and Scragg (1993). 


These approximate formulae are based on experience and 
represent an average of many configurations studied. They 
should be used with care as the values relative to any particular 
geometry can deviate therefrom. 


4.3.6 The Effect of the Free Surface on the Effective 
Aspect Ratio 


When the top of a lifting surface pierces the water surface or 
is only slightly immersed (in terms of the span), an appreciable 
influence of the Froude number on the lift is found. This 
situation arises for example in the case of rudders situated well 
aft on the canoe body and in the case of leeboards such as, for 
example, found on traditional Dutch sailing yachts. This effect 
of speed is caused by the reduction in the pressure differential 
between low and high pressure side at the top of the lifting 
surface due to the change in the free surface elevation and the 
associated waves formed at, and alongside, the (partially) 
submerged lifting body. 


Using the results of studies by Hooft (1988), Chey and 
Kowalski (1963) and others, it is possible to derive approxi- 
mate formulae for the effective aspect ratio as influenced by 
Froude number. For the rudder we can write: 


fir = 1 - Ayp/(CR/4) 


where f,,, = free-surface influence factor of the rudder, 
with permissible values between 0 and 1.0; 
in which h, - = depth of top of rudder below water surface at 


the speed of interest, including effects as- 
sociated with the wave formation of the hull, 
the angle of heel, etc. 


When f;,, is greater than 0, i.e. when the top of the rudder is 
situated above a depth equal to cp/4 in meters, the free surface 
can be considered to influence the value of AR,. In that case 
the effect of forward speed is approximately determined from: 


AR. = fisr-Fare-AR 


where Fyp, = 1 + 0.422/Fy,? = (for Fyp = > 0.6413) 
= 4.246Fy,p - 0.1230 (for 0.5 < Fyp < 0.6413) 
=2 (for Fup <= 0.5) 
Fyr = Vp/(g-cp)” 
and AR = geometric aspect ratio, equal to the height of 


the submerged part of the rudder at the loca- 
tion of the quarter-chord line, divided by the 
average chord length of the submerged part. 


An useful study on the effect of the free surface on yawed 
surface-piercing lifting surfaces was recently carried out by 
Kuhn and Scragg (1993). A comparison of the results of the 
above formulae for the lift, the effect of tip geometry, the lift- 
curve slope, the influence of viscosity on the lift-curve slope, 
and the value of the effective aspect ratio, with the measured 
data published therein, is given in Figure 16, for which the 
necessary calculations are given in Appendix 5. 


4.3.7 Canoe Body-Lifting Surface and Bulb-Lifting Surface 
Interaction Effects 


The canoe body factor Leg f accounts for the fact that the 
canoe body does not constitute an orthogonal plane from which 
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the lifting surface protrudes (see Paragraph 4.3.5), and the fact 
that when the canoe body has appreciable deadrise, a certain 
amount of additional side force is induced by the keel on the 
canoe body. Various studies have revealed that the side force 
induced on the canoe body is significant in certain cases, 
sometimes as much as the side force developed by the keel 
itself. 


On the basis of the results of many model tests carried out, 
the influence of the canoe body on keel lift, and the influence 
of the keel on the lift carry-over on the canoe body, has been 
deduced to be approximately as follows: 


An effective "draft" of the canoe body, in relation to the area 
on the canoe body influenced by the flow around the lifting 
body, is defined by: 


Tepes _ 0.53c.sin(acp) 


where ac, = effective deadrise angle of the canoe body in 
way of the quarter-chord line of the lifting 
surface, relative to a point on the surface a 
distance equal to 0.5c up from where the top 
of the quarter chord line intersects the canoe 
body, measured along the girth thereof. 


This effective draft is to be compared with the local "draft" 
of the canoe body, corrected for the depth of the wave trough 
at that location. An approximate expression for this second 
effective "draft" at, for example, the keel is: 


Teper = Tes x ~ 0.75Fx*.(Lcp-cos(27.xq/Leg)) 


= draft to where the quarter-chord line of keel 
intersects the canoe body, for Fy and @ = 0; 
F, | = Froude number based on boat speed and the 
length on the waterline of the canoe body 

Leg, with a maximum value equal to 0.4; 
and Xk = distance of quarter-chord of keel at the root, 
from forward ending of waterline, at zero 

speed. 


where Tog y 


If Tope > Teper, then Tope = Teper, else Tepe = Tere 


The effective span of the lifting surface is then defined to be 
b + Tp, and the ratio of this effective span to the geometric 
span is yc3 = (b + Tep,)/b. The side force increment factor is 
then calculated, defined as: 


Lop tar = 0.3125/Fy? + (1++ep)" - 4 


In this formula the minimum value of the Froude number is 
0.25 and the maximum value is 0.40. 


Later work carried out revealed that the length-beam and 
beam-draft ratio values of the canoe body have a significant in- 
fluence on canoe body-lifting surface interaction. Further 
refinement resulted in the following expressions, with which 
Lp fact Needs to be multiplied, to account for the L¢,/By, and 
Byi/Tcg values differing from the geometry for which the 
above relationship was first derived. 


Lice fuci(geom)! iF F,7(5.935/(Bwi/Ten) + 0.01357(Bw/Tcp) 
- 11.02KBy;/Tep)”) - 0.003517 (Lep/Bw.)? + 0.6698/(By,/T cp) 


+ (8.526 x 10°). (Leg/Bw,)?.(Bw1/Tca)* - 0.08376 


Lec factigeom2 = 1.165 - 0.1063(Lep/Byy)-Fy + 2.639Fy.¢? 

+ 0.0001 885(Bwi/Tcp)*.Fy’.d? - 6.0953Fy.¢? 

- 0.098914By,/T cp). 

In the resulting expression for Leg tacts VIZ 

Teese COLS1QSIR SS 420 + en) D) Len needa Leatasean 
Fy, = Froude number based on boat speed and the length on 


the waterline of the canoe body L,,, with a minimum 
value of 0.25 and a maximum value equal to 0.4; 


By, = maximum waterline beam of canoe body; 
Tc, = maximum draft of the canoe body; 
o@ = heel angle in radians. 


Figure 17 shows the significance of this interaction for the keel 
of a typical IACC yacht for a Froude number of 0.35, a heel 
angle of 20 degrees and a range of beam-draft and length-beam 
ratio values. 


The effect of viscosity on the lift carry-over between canoe 
body and an appendage is not ‘significant. Accordingly, with 
the advent of CFD codes that can account for the effect of the 
free surface on lift and induced drag in complex three-dimen- 
sional geometries, it is now feasible to determine the above 
interaction factors for a particular configuration more accu- 
rately. This is discussed in Paragraph 5. 


tcp-K fact. 
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Figure 17. Significance of the lift-interaction between the 
canoe body and the keel (fin) of a typical IACC yacht, for a 
range of beam-draft and length-beam ratio values, for Fy = 
0.35 and a heel angle of 20 degrees. 
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The interaction between the bulb and a lifting appendage, 
when the bulb is situated at the tip thereof, is simpler to ac- 
count for than that of the canoe body described above. This is 
because the influence of the free surface in the case of the 
bulb, when submerged deeply enough, can be neglected. Many 
of the formulae derived by Hoerner (1965), for example, can 
be applied in this instance. 


Approximate results are obtained with the following formula: 
Lu tac = (1 + 0.4by/b)? 
4.3.8 The Effect of Appendage Flap Angle on Side Force 
The side force associated with angle-of-attack of a flap, in the 
case of a trim tab (or a rudder attached to a skeg) for example, 


is accounted for as follows: 


ASF = YpV,".A,.cosd.L cp fac(OC,/Aa).(0a/95)5c0sd 


where 6 = angle of the flap, relative to the centerline 
of the yacht in radians; 

and 00/05 = flap-effectiveness factor, or slope of the 
angle-of-attack curve with respect to the 
flap angle. 


The flap effectiveness factor can be expressed as (Hoerner, 
1975): 


60/05 = (00x/06)y-Ksp-COSA, 


with (da/05), theoretical value of the two-dimensional, 


flap-effectiveness factor; 


Kgr = factor to account for viscosity effects on 
flap-effectiveness factor; 
and Ay = average sweep angle of hinge line of flap. 


The theoretical value of the two-dimensional flap effectiveness 
factor, for infinitely-thin foils, can be determined from 
Schlichting and Truckenbrodt (1979, page 63), viz: 


(0a/85),, = 2/7((EC. - E))? + arcsin(E!?)) 


where E = ratio of flap area to total area of lifting surface, 
of which the flap is a part. 


The factor which accounts for viscosity is (Van Oossanen, 
1981): 


ky, = 0.75 + 0.25E°! - 7.35(1 - E®!).tanz, 
4.3.9 The Effect Of Additional Angle-of-Attack or Camber 


The side force associated with additional lifting surface angle- 
of-attack and camber, such as incorporated in winglets or 
leeboards for example, can be determined by using a corrected 
angle-of-attack, obtained by adding the additional angle-of- 
attack and subtracting the zero-lift angle-of-attack due to 
camber, from the leeway angle, viz: 


B, = B + ayy - % 


where 6, = effective value of the leeway angle to be used 
in the formula for SF and ASF; 
yr = additional angle-of-attack of the nose-tail line 
of the average section of the lifting surface; 
and OQ = zero-lift angle due to camber. 


The above formulae can be used for all appendages such as 
keels with or without a trim tab, rudders with or without a 
skeg, centerboards and leeboards, winglets, etc. 


4.3.10 Effect of Keel Downwash on Rudder Angle-of- 
Attack 


The downwash from a lifting surface such as the keel is 
usually significant enough to reduce the angle of incidence of 
the flow to an appendage situated downstream thereof by as 
much as one-half (Scott, 1974). Here too however, model test 
results have revealed that effects of Froude number and heel 
angle are present, making it difficult to develop suitable for- 
mulations for this effect. In addition, the top part of any ap- 
pendage situated well aft on the canoe body such as a rudder is 
affected by the retarded flow due to wake effects, resulting in 
not precisely knowing on what value of the dynamic pressure 
(4pV;,’) the calculation of the lift coefficient is to be based in 
determining the angle-of-attack, on having measured the total 
force on the rudder for example. 


The work to date carried out on this topic has lead to a set of 
formulae which, although not sufficient for general application, 
have worked well in some recent projects. 


The effect of the wake of the canoe body in reducing the 
dynamic pressure at the location of a rudder, when positioned 
well aft, can be approximately accounted for by modifying a 
formula given by Schoenherr and reported in Principles of 
Naval Architecture Vol. 2 (Van Manen, Van Oossanen and 
Vorus, 1988, page 159), viz: 


w = 0.10 + 4.5(Top/hyg + bp))-Cp-Cp-Byr/(Lcn-Cwr- 
(7 - 6Cg/Cwp).(2.8 - 1.8Cp)) 


where w = Taylor wake fraction; 
Cz, = block coefficient of canoe body; 
Cp = prismatic coefficient of canoe body; 
Cyp = waterline area coefficient of canoe body; 


hyp height of top of rudder below water surface at 
the speed of interest, including effects as- 
sociated with the wave formation of the hull, 
the angle of heel, etc; 

span of submerged part of rudder. 


tI 


and bp 


The induced flow angle at the quarter-chord location of a 
lifting surface, here assumed to be the keel, is (see Paragraph 
2.3.6.1): 


Aix = Lyo(] + 0x)/(2pV?.Arx)/(@AR&) 


where Ly = total lift on keel, including that induced by 
canoe body and bulb and the lift associated 
with a trim tab, but excluding the carry-over 
lift on the canoe body and bulb (if fitted). 


For moderate to high values of the effective aspect ratio, this 
induced flow angle is about half of the value reached well 
downstream. Various relationships are available in the liter- 
ature relating the induced angle of the downwash to the down- 
stream distance and pertinent geometric variables (Stanton, 
1983), according to which the following relation can be deve- 
loped for the induced flow angle at the downstream lifting 
surface, here considered to be the rudder: 


Aix R = 15x - Xp)/(3cx))"* Lxo/(4PV?.Ark)/ (TAR ex) 
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where = induced flow angle at the rudder due to the 

downwash of the keel; 
and Xx - Xp = horizontal distance between the centers of 
lift on the quarter-chord lines of keel and 


rudder; 


OK R 


It is important to use the correct value of the lift of the 
forward lifting surface in this calculation. The correct value 
thereof for the keel, fitted with a bulb, is approximately: 


Lo = LM(2 - 1/¥cp_x)-(1 + Opun_x)) 
with Ly = ’pV_?-Arx-Lep ktact-Lputr_fact-CLaxitip) (OC x/da)Bcosd 


Here the non-linear component of lift is considered not to exist 
as the cross-over flow (i.e. the flow separation at the tip), 
causing its existence, is part of the flow associated with the 
bulb and needs to be accounted for in the calculation of the 
non-linear lift on the bulb (see Paragraph 4.3.11). 


The angle a,x p reduces the angle-of-incidence of the flow to 
the rudder. This must be taken into account when calculating 
the side force thereof. 


Further work is presently being carried out involving model 
tests and mathematical analyses of the additional effect of 
Froude number and heel angle on this downwash angle. 


4.3.11 The Side Force of Canoe Body and Bulb 


The side force developed by the canoe body itself, excluding 
induced effects associated with appendages, cannot be cal- 
culated with the above formulae. This is because the canoe 
body, when considered as a lifting body, has a very small 
aspect ratio for which the lift can no longer be considered to 
be linear with angle-of-attack. The following approximate 
formula was derived from a study of the results of model tests: 


SFop = 4eVeCAn.T cp’ + 1.8Arce | Bcosd - (Bcosd), | ). 


(Bcos@ - (Bcos¢),) 


where Toy, = maximum draught of the canoe body (without 
appendages); 
and Aicgh = lateral area of submerged part of canoe body 


at zero heel angle. 


Work is presently being carried out to determine the relation- 
ship of the zero lift angle-of-attack with heel angle, Froude 
number and canoe body geometry. Preliminary work reveals 
that the following relationship has some significance: 


(BcoS¢)) = a.G7(1 - b.Fy?)) 


The values of a and b have not yet been established accurately 
enough to give them here. 


In a study presently being carried out of the side force deve- 
loped by the canoe body itself with and without a rudder fitted, 
for which a special series of model tests were performed, it 
was found that the above formula for canoe body side force 
can only be used when a rudder or large skeg is fitted aft. It 
would appear that without any appendages fitted aft, the for- 
ward part of the canoe body produces a side force linear in 
angle-of-attack as can be expected but that the aft part of the 
canoe body produces a side force in the opposite direction, the 
significance of which depends on local flow separation. A 
strong dependence on the length-beam and beam-draft ratio of 


the canoe body is suspected. As soon as the rudder is fitted, 
however, the force on the aft part of the canoe body is drasti- 
cally altered. It follows that an accurate prediction of the canoe 
body side force can only be made in presence of the rudder, 
which significantly limits the degree of flow separation. This 
work will be reported on in due course. 


When a large bulb is fitted at the bottom (tip) of the keel fin, 
such as for LACC yachts, the side force and resistance charac- 
teristics thereof need to be taken into account explicitly. The 
side force can again be approximately deduced from: 


SFyn = “2eVeCAT-Hyy? + 1-8Apuy- | cose | ).Bcosd 
4.4 Calculation of Induced Resistance 


For each of the appendages, the following equation is adop- 
ted for the induced resistance: 


R, = (SF/cosd)?.(1 +0)/('4pV_2.A,-TAR,) 


where SF = side force of each individual appendage as 
; calculated from the equations given in Para- 
graph 4.3. Note that this has to be divided by 
cos@ to obtain the value of the total (lift) 
force produced by the appendage; 
and i+o = factor accounting for the effect of the plan- 
form. 


The effect of the planform on the induced drag is basically 
that associated with taper. A suitable formula for this effect, 
for moderate values of the taper and aspect ratio, is (Van 
Oossanen, 1981): 


l+o = 1 + (0.012 - 0.057TR + 0.095TR? - 0.04TR‘)AR, 


where TR = taper ratio of the lifting surface, defined as the 
ratio of the chord length at the tip of the lifting 
surface to that at the root of the lifting surface. 


As stated before, care needs to be exercised to determine the 
correct lift when calculating the induced resistance of the ap- 
pendages. The lift to be considered is that associated with the 
induced angle of the flow across the appendage only. The 
additional lift induced on the canoe body and the bulb (if fitted) 
by the lifting surface should not be included. Accordingly, the 
same value of the lift needs to be taken into consideration as 
that defined in Paragraph 4.3.10, where the induced flow angle 
was discussed in connection with the downwash angle. 


For the canoe body a slightly different approach is required. 
Here, it is difficult to define the appropriate value of the aspect 
and taper ratios, so that it is necessary to introduce some 
assumptions based on model test results. Reasonable results, 
for various types of yachts, can be obtained on adopting the 
following expressions for the aspect and taper ratios, respec- 
tively: 


ARcg = 2Tcp/(0.75Lep) 
and TRep = 0.3 


It follows that the induced resistance of the canoe body can 
then be expressed as: 


Rica = (SFcp/cos¢)*.(1 + ocp)/(2PVp?-Arcy- TAR cp) 


Present research being carried out indicates that this expres- 
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sion is not accurate enough for most applications. Further 
refinements are presently being pursued. 


4.5 Calculation of Viscous Resistance 


As discussed in Paragraph 2.3.2, the viscous resistance con- 
sists of skin friction resistance and viscous pressure resistance 


(also called form resistance). In the case of lifting surfaces (the - 


appendages), the viscous pressure resistance is well documen- 
ted and it is not difficult to calculate the total viscous resistance 
thereof accurately enough for most purposes. In the case of the 
canoe body, however, the viscous pressure resistance is dif- 
ficult to determine accurately, although use of an average value 
of the form factor 1+k eg, = 1.09 will not lead to large errors 
in most cases. 


Various approximations are available for the form factor, as 
reported by Van Manen and Van Oossanen (1988). In most 
cases these form factor values have been derived from analyses 
of tests with models of merchant-type vessels and are not ac- 
curate for slender, yacht-like hull forms. Because of this situa- 
tion, it is in some cases more appropriate to not attempt to 
model the total viscous resistance of the canoe body but, 
‘instead, to calculate only the skin friction resistance of the 
equivalent flat plate and to consider the excess viscous resis- 
tance as part of the so-called residual resistance. The residual 
resistance is defined as the total (measured) resistance less the 
frictional resistance. This concept has been used by the Delft 
University of Technology in developing their well-known resis- 
tance polynomials. 


The viscous resistance model for the keel and similar ap- 
pendages is: 


Rvx = pV sp. Cex(1 +kx)Sx 


where. Cx. = coefficient of friction of the keel; 
1+k, = form factor of the keel; 
and Sx = wetted surface area of the keel. 


The coefficient of friction can be calculated from: 
Cre = 0.075/(log(Ryx) - 2)? - 1800/Ryx 


where R,x = Reynolds number based on the average chord 
length of the keel. 


This relationship for the coefficient of friction is valid for 
Reynolds numbers in excess of 500,000 and approximately 
describes the effect of (natural) laminar-to-turbulent transition 
of the boundary layer on the friction coefficient. In VPP 
models various options are usually included for choosing the 
appropriate state of the boundary layers (allowing for natural 
laminar-turbulent transition or the fully turbulent case). When 
the canoe body or appendages are rough, as is often the case 
for cruising yachts, a turbulent boundary layer should be con- 
sidered, without the effect of the 1800/R, term. In the case of 
very smooth appendages, when the Reynolds number is less 
than 500,000, a laminar formulation can be adopted (see also 
Appendix 2). 


The form factor 1+k for the keel and rudder is dependent on 
the type of thickness distribution used. For foils with the 
position of maximum thickness well forward (such as the 
NACA 4-digit aerofoils), the following relationship can be 
used (Hoerner, 1965): 


1+k = 1 + 2(t/c) + 60{(t/c)* 


For foil shapes with the position of maximum thickness further 
aft (such as the NACA 63, 64, 65 and 66 aerofoils), the fol- 
lowing relationship is used (Hoerner, 1965): 


1+k = 1 + 1.2(vc) + 70(t/c)* 


When the value of the form factor is known, the viscous 
resistance of the canoe body follows from: - 


Rycs = 4pV5"-Cecp-(1 +Kep)-Scp 
coefficient of friction of the canoe body; 


form factor of the canoe body; 
wetted surface area of the canoe body. 


where Crcp 
1+kcp 
and Sc 


teow dl 


For the coefficient of friction of the canoe body the same 
formula can be used as for the appendages, given above, using 
a length value equal to Leg. 


A suitable formula for the wetted area, derived by Holtrop 
(1988) is: 


Sce = Scptat-Lep-(2T cg + Bwr)-Cu'? 
with Sopp = 0.616C,, + 0.11108 + 0.245(C,,/C,,) - 0.0228 
and Cy; = 0.453 + 0.443(Cp.Cy) - 0.286Cy, 
- 0.00347(By,/T eg) + 0.370Cyp 
In these formulae, 


Cy = area coefficient of the transverse section of the canoe 
body with the greatest area; 


Alternatively, the formula derived by Gerritsma et al 
(1992/1993) can be used: 


Scp = (1.97 + 0.171(Bw1/Tep))-(0.65/Cy)!?(V cp-Lica) 
4.6 Calculation of Wave Resistance 


In practical terms, the wave resistance of the yacht at zero 
heel (i.e. at zero side force) is defined to be the total measured 
resistance in calm water, as measured during a model test in a 
towing tank for example, less the viscous resistance of canoe 
body and appendages. As stated above however, the viscous 
resistance of the canoe body requires knowledge about the 
value of the form factor. If the value of the form factor cannot 
be determined accurately, use is often made of the skin friction 
resistance of canoe body and appendages and the residual 
resistance is evaluated instead. That is, the viscous pressure 
resistance of canoe body and appendages and the excess resis- 
tance associated with the fact that the C, value of the equiva- 
lent flat plate is not equal to that of the canoe body or the ap- 
pendages, is neglected. 


Gerritsma and co-workers of the Hydrodynamics Laboratory 
of the Delft University of Technology, when developing their 
resistance polynomials, chose to only subtract the friction resis- 
tance of the canoe body (using 0.7L¢, in Rycg) and the ap- 
pendages from the measured resistance (Gerritsma et al, 1981, 
1988, 1991 and 1992/1993). Although enough low-speed resis- 
tance data was collected for all of the 39 models, they elected 
not to carry out a Prohaska analysis of the upright resistance 
measurements of these 39 models and not utilize the form 
factor values for the standard keel and rudder used (1.215 and 
1.252 respectively). Nevertheless, the polynomials developed 
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on this basis still constitute a sufficiently accurate method for 
the calculation of the residual resistance for most design stu- 
dies. 


The results of the upright resistance tests were re-analyzed 
when the results of the series II] models (models 28 to 39) 
became available, early in 1990. The final polynomials were 
recently published by Gerritsma et al (1992/1993). One set of 
polynomials are now provided for the speed range up to and 
including F,=0.45, based on all 39 models, while only the 
light displacement models were used to develop a second 
polynomial for the higher speed range from F,=0.475 to, and 
including, Fy=0.75. 


The low-speed polynomial expresses the residual resistance as 
a fraction of the weight of displacement of the canoe body, as 
a function of the prismatic coefficient C,, the longitudinal 
position of the center of buoyancy LCB, the waterline beam- 
draft ratio By./Tcpg, and the length-displacement ratio 
Lop/Vcg, as follows: 


Rg = 1000Ac,.(a) + a;-Cp + a).LCB + a,(By,/Tep) 


+ ay(Lep/V cp") + as.C,” + ag(Cp.Lep/V cg!) 
+ a,.LCB?’ + ag(Lcp/Vep'”)? + ag(Lep/Ven)’) 


The coefficients ay through a, are given in Table 2 for Froude 
number values of 0.125 to 0.45 with ‘steps of 0.025. The 
Froude number is defined-as V,/(g-L¢,)!”. 


For the higher speed range the polynomials are functions of 
the length-beam ratio L¢p/By,. the waterplane area-displace- 
ment ratio Ayp/Vcp"", and LCB, as follows: 


Rg = 1000Ac3.(Co + €)(Lcp/Bwi) + Co(Awp/Vc3”?) + c3-LCB 
+ ¢,(Lcp/By,)? + ¢5(Lep/Bw.)-(Awe/Vcp27)?) 


The coefficients cy through c, are given in Table 3 for Froude 
number values of 0.475 to 0.75 with steps of 0.025. If the 
value of the waterplane area-displacement ratio is not readily 
available, Gerritsma et al (1992/1993) provide a useful for- 
mula, viz: 


Awp!(Lep-Bw,) = 1.313Cp + 0.0371 Lep/Veq!) 
- 0.0857Cp(Lep/Ven") 


No guidance is given by the authors as to the range of values 
of the independent parameters for which these polynomials are 
valid. This can be approximately deduced however from the 
fact that the low speed polynomial is based on all 39 models 
while the high speed polynomial up to and including F,=0.60 
is based on models 23 to 39, while for the speed range from 
Fy=0.625 up to and including 0.70 use was made of models 
24, 25 and 26 and models 28 to 39. At the highest Froude 
numbers of 0.725 and 0.75 the polynomials are based on 
models 26, and 28 to 39. 


Because the resistance formulation in the IMS VPP is based 
on these polynomials, the Delft experimental data has recently 
been the subject of further analysis. As a consequence thereof, 
concern has been raised about the accuracy of these polyno- 
mials in connection with the fact that in addition to the viscous 
pressure resistance of canoe body, keel and rudder, the wave 
resistance of keel and rudder is also included in the data on 
which these polynomials were developed. For predictions for 


heavy displacement yachts this would appear to not be of great 
significance. For light displacement yachts however this effect 
is not negligible, particularly since the volume of the standard 
keel and rudder used is relatively large in relation to the 
volume of the canoe body. For that reason, upright resistance 
tests with more than half of the 39 models were recently 
carried out, without keel and rudder fitted, in an attempt to 
deduce a correction formula for the excess residual resistance 
of the standard keel and rudder when compared to more sien- 
der appendages. Although further work is required in this area, 
a preliminary formula has been prepared which has been used 
in the newest version of the IMS VPP. 


More accurate utilization of these polynomials for a par- 
ticular design can be realized by subtracting from the calcu- 


Table 2. Coefficients of the terms in the DUT polynomials for 
the residuary resistance in the speed range from Fy=0.125 to 
Fy=0.450. 














Fy a a; a, a ay 

as ae a, ag ay 

0.125 -6.7356 38.368 -0.008193 0.055234 -1.9972 
-38.861 0.95659 -0.002171 0.27290 -0.017516 

0.150 -0.38287 38.173 0.007243 0.026644 -5.2953 
-39.550 1.2196 0.000052 0.82457 -0.047842 

0.175 -1.5035 24.408 0.012200 0.067221 -2.4486 
31.914 2.2161 0.000074 0.24435 -0.015887 
0.200 11.292 -14.519 0.047182 0.085176 -2.67302 
-11.418 5.6541 0.007021 -0.094934 0.006325 
0.225 22.179 -49.168 0.086000 0.15073 -2.87868 
7.1671 8.6003 0.012981 -0.32709 0.018271 
0.250 25.909 -74.757 0.15352 0.18857 -0.88947 
24.121 10.485 0.025348 -0.85494 0.048449 

0.275 40.976 -114.29 0.20723 0.25083 -3.0727 
53.016 13.022 0.035934 -0.71546 0.039874 

0.300 45.838 -184.76 0.35703 0.33834 3.8717 
132.26 10.861 0.066809 -1.7192 0.095977 

0.325 89.204 -393.01 0.61747 0.46047 11.543 
331.12 8.5981 0.10407 -2.8152 0.15596 

0.350 212.68 -801.79 1.0873 0.53894 10.803 
667.64 12.398 0.16647 -3.0261 0.16506 

0.375 - 336.24 -1085.1 1.6442 0.53270 -1.2242 
831.14 26.183 0.23880 -2.4505 0.13915 

0.400 566.55 -1609.6 2.0161 0.26572 -29.244 
1154.1 51.462 0.28805 -0.17835 0.018446 

0.425 743.41 -1708.3 2.4358 0.013553 -81.162 
937.40 115.60 0.36507 1.8390 -~0.062023 

0.450 1200.6 -2751.7 3.2086 0.25492 -132.04 

1489.3 196.34 0.52823 1.3791 0.013577 
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Table 3. Coefficients of the terms in the DUT polynomials for 
the residuary resistance in the speed range from F,=0.475 to 
Fy=0.750. 


Fy Co Cy Cy Cy Cy Cs 
0.475 180.10 -31.503 -7.4511 2.1950 2.6896 0.006480 
0.500 244.00 -44.526 -11.155 2.1790 3.8574 0.009676 
0.525 282.99 -51.520 -12.973 2.2745 4.3437 0.011066 
0.550 313.41 -56.583 -14.420 2.3261 4.6904 0.012147 
0.575 337.00 -59.190 -16.070 2.4192 4.7668 0.014147 
0.600 356.46 -62.854 -16.851 2.4371 5.0788 0.014980 
0.625 324.74 -51.313 -15.346 2.3341 3.8554 0.013695 
0.650 301.13 -39.796 -15.023 2.0597 2.5457 0.013588 
0.675 292.06 -31.853 -15.585 1.8479 1.5699 0.014014 
0.700 284.46 -25.146 -16.154 1.7040 0.8179 0.014575 
0.725 256.64 -19.319 -13.085 2.1528 0.3483 0.011343 
0.750 304.18 -30.115 -15.854 2.8632 1.5244 0.014031 


lated residual resistance of the design of interest, at model 
scale, the viscous pressure resistance of canoe body, keel and 
rudder pertaining to the equivalent Delft models, after which 
the total viscous resistance of the subject configuration for 
which the prediction has to be made, can be added. This is 
executed as follows: 


The viscous pressure resistance of the standard keel and 
rudder of the Delft models is approximately as follows: 
Cypxim) = Yop V amy -Crxiny-kx-Sxap 


with ky = 0.215, Sxgy = 0.1539 m’, cx = 0.338 m and an 
average value for 44, equal to 1.05 x 10° m’/s, we obtain: 


Cypxa = 1-241 Vgqny’/og(322000V ay) - 2)? 


Also, Cyprm = “eo V amy -Crrmy-Ke-Seemy 
and with kp = 0.252, Spm) = 0.0550 m’, cp = 0.110 m and 
again an average value for ), equal to 1.05 x 10° m*/s, we 
obtain: 

Cypremy = 0.520 Vay’ /(log(105000V gay) - 2)? 

The viscous pressure resistance of the canoe body of the Delft 
models can be approximately determined by adopting an 
average value for the form factor of, say, keg = 0.10. The 
value thus obtained is: 


Cyvpcam) = Ye pmmyV qm) -Crcamy-Kee-Scam 
i.e. Cypcam = 3.75 V gay -Scamy/(log(1905000V pay) - 2 


which is obtained using kcg = 0.10, Leg = 2.0 m and again an 
average value for » 4) equal to 1.05 x 10° m’/s. It should be 
noted that the first 22 models have a length on the waterline of 
1.6 m while models 23 through 39 have a length on the water- 
line of 2.0 m. The wetted area of these models follows from: 


Scam = (1.97 + 0.171 By,/Teg).(0.65/Cy4)".(2V cg)” 
The complete procedure starts by calculating Cypximy, Cvprm) 


and Cypcam) for the appropriate canoe body of interest, re- 
quiring the calculation of Scgq4) and, for the Froude number of 


interest, determination of Vgm) = 4.429Fy in m/sec. For this 
same Froude number, and the displacement of the model cor- 
responding to a waterline length of 2.0 m, the residual resis- 
tance is then determined from the polynomials given above 
using the coefficients given in Table 2 or 3. From this value of 
the residual resistance of the model is then subtracted the sum 
of the viscous pressure resistance of canoe body, keel and rud- 
der, yielding the equivalent wave resistance of the canoe body 
model and the standard keel and rudder. Extrapolation of this 
resistance to the full. scale by multiplying by the cube of the 
scale ratio factor and correcting for the differences in fluid 
density and viscosity, then adding the total viscous resistance 
of the full-scale canoe body and the actual appendage con- 
figuration, will yield a more accurate determination of the 
upright resistance. 


4.7 Calculation of the Resistance Due to Heel 


The resistance due to heel is dependent on the form of the 
canoe body. Model test results reveal that some yachts have 
almost negligible resistance increase with increasing heel 
angle, everything else remaining the same, while other yachts 
show a marked increase. In part this is due to the fact that the 
wetted surface area of some canoe body forms decreases as the 
heel angle increases, up to about 20 degrees or so, thereby 
offsetting the increase in wave and viscous resistance asso- 
ciated with the less efficient streamlined form of the yacht as it 
heels. It is found that the resistance due to heel increases 
particularly when the beam-draft ratio increases. 


As discussed in Paragraph 3.5.5, the nature of the resistance 
component due to heel is usually considered as inviscid. When 
wetted surface area changes with heel however part of the 
associated change in resistance will have a viscous origin. 
Without knowing approximately how much of this resistance 
component is viscous or inviscid, no division can be made and 
different scaling principles cannot be adopted. 


Only model tests can presently yield reliable values for the 
resistance due to heel. A study of available experimental data 
available lead to the following expression using regression 
analysis: 


R, = “’pVs?.Scp-d.(0.000891Fy.(Bwi/T cp) 
+ 0.004267¢.By:/Tcp) - 0.003142) 


in which @ is the heel angle in radians. The minimum and 
maximum values for the Froude number in this formula are 
0.25 and 0.40 respectively, because of the limited data used. 


The linear term in Froude number herein is associated with 
an increase in wave resistance. The remaining terms describe 
effects associated with the wetted surface area and, as such, 
are subject to Reynolds scaling. Figure 18 gives the non- 
dimensional resistance due to heel as a function of heel angle 
and Froude number, for a length-beam ratio equal to 4.0 and a 
range of beam-draft ratios. 


The analysis of the results of the Delft Series in this respect 
revealed the following formula (Gerritsma et al, 1992/1993): 


Ry = “4pVp2-Scp-d-Fy?(0.006747(Tep/T max) 
+ 0.002517(By,/Tep) + 0.003710(By,/Tep)-(Tep/T mnax)) 


Note that these 2 formulae are not in agreement, as can be 
seen from Figure 18, where curves for Lep/By, = 4 and 
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By./Tceg = 6, for both formulae, are given. It is not clear why 
the ratio T,,/T,,,, in the Delft formula should influence the 
resistance due to heel. A change in the span of the keel alone 
will not change this component of resistance. Presumably, this 
parameter describes some feature associated with the canoe 
body of the Delft models not accounted for in the selection of 
parameters used for the regression carried out. 


4.8 Calculation of Equilibrium Heel Angle 


4.8.1 Calculation of Righting Moment 


4.8.1.1 Hydrostatic Calculations 

Since a linesplan or a table of canoe body offsets is not 
normally available when VPP calculations are carried out in 
preliminary design work, an accurate righting moment or GZ 
curve, will not be available. The stability characteristics of the 
subject vessel must therefore be calculated by approximate 
means. This will now be discussed. 


The stability moment is GZ.A, where A is the total weight of 
displacement and GZ is the stability lever at the heel angle of 
interest. This stability lever can be written as: 


GZ = (KB + BM - KG).sind + MN.sing 


where KB = vertical position of the center of 
buoyancy, for the upright condition; 


BM = metacentric radius for the upright con- 
dition; 
KG = vertical center of gravity; 
and MN.sin@ = residual stability lever. 


The residual stability lever is associated with the change in 
KB and BM from the zero heel condition to the angle of heel 
to be considered, according to the concept introduced by 
Prohaska (1947). It can be considered to be a function of the 
metacentric radius BM, as follows: 


MN = Cys.BM 


where Cps is the coefficient of residuary stability. Note that 
here the residual stability lever is defined as MN.sing and not 
as MN as is the case in the original definition of residual 
Stability by Prohaska. 


For each of the entities KB, BM and Cy., approximate ex- 
pressions were derived, as follows: 


KB = (KBog.Veg + KBagoy- Vapi + KBappo-Vapp2 + ----)/V 


in which KBcg is the vertical position of the center of buoy- 
ancy of the canoe body above the bottom of the deepest ap- 
pendage (usually the keel), Vc, is the displacement volume of 
the canoe body, KB,,,; is the vertical position of the center of 
buoyancy of appendage number i above the bottom of the 
deepest appendage, V,,,, is the volume of appendage number i, 
and V is the total volume of displacement. 


The following formulae are good approximations for the ver- 
tical position of the center of buoyancy of the canoe body, 
relative to the bottom thereof, and the metacentric radius in the 
upright position: 

KBep = (0.8415 - 0.333C,/Cwp)-Tep 


BM = BMy.Veu/¥ 


with BMgg = (0.07032 + 0.01039Cp.Cy,2).(Cp/C,)-By,2/Tep 


Expressions for the residual stability coefficient Cp, were 
determined using regression analysis for a large series of canoe 
body forms for angles of heel of 10, 20, 30 and 40 degrees. 
The above-water shape of the canoe body was found to be the 
determining aspect of the hull form in these expressions. The 
appropriate value of MN at the angle of heel of interest is 
found by interpolation. 


The results of the Delft standard model tests have also yiel- 
ded stability information. This is given in the next paragraph. 
The simple expression for the hydrostatic part of MN.sind, 
without the effect of forward speed, is: 


MN. sing = (0.0636 - 0.0196(By,/Tc,))-¢7-Leg 
4.8.1.2 Effect of Forward Speed on Stability 


On obtaining the value for GZ using the above formulae, the 
righting moment as would normally follow from hydrostatic 
considerations is determined. At speed, however, most sailing 
yachts experience important changes in sinkage and trim and, 
perhaps more importantly, an appreciable change in the shape 
of the waterline due to the wave through amidships. The water- 
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Figure 18. Example of resistance due to heel, at zero side 
force, as found from regression analysis of model test results 
for a waterline length-beam ratio of 4.0 and a range of beam- 
draft ratio values. 
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line obtained by intersecting the body plan with the wave 
profile along the hull should really be used for stability -cal- 
culations whenever the trim, sinkage or this wave profile 
differs markedly from the static situation valid for zero speed. 
When this is done and the actual waterline is found to be 
narrower than that at zero speed (which is the case when the 
hull has appreciable deadrise and flare in the region of the 
design waterline amidships) the moment of inertia, and thus the 
stability, is reduced. 


Gerritsma et al (1992/1993) have studied this loss in trans- 
verse stability using the test results of the Delft standard series, 
representative of all types of ballasted, mono-hull sailing 
yachts. For moderate beam-draft ratios the reduction in 
righting moment lever was found to be relatively small. For 
beam-draft ratios approaching 10 or more, however, this 
reduction is shown to be appreciable. For Delft Model No. 31, 
with a waterline beam-canoe body draft ratio of 15.8, a water- 
line length of 10 meters and a metacentric height GM = 2.78 
meters, the reduction in transverse stability was found to be 
29% at Fy=0.45. 


The following formula for the residual stability lever was 
derived using regression analysis: 


MNsing = (D,.¢.Fy + D3.¢7).Lep 
with D, = - 0.0406 + 0.0109(By:/T¢,) - 0.00105(By,/T cp) 
and D, = 0.0636 - 0.0196(Bwi/T cp) 


In this formula the term D,.¢? is associated with the residual 
stability of the hull at zero speed discussed above, while the 
term D,.¢.Fy is associated with forward speed. 


It is recommended that this effect of forward speed be in- 
cluded in VPP calculations for the equilibrium heel angle since 
the influence thereof can be significant. It should be noted that 
this formula is only valid up to a heel angle of 30 degrees. 


4.8.2 Calculation of Heeling Moment 


The heeling moment due to the wind on the sails at the boat 
speed, leeway angle and heel angle considered, can be deter- 
mined from the fact that for equilibrium to exist this is equal to 
the product of the hydrodynamic force acting normal to the 
plane of symmetry of the yacht and the distance between the 
point of application of this hydrodynamic force and the point 
of application of the (reaction) force on the sails. For equi- 
librium to exist, the resultant aerodynamic force must be equal 
and opposite to the resultant hydrodynamic force, constituting 
the heeling couple. The point of application of the resultant 
hydrodynamic force can be found from: 


Zor = CZsrappi-SFapp + Zgrapp2-SFapp2 + ---)/SFry 


where Zcg, = vertical center of the resultant hydrodynamic 
force above the bottom of the deepest appen- 

dage; 
Zsrappi = Vertical center of the side force on appendage 


i above the bottom of the deepest appendage; 
SF, = Side force on appendage i; 
and SFyy = total hydrodynamic side force. 


The contribution of the resistance forces to Zcg, has been 
neglected here. In the program this is properly taken into 
account by utilizing the true resultant heeling force normal to 
the plane of symmetry of the yacht. Here we assume that the 


heeling force is simply SF,,,/cos¢. 


The vertical position of the center of the aerodynamic force 
on the sails is essentially an input to the program. For each sail 
the approximate center of effort has to be input after which a 
correction is made to these values to take into account a stan- 
dard wind gradient, leading to a final value for zcg, the height 
of the center of effort of the sail plan. 


The arm of the couple constituted by the hydrodynamic and 
aerodynamic heeling forces is zcg - Zcpz. The angle of heel fol- 
lows from the solution of ¢ from the following equation: 


(Zcg - Zc_z)-SFpy/cosé = A.(KB + BM - KG + MN)sing 


Note that in this moment equilibrium equation the vertical 
center of effort of the wind on the sails is the only aero- 
dynamic parameter. 


If the righting moment of the crew has to be included, which 
is often the case when the crew weight is an important propor- 
tion of the total displacement weight (and always the case when 
the project involves a racing yacht), the following additional 
Tighting moment is to be added to the right-hand side of the | 
above equation: 


Meew = Megy(0-475B,,9, - 0.305).cosd 


where m,,. = mass of the total crew sitting on the wind- 
ward rail; 
and Byax = maximum beam of deck. 


Solution of this moment equation yields a new value for the 
heel angle in the heel angle iteration loop. If the resultant 
equilibrium heel angle found is greater then a previously 
specified critical hee] angle at which the sails are to be flat- 
tened, reefed or the foot (of the genoa or jib) shortened, the 
heel angle is then set equal to this critical value. If the sails are 
to be reefed in the conventional way the moment equation 
above is solved iteratively for the value of Zcg - Zcip- If the 
sails are only to be flattened, or if the sail area is to be re- 
duced without a change in the vertical position of the center of 
effort, this moment equation is iteratively solved for a new 
value of the heeling force. This is achieved by iteratively 
solving for a new (reduced) value of the leeway angle. 


4.9 Calculation of Forces Developed by the Sails 


Once convergence for the heel angle has been attained, the 
final value for the total side force and the hydrodynamic 
resistance at a specific boat speed and leeway angle are known. 
At this stage, a relation between the apparent wind angle and 
the aerodynamic drag angle in the plane normal to the mast 
(with the mast oriented vertically in the zero heel condition), 
can be utilized, viz: 


arctan(Ryy/(SFpy/cos¢)) + arctan(Cyz,/Ci,)= Baw. + Bcosd 


where Ryy = total hydrodynamic resistance force; 
SFr; = total side force, in the horizontal plane; 
Cora/Ciza = ratio of total aerodynamic drag force to 
total aerodynamic lift force in the plane 
normal to the mast; 
and Baw. = apparent wind angle in the plane normal 


to the mast. 


Since the values of Cpra/Cyr, and Baw, are the only un- 
knowns, an assumed value for one of these will yield the value 
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Table 4. IMS sail coefficients as a function of the apparent 
wind angle (@aw_) for spinnaker, jib or genoa, and mainsail. 
Values for a mizzen are set equal to those of the mainsail. 
Values for a mizzen spinnaker are set equal to those of the 
spinnaker and values for a mizzen staysail are considered to be 
equal to those of the jib (genoa). 

Parasitic 


Apparent wind angle Lift coefficient 


in degrees of spinnaker drag coefficient 
of spinnaker 

Baw) = 28 C,,(1) = 0.00 Copp(1) = 0.10 
Baw(2) = 41 C.(2) = 1.31 Copa(2) = 0.15 
Baw.3) = 46 C..,3) = 1.56 Copsp(3) = 0.20 
Baw (4) = 60 C4) = 1.71 Copp(4) = 0.40 
Baw(5) = 75 C..(5) = 1.69 Cops(5) = 0.70 
Bawi_(6) = 100 C,,,(6) = 1.40 Cop (6) = 1.00 
Bawi.(7) = 130 Cyy(7) = 0.83 Cops (7) = 1.10 
Baw_(8) = 150 C,,(8) = 0.50 Copy (8) = 1.10 
Baw (9) = 180 C,,,(9) = 0.00 Cres(9) = 1.10 
Apparent wind angle Lift coefficient Parasitic 

in degrees of jib/genoa drag coefficient 

of jib/genoa 

Baw) = 7 C1) = 0.000 Cyp(1) = 0.050 
Baw (2) = 15 C,(2) = 1.000 Cpp(2) = 0.023 
BawG3) = 20 C,,3) = 1.375 Cpp(3) = 0.031 
Baw.(4) = 27 C4) = 1.450 Cop(4) = 0.037 
Bawi(5) = 50 C5) = 1.430 Cop{5) = 0.250 
Bawi(6) = °60 C,(6) = 1.250 Cpp(6) = 0.350 
Bawt_(7) = 100 C,(7) = 0.400 Cop(7) = 0.730 
Bawt9) = 180 C9) = -0.100 Cpp(9) = 0.900 
Apparent wind angle Lift coefficient Parasitic 


in degrees of mainsail drag coefficient 
of mainsail 
Baw(1) = 0 Cim(1) = 0.00 Copm(1) = 0.050 
Baw(2)= 7 Cyp(2) = 1.00 — Crop(2) = 0.030 
Bouwi.3) = 9 C.(3) = 1.22 Cppm(3) = 0.027 
Byun 4) = 12 Cin(4) = 1.35 Copm(4) = 0.027 
Baw.) = 60 Cin(S) = 1.25 Copm(5) = 0.114 
Baw(6) = 90 C.n(6) = 0.96 Copm(6) = 0.306 
Baw(7) = 120 Cun(7) = 0.58 Cppm(7) = 0.671 
Bawi.(8) = 150 Cin(8) = 0.25 Copm(8) = 1.110 
Baw(9) = 180 Cim(9) = -0.10 Copm(9) = 1.200 
Apparent wind angle _ Lift coefficient Parasitic 


in degrees of mainsail drag coefficient 
of mainsail 
(fully-battened) (fully-battened) 

Baw i) = 0 Cin) = 0.00 Copm(1) = 0.027 
Baw f2) = 7 C,,(2) = 1.15 Copri(2) = 0.027 
Baw.3) = 9 Cin) = 1.40 Copm(3) = 0.027 
Baw (4) = 12 C4) = 1.55 Copm(4) = 0.027 
Baw (5) = 60 C5) = 1.44 Copm(S) = 0.103 
Bawi(6) = 90 Ci. n(6) = 0.96 Copm(6) = 0.275 
Baw (7) = 120 Cin(7) = 0.58 Copm(7) = 0.671 
Baw 48) = 150 Cun(8) = 0.25 Copm(8) = 1.110 
Bawi(9) = 180 Cin(9) = 0.10 Copm(9) = 1.200 


of the other. In the program, a first value for the Cpy,/Cir, 
ratio is adopted, with which a first value for the apparent wind 
angle becomes known. With this value, the lift and drag coef- 
ficients for each sail are determined by interpolation from the 
IMS sail coefficients. These sail coefficients are given in Table 
4. 


In accordance with the IMS model (Poor, 1986), the effect on 
the lift and drag coefficients of sail blanketing, exposed masts 
(for fractional rigs) and the easing of the jib/genoa is con- 
sidered. 


The resultant, normalized lift and drag coefficients are ob- 
tained by adding the components for the various sails. For the 


upwind case when only the jib/genoa and mainsail are used, 
we obtain: 


Cy, = (CyBKT yA + Cup-BKT y-Am-SML)/Ay 
and Cpp, = (Cpp-BKT yA; + Cppm-BKTyy-Am-SMD)/Ay 


Similar terms are added for the mizzen and mizzen staysail 
when applicable. The symbols used are as follows: 


C,; = lift coefficient of jib/genoa; 

BKT,, = upwind blanketing factor of 
jib/genoa; 

A, = area of jib/genoa; 

Ci. = Lift coefficient of mainsail; 

BKT nu = upwind blanketing factor of mainsail; 

Aa = area of mainsail; 

SML = factor for reducing lift of mainsail — 
for fractional rigs; 

Ay = reference sail area; 

Cop; = parasitic drag coefficient of 
jib/genoa; 

Com = parasitic drag coefficient of mainsail; 

and SMD = additional drag factor for exposed 


mast in the case of fractional rigs. 
The value of the reference sail area is: 
Ay = 1.3/2 + A,/1.16 


height of fore triangle; 
base of foretriangle. 


where | 
J 


This reference sail area was used to normalize the lift and 
parasitic drag coefficients of the sails in Table 4. Accordingly, 
this value is also used here. ; 


Similar expressions are used for the downwind case when the 
spinnaker is used, viz: 


Cus = (Cisp-Asp + Cum-BKT yg-Am-SML)/Ay 


and Copy = (Copsp-Asp + Copm-BKTpg-Am-SMD)/Ay 


where C.,, = lift coefficient of spinnaker if present 
(otherwise the jib/genoa is included in- 
Stead); 
Ay = area of spinnaker; 
BKT,,; = downwind blanketing factor of mainsail; 


parasitic drag coefficient of spinnaker if 
present (otherwise the jib/genoa is included 
instead). 


and Copsp 


Again, similar terms are added for the mizzen sails. 
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The induced drag of the sails is substantial. This is not in- 
cluded in the parasitic drag and must be calculated separately. 
This component of the drag is proportional to the square of the 
lift. coefficient. In addition, as the sails are usually set to 
operate at (or close to) their maximum lift coefficient, a certain 
amount of flow separation will be present which must also be 
accounted for. This component of the drag is also assumed to 
be proportional to the square of the lift coefficient of each sail. 
For the upwind case the respective, normalized induced and 
separation drag coefficients are: 


Cow = (Cy2-BKT-A; + Cun?-BKTy-An)/(@AR,s-An) 
and Cys, = (C,2.BKT,,.A;-KP, + Cy,?,BKT,,-An-KP,,)/An 


effective aspect ratio of the sail plan; 
Separation drag coefficient for jib/genoa, for 
which 0.016 is adopted; 

separation drag coefficient for mainsail, for 
which 0.016 is also adopted; 


where AR,s 
KP. 


J 


iil 


and KP, 


For the downwind case this is: 


io) 
g 
| 


= (C,2.BKTy.A, + Cin? BKT yg-A,,)/(TAR,s.An) 
and Cysy = (C.2-BKT,,.A,-KP, + C.,?-BKTyy.-Ap-KP,,)/Ay 


where KP, = separation drag coefficient of spinnaker, for 
which 0.019 is adopted. 


Equations for the calculation of the blanketing factors and the 
effective aspect ratio are included in the program. 


When the apparent wind angle is greater than 100 degrees, 
the downwind coefficients are adopted while the upwind coef- 
ficients are adopted when the apparent wind angle is less than 
70 degrees. In the range of apparent wind angles between 70 
and 100 degrees, the thrust and heeling components of the 
resultant aerodynamic force for the upwind and downwind 
cases are compared. A special performance algorithm invol- 
ving these is used for sail selection in that case. : 


The final lift and drag coefficients are: 
Cir, = C.r’f 
and Cora = Cpo + Cpp.f.r? + (Cp, + Cps).f?.1? 


where Cp) = parasitic drag coefficient of the above-water 
hull, rig and crew exposed to the wind; 


f = factor which accounts for the flattening of the 
sails with which the lift coefficient is reduced 
when the critical heel angle is reached, before 
sail area reduction commences; 

and r = fraction of sail area being used. 


When only the foot of the sails is reduced and not the hoist 
and foot together as in conventional reefing, the reduction in 
lift and drag components due to lift are reduced linearly with r 
and not with r*. It should also be noted that the parasitic drag 
coefficient Cp, is only reduced when flattening for the upwind 
case. 


The parasitic drag coefficient of the above-water hull, rig and 
crew exposed to the wind is approximately: 


Coo = (0.2 + 0.7sinBawi_)-Apce + Aprig + O-SApcrew)/An 


where Ap, = effective projected area of hull exposed to 
the wind; 
Apig = effective projected area of mast(s) and stan- 


ding rigging exposed to the wind; 


and Apcew = effective projected area of total crew. 


The effective projected area of the hull exposed to the wind is 
approximately: 


Agcy = Apepr + (Ape2 - Apeps)-SInB awh. 


where Apa = Haut+superse(Bwr + Byax)/2 
Apa = Hiyu+supersir “Cece + Lo,)/2 
Here, Hynssuperssr = representative height of above-water hull 


and superstructure. 


Approximations for Hyun+supersis> Aprig 29d Aperew are included in 
the program. 


On having obtained these values for C,,;, and Cpy,, the 
aerodynamic drag angle can be obtained, from which a new 
apparent wind angle can be determined as follows: 


€, = arctan(Cpz,/Cyr,) 
and Bawa. = arctan(Ry,/(SFy,4/cosd)) + €, - Bcosd 


With this new value for the apparent wind angle the complete 
calculation for the sail lift and drag coefficients is repeated and 
a new value for the apparent wind angle again obtained. This 
iteration procedure is continued until convergence is reached. 


4.10 Calculation of Apparent and True Wind Speed and 
Angle 


Once convergence has been obtained for the apparent wind 
angle and the aerodynamic force coefficients, it becomes a 
simple matter to calculate the remaining entities relating to the 
performance of the yacht. With the apparent wind angle 
known, the apparent wind speed at the center of effort of the 
sails, in the plane normal to the mast, is determined from: 


Vaw = (Fresta!(204Crresta-An)) 
in which the resultant force is given by: 

Fre = (Ry? + (SFry/cos¢)?)!? 
and Crresta = (Cota? + Cyry’)!? 


The true wind speed and angle can be obtained from the ap- 
parent wind speed and angle and the boat speed, from the so- 
called velocity triangle. In the plane normal to the mast this is 
as follows: 


Vow. = CV awn + Ve - 2V awi.-Vp-COS(B awh. a Bcos¢))'? 
and Bry = arcsin(V,.sin(Baw, + Bcosd)/Viw) + Baw. 
where Vry, = true wind speed in the plane normal to the 

mast at the center of effort of the sails; 


true wind angle, in the plane normal to the 
mast. 


and Pow. 


It should be noted that the apparent and true wind angles here 
are relative to the symmetry plane of the yacht and not relative 
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to the boat speed vector, i.e. not relative to the track of the 
yacht through the water. 


It is important to determine the apparent wind speed and 
angle and the true wind speed and angle at the mast-head, so 
as to be able to correlate these with measurements made with 
the wind-sensing instruments located there. In the program this 
is executed utilizing the following relation for the wind gra- 
dient with height: 


Vawmn = Vawt-(Ho/Zce)? 


where Viy,, = true wind speed at the mast-head wind 
sensing gear; 
Hy = height of the mast-head wind-sensing gear 
above the water surface; 
ZcE = height of the resultant center of effort of 


the sails above the water surface. 


The true wind angle is assumed not to change with this change 
in wind speed up the mast. The speed-made-good to windward, 
denoted as Vyc, which is the all-important performance entity 
when sailing to windward, is then calculated from: 


Vac = Vp-cos(Brw. + Bcosd) 
4.11 Calculation of Added Resistance in Waves 


In the calculation of the added resistance in waves it is as- 
sumed that the wave direction coincides with the true wind 
direction. Accordingly, this component of resistance can only 
be considered once this wind direction is known. When the 
added resistance due to waves has been calculated, a new value 
of the total hydrodynamic resistance is obtained with which the 
complete calculation of the resultant hydrodynamic force, the 
heel angle, the sail coefficients etc, have to be repeated, 
leading to an iteration procedure with respect to the true wind 
direction and wind speed. Usually only 2 to 3 steps are re- 
quired to find convergence, yielding the total performance of 
the yacht in waves. 


AS it is not possible to include a full, so-called strip-theory 
method for the calculation of the added resistance due to waves 
in a VPP, an approximate method has to be adopted. A method 
similar to that by Gerritsma et al (1992/1993) is used for this 
purpose. 


In this method, the added resistance due to waves is ex- 
pressed in a formula involving the length L.,, the significant 
wave height H,,,, the length-displacement ratio L¢p/Vc,'3, the 
non-dimensional value of the longitudinal radius of inertia 
k,,/Lcg, and the waterline beam-length ratio By,/L¢,, in ad- 
dition to environmental conditions such as wave direction, 
wave period and wave height, fetch and wind duration over 
fetch. 


The added resistance in waves is expressed as: 
Raw = p-8-Leg-Hy3"-(a a b.(Lep/V ca). (Kyy/Lcp)-(Bwi/Les)’) 


where g = acceleration due to gravity; 
and a and b = values dependent on Froude number, wave 
direction, and wave period. 


The values for a and b were calculated for 7 of the Delft series 
I systematic models for zero heel, for a matrix of Froude 
numbers (0.15, 0.25, 0.35, 0.45 and 0.60), wave periods (1, 
1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5 and 6 seconds), wave direc- 


tions (35, 50, 65 and 80 degrees relative to the centerline of 
the yacht), and non-dimensional pitch gyradius values of 0.23, 
0.25, 0.27 and 0.31, for a Bretschneider spectrum formulation. 
These calculations were carried out for a wave height of 1 
meter only, as the added resistance is proportional to the 
square of the wave height. 


All calculations were carried out for Lp, = 10 meter which 
necessitates that for other values of the length, the wave period 
T, used for determining the values of a and b needs to be 
modified according to: 


T = T,.(g/Lcg)!” 


In the program a three-dimensional interpolation scheme is 
utilized to find the appropriate values of a and b. Also, the 
option is provided to input either a given wave height and 
wave period or to calculate the wave height and wave period. 
To use this latter option the program uses the calculated true 
wind speed and direction, and additional information on the 
length of fetch, and duration of the wind speed over the fetch, 
which have to be input. 


When the waves are aft of abeam, most yachts reveal some 
tendency to surf. Research into the speed increase due to 
surfing action in waves is presently being conducted in various 
laboratories. It is hoped that an approximate method for the 
calculation of this important speed increase, and its dependency 
on the main design variables, can soon be included in VPP 
models to render the outcome of preliminary design analyses 
more realistic in this sense. 


5 COMPUTATIONAL FLUID DYNAMICS 
PROGRAMS 


(CFD) 


5.1 Introduction 


Great strides have been made in recent years in the field of 
numerical aero- and hydrodynamics. The motivation for the 
significant development that has taken place is basically the. 


’ availability, for the first time, of adequate computing power. 


Access to super computers is now no longer a problem, while 
the run-time cost thereof is becoming acceptable to most 
organizations. 


Significant progress has in particular been made in the field 
of potential flow calculations. Here the influence of viscosity is 
neglected, requiring that the flow around the subject body is 
well-behaved since flow separation cannot be modelled. On 
neglecting viscosity the prediction of viscous resistance cannot 
be addressed. All of the properties associated with the pressure 
distribution over the body can however, leading to the use of 
this type of numerical model for the prediction of lift and 
induced resistance. 


Progress in recent years in the calculation of the effect of the 
free surface on the potential flow around ship and ship-like 
bodies is also significant. This problem is appreciably more 
difficult then the calculation of the flow in absence of the free 
surface such as in the case of a deeply submerged body. Some 
codes with this capability are now becoming commercially 
available. 


Although significant work is being carried out in the area of 
numerical viscous flow modelling, practical tools have so far 
only become commercially available for the two-dimensional 
case and simplified three-dimensional cases. This is primarily 
due to the fact that the dicretization of the flow must be pur- 
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sued with considerable finer detail then in the case of potential 
flow. Because of this problem, solution of the all-important 
Navier-Stokes equations, governing ali viscous flows, will 
require more computing power then is presently available on 
an affordable basis. In a few years, however, the availability of 
codes for the complete prediction of the three-dimensional 
viscous flow around streamlined and less streamlined bodies, 
including the calculation of separated and wake flow, should 
become available. 


The various classes of CFD tools now available to the prac- 
ticing naval architect, in the field of sailing yacht design and 
performance analysis, is as follows: 


- Programs for the calculation of the wave pattern around a 
moving, streamlined hull and the associated prediction of 
the wave resistance; 


- Programs for the calculation of the flow field around a 
lifting body and the associated prediction of the lift and 
induced resistance; 


- Programs which combine the capabilities of the above 2 
types of codes, yielding the ability to calculate the wave 
pattern, wave resistance, the lift and induced resistance and, 
for sailing yachts, the all-important coupling between the 
effect of the free-surface on the lift and induced drag. 


- Programs for the calculation of the viscous flow around 
specific streamlined, lifting and non-lifting bodies, yielding 
important insight into boundary layer behaviour. 


- Programs that predict the total, viscous flow in the presence 
of the free surface, for lifting bodies. This type of program 
is what is ultimately required to solve the various design 
problems. Although not yet available in the form of a prac- 
tical tool, development work in this area is particularly 
important and will be reviewed, together with the above- 
mentioned class of CFD tools, below. 


5.2 DAWSON-Type Models for Wave Resistance 


Although various theories for the calculation of the wave 
pattern and wave resistance have been developed over the 
years, it is the work of Dawson (1977) that has been shown to 
yield the most accurate results for steady, free-surface flow. 
Most of the effort in the development of numerical wave resis- 
tance is now carried out in further developing this model. A 
description of the state-of-the-art with respect to the underlying 
theory was recently presented by Letcher (1993). From this 
study in particular it has become clear that much work still 
remains to be carried out to perfect this approach, although 
results to date indicate a good agreement in many instances. 
Some codes utilizing this theory have also been developed to 
deal with lifting bodies. These latter codes will be described in 
Paragraph 5.4. 


Various other programs are available for the prediction of the 
wave resistance of ship and yacht hulls. Most of these utilize a 
less refined description of the hull or body shape and are 
suitable for running on smaller computers. Mention can be 
made in this respect of the programs based on the slender ship 
theory developed by Noblesse (1983). 


For sailing yachts it is essential that the program properly 
takes into account sinkage and trim. Rather then input various 
sinkage and trim values manually to study the effect thereof, 
the program should be set-up to determine and evaluate auto- 


matically the required equilibrium sinkage and trim condition 
by integration of the pressure distribution over the body. This 
necessitates that certain parts of the above-water hull form be 
input as well, in particular the bow and stern overhang, which 
become immersed at higher Froude numbers. Recent results 
obtained with this class of program reveal a good correlation 
with experimental results in a qualitative sense, although some 
bad results in some cases persist. 


As yet this class of CFD code can still not be used with the 
aim of obtaining a quantitative prediction of the wave resis- 
tance for a particular design. Although the moment in time 
when this can be done is fast approaching, these codes are 
presently best utilized by first running the program for a 
known case, for which model test results are available, and 
then to run the program for a series of similar designs for 
which the effect of different design variables has to be studied 
in a qualitative sense. On having found a particularly good 
design in this way a final check of the results through model 
testing is advisable to confirm the numerical results, particular- 
ly if the project requires accurate performance estimates. 


5.3 VSAERO-Type Models for Lift and Induced Resistance 


Numerical models for the calculation of the potential flow 
around complex three-dimensional bodies have been developed 
by workers in the field of aerodynamics intensively during the 
last 20 years. Various codes are now available that have been 
correlated extensively with wind tunnel test results yielding 
very satisfactory results for most engineering design purposes. 
Most of the codes available in the public domain are based on 
the well-known VSAERO code originally made available by 
NASA. The underlying theory of this code is that of Morino 
(1975). This theory applies the Green function to the equation 
of the velocity potential. The subject body is divided into many 
hundreds (or thousands) of combined source and doublet 
panels, depending on the complexity of the body. Once the 
source and doublet strengths of all the panels are known, the 
velocity and the pressure at each panel can be determined. The 
integration of the surface pressure over the entire surface 
yields the resultant forces and moments in each of the axis 
directions. this constitutes the so-called pressure analysis 
method. For it to yield accurate results a relatively large 
number of panels is required. The so-called far field analysis 
method, usually available as an additional option, allows for an 
accurate determination of the resultant forces on using a rela- 
tively small number of panels. Particularly the induced drag is 
determined accurately in this way. 


Some of the commercially available VSAERO codes have 
been modified extensively relative to the original NASA code. 
The effect thereof on the accuracy is often insufficiently 
documented. A noteworthy exception is the so-called DSA 
code developed from the available NASA source by Letcher 
(1989). 


Application of this type of CFD code in sailing yacht design 
and performance analysis is extremely rewarding. Complicated 
keel configurations can be analyzed providing lift, induced 
drag, streamline and pressure distribution information. A 
necessary requirement for efficient use of this type of code for 
the analysis of different designs and design variables, is that 
panel generation software be available which generates the 
panel input and panel neighbour information files. If the work 
of panel-input generation is to be carried out manually, the 
process becomes tedious and very time-consuming. 


Even though the keel and other appendages of the yacht are 
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sometimes relatively deeply submerged, there will always be 
an influence of the free surface on the lift or on the drag or 
both, particularly at relatively high values of the angle-of- 
attack (leeway angle or rotation angle of the appendages). It 
follows that this type of code, in not being able to take the free 
surface into account, will always be of limited use only. Care 
needs to be taken in performing systematic studies in which 
small variations in geometry are investigated, to not arrive at 
incorrect conclusions. The effect of the free surface in some 
instances can be such so as to yield a completely reversed 
trend to that found when assuming the free surface has no 
influence. Accordingly, this tool is best used when studying 
new concepts or a wide range of geometry parameters. An 
example of the way this code is best used is the case of the 
development of the keel for "Australia II”, in which widely- 
varying keel geometries were studied (Van Oossanen, 1985). 


5.4 HYDROPAN-Type Models for Wave Resistance, Lift 
and induced Resistance 


A new class of CFD code has recently become available. 
This code combines the capability of the Dawson type of 
theory with that of the Morino combined source and doublet 
panels distributed on the surface of the body and the free 
surface. This was first carried out by Van Beek et al (1985) at 
the Netherlands Aerospace Laboratory, leading to the so-called 
HYDROPAN code, which work has since been continued at 
MARIN in the Netherlands. A noteworthy development along 
similar lines has been carried out by Boppe et al (1987), 
leading to the code called SPLASH in the USA (Rosen et al, 
1993). Both codes have produced very promising results with 
respect to wave resistance, the lift on the appendages, and the 
associated induced resistance, for yacht-like configurations. 
More importantly, the effect of the free surface, as predicted 
by these codes, is as found from model tests. 


Up till now, these codes have also not yet been developed to 
the stage where the results obtained can always be used in a 
quantitative sense. Their value for the time being lies in the 
qualitative prediction of performance and the associated ability 
to carry out exploratory analysis. It is this program in par- 
ticular that most serious workers in the field of yacht design 
and performance analysis will be wanting to use in the near 
future. 


5.5 Viscous Flow Models 


Programs for the evaluation of the two-dimensional boundary 
layer behaviour of the sections used for keels, winglets and 
tudders have been available for some time. The estimation of 
laminar-turbulent transition behaviour, the point where flow 
separation occurs based on pressure distribution information 
derived from potential flow computations, etc, are common- 
place and used by most naval architects with good results. On 
the contrary, the use of available codes for the estimation of 
simplified three-dimensional viscous flow on less complicated 
bodies such as the canoe body of sailing yachts, a ballast bulb 
fitted to the tip of the keel, etc, is as yet not carried out to any 
great extent. This is only partly a result of the fact that only a 
few codes are commercially available (and then only at a 
relatively high cost) for this task. In part this is also because 
the interest of the designer of sailing yachts is mainly focussed 
on the design of the appendages and the wave resistance of the 
canoe body. The possible gains in reducing the viscous resis- 
tance to a minimum in sailing yachts are sometimes only small 
in comparison to the gains that are possible on minimizing the 
wave resistance, for example. This is likely to remain so for 
some time, particularly since the move away from the TOR (the 


International Offshore Rule) to IMS and other yacht classes, 
yielding a finer, fairer canoe body shape, requires even less 
capability to design an efficient hull form from a viscous flow 
point of view, free of any flow separation problem. The only 
significant requirement for three-dimensional viscous flow 
codes at the present time would seem to lie in the need of 
designing efficient bulb shapes. 


5.6 Seakeeping Models 


Numerical codes for the calculation of the motions and added 
resistance of ships and sailing yachts with complicated ap- 
pendage configurations are also becoming available. A code 
that should be mentioned in this regard, recently used for the 
performance analysis of LACC yachts, is that developed at the 
Massachusetts Institute of technology by Sclavounos and Nakos 
(1993). A panel method is utilized as in other codes. Although 
the published results are for the canoe body only, there is no 
reason why this code cannot be utilized for the complete 
configuration at heel and leeway. It is this capability of such a 
code that will generate the most interest amongst workers in 
the field of sailing yachts since no other type of program is 
available that will allow the study of such a complex con- 


' figuration at the combinations of heel and leeway of interest, in 


waves. 
5.7 Total Flow Models 


Considerable effort is presently being placed on the develop- 
ment of so-called total flow codes for the complete calculation 
of the viscous flow around hull forms with lifting appendages, 
including effects associated with the free surface. To this end 
various research organizations have coupled the various in- 
dividual codes and organized the resulting software to run in 
an iterative way, to deal with the problem of first having to 
know the pressure distribution on the body before the viscous 
flow can be calculated. Insufficient information is as yet avai- 
lable on the required computing power and cost of such soft- 
ware, although it is known that at least one such software 
system is commercially available. 


6 FINAL REMARKS 


Within the confines of time and space afforded in the writing 
of this paper, an attempt has been made to address the topic of 
speed prediction of sailing yachts in a concise, yet all-encom- 
passing manner. The three different methods available for this 
task have been described, viz: model testing, the use of semi- 
empirical VPP models, and computational fluid dynamics. At 
this stage none of these methods alone can suffice in any major 
project involving the optimization of speed potential. Model 
test results need to be analyzed with the aerodynamic part of a 
VPP model to enable extrapolation of the measured perfor- 
mance in determining boat spéed at the various points of sail. 
Velocity prediction programs, required to ascertain the impact 
on performance of overall hull and sailplan design variables, 
need constant improvement, for which experimental and CFD 
analyses are necessary. Computational fluid dynamics, in not 
yet being able to calculate boat speed at the various points of 
sail, requires the use of VPP models to manipulate the detailed 
results for parts of the performance equation such as appen- 
dage lift, induced drag, wave resistance, added resistance in 
waves, etc. 


In a way the "past - present - future" theme of this SNAME 
Centennial Annual Meeting is well served by the exposition of 
the three methods outlined above for speed prediction. The 
method based on model testing is representative of the past. 
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The method based on VPP modeling represents the present, 
while the promise of CFD codes represents the future. This is 
not to uphold the assertion of some that model testing is now 
redundant, or that VPP models represent a better methodology 
than model testing. At this time both methods are complimen- 
tary as stated more than once in the body of this paper. 


As already hinted at in the introduction (Paragraph 1), this 
paper is presented in the hope that the various workers in the 
field will become more conscious of the fact that a lack of 
knowledge about the details of sailing yacht performance still 
exists. This paper has hopefully illustrated that it is incorrect to 
always assume, for example, that the total hydrodynamic side 
force is a linear function of leeway angle, that the total resis- 
tance due to side force is a function of the square of that side 
force, that at large angles of heel the side force is zero for 
zero leeway angle, etc. These aspects need to be studied from 
case to case. 
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APPENDICES 


Appendix 1 Example of the Analysis of Upright Resistance 
Test Results with a Model of the Canoe Body 
Only, Using the Prohaska Method 


A 1:3 scale model of the canoe body of the 12 Meter yacht 
“Australia II", winner of the America’s Cup in 1983, was 
tested in 1982 for the purpose of determining the form factor 
of the canoe body and the base resistance thereof, without the 
effects associated with wave-making of the appendages playing 
a role. The model was altered slightly to obtain the wanted 
flow in absence of the rudder. No turbulence stimulation 
means were adopted in these tests because of the large size of 
this model, the length on the test waterline being 4.677 meters. 


The calculation of Cerca Cream and F,‘/Crcpm), a8 des- 
cribed in Paragraph 3.5.1, is given in Table 5 for each of the 
19 speeds comprising this particular test. 


The results given in Figure 9 and Table 5 clearly reveal lami- 
nar flow to be present on the model at speeds up to about 0.9 
m/sec (for which Ry = 4.08 x 10°). The value of the form 
factor was determined to be 1+Kcg = 1.118, a representative 
value for good 12 Meter yachts. 


Table 5. Results of upright resistance measurements and cal- 
culations of Cercpam/Crcagy and Fy/Cecgmy for 1:3 scale model 
of "Australia II” without appendages. 





Vem) Fy Ceca) Cerceim/Crcsem 
TCB(M) Ryx107 RTCB(M) FCB(M) 
0.300 0.41 0.044 0.136 4.389 1.808 0.412 0.001 
0.450 1.57 0.066 0.204 4.038 3.062 0.758 0.005 
0.603 3.36 0.089 0.274 3.810 3.645 0.957 0.017 
0.751 5.60 0.111 0.341 3.651 3.923 1.075 0.041 
0.901 8.19 0.133 0.409 3.527 3.986 1.130 0.089 
1.050 10.91 0.155 0.476 3.427 3.908 1.140 0.169 
1.200 14.19 0.177 0.544 3.344 3.891 1.164 0.295 
1.350 17.81 0.199 0.612 3.273 3.861 1.180 0.482 
1.498 22.13 0.221 0.680 3.212 3.896 1.213 0.745 
1.652 27.62 0.244 0.749 3.156 3.998 1.267 1.121 
1.800 33.79 0.266 0.817 3.109 4.119 1.325 1.604 
1.953 41.66 0.288 0.886 3.064 4.315 1.408 2.255 
2.101 51.31 0.310 0.953 3.025 4.591 1.518 3.060 
2.251 62.88 0.332 1.021 2.989 4.902 1.640 4.080 
2.401 79.67 0.355 1.089 2.956 5.459 1.847 5.340 
2.550 112.54 0.377 1.157 2.926 6.836 2.337 6.865 
2.701 165.42 0.399 1.225 2.897 8.956 3.092 8.727 
2.852 236.20 0.421 1.294 2.870 11.470 3.996 10.949 
3.000 324.76 0.443 1.361 2.846 14.253 5.009 13.521 
Vem) = Speed of model in m/sec; 
+cam) = Jotal resistance of model in Newton; 
Fy = Froude number based on Logm) (bcamy = 4-677 
m);. 
Ry = Reynolds number based on Vay) and oe with 
Ym = 1.031 x 10°m?*/sec (for 19 °C); 
Crcam) = Friction coefficient of canoe body (x 10°); 
Cerca) = Real AP nV pa scan) with pay) = 1000 kg/m? 
and Scam) = 





Appendix 2 Examples of Prohaska Analyses of Upright 
Resistance Measurements on Model Fitted with 
Keel, Winglets and Rudder 


This example considers the same model as in Appendix 1, but 
now fitted with keel, winglets and rudder. Although the 
Prohaska analysis is not required when the results are available 
of tests with the canoe body alone, the technique of subtracting 
the viscous resistance of appendages to derive the value of the 
form factor for the canoe body in that case provides important 
information on the nature of the viscous resistance of those 
appendages. Specifically, it is possible to determine the degree 
of laminar flow present on keel or rudder or both, with speed, 
by comparing the measured resistance values of the canoe body 
to those obtained by subtracting the assumed viscous resistance 
of the keel and rudder from the total measured resistance. In 
this way various formulations for appendage resistance can be 
studied, yielding considerable insight in the build-up of the 
resistance of those appendages. 


In the case of tests with large models, in which no turbulent 
stimulation means are adopted, a good correlation with the 
results of upright resistance tests with the canoe body alone is 
usually obtained when the viscous resistance of the keel (and 
winglets) is based on a Jaminar-turbulent transition formulation 
for the value of the friction coefficient, and a fully turbulent 
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flow is assumed for the rudder. This result indicates that the 
rudder, when situated exactly downstream of the keel such as 
in the upright resistance test, probably experiences fully tur- 
bulent flow in most cases. In addition to these corrections, an 
account of the viscous resistance associated with boundary- 
layer interaction between canoe body and keel (and between 
canoe body and rudder-skeg arrangement) usually needs to be 
made as well. This will now be demonstrated. 


A formulation often used by the author for modelling the 
value of the (flat plate) friction coefficient in the laminar-tur- 
bulent transition region, is based on the formulation by Prandtl 
(1963), with the coefficient 1800 instead of 1700, to model 
transition to commence at Ry = 5 x 10° and terminate at Ry = 
5 x 10%, viz: 


C, = 0.075/(Log(R,) - 2)? - 1800/Ry 


This has been used in Table 6 for the calculation of Ryxuy (and 
Rvwim). AS stated above, the calculation of Rypay) usually 
needs to be based on a completely turbulent value for Crp. 


A formula for the approximate viscous resistance associated 
with canoe body-keel boundary layer interaction has been 
developed by the author from a basic expression given by 
Hoerner (1965). This formula takes into account the radius of 
the fillet at the junction of the canoe body and the keel, in 
addition to the thickness parameters identified by Hoerner, viz: 


Rintk = 2p Vp' Atyx)’O.75(t/c),x)-(1 - 8((r7/C,)x 


+ 422(rc,))* - (Uc) x’) 


where t = maximum thickness of keel, without fillets, 
at canoe body-keel interface; 
(t/c),, = thickness-chord length ratio at the canoe 


body-keel interface, without fillets; 

= minimum radius of keel fillet in transverse 
plane, divided by chord length at root of 
keel. 


and (r/c), 


When the canoe body has considerable deadrise, as in this 
case, the angle between the surface of the canoe body and the 
keel, as seen in the transverse plane, is significantly greater 
than 90 degrees. The value of R,y; should then be multiplied 
by the sin of that angle to approximately take into account the 
lesser extent of the boundary layer interaction in that case. 


Table 6 gives all of the required calculations in this case. The 
resulting Prohaska plot is also shown in Figure 9. 


Appendix 3 Example of Analysis of Model Test Results with 
Respect to Side Force 


The results of the side force measurements on the 1:7 scale 
model of the IACC yacht shown in Figure 2 are here analyzed 
according to the method outlined in Paragraph 3.5.4. The 
measurement values are first of all given and re-worked to 
yield lift coefficients as given in Table 7. 


The result of applying regression analysis to the measured 
data, using the mathematical model given in Paragraph 3.5.4 
yielded the following result: 


(Cyrmaw)man = ~ 0-07092¢? + 1.649(1 - 0.2321¢7).(Bcosd) 


+ 0.4534(Bcosd)* - 39.31¢6.Fy(Bcose)* 


The accuracy of this simple mathematical model for the total 
lift coefficient can be investigated by comparing the values of 
(Coropexp With (Cityemy)man in Table 7. 





Table 6. Results of upright resistance measurements and cal- 
culations of Cyay/Crcay and Fy*/Cecgy for 1:3 model of 
“Australia II” fitted with keel, winglets and rudder. 








Rvappawy Rr Ryxl0” Cray Cran/Crcean 
Vow Riwtkimy Fy Cream Fy‘/Crcpomy 
0.297 0.30 0.02 1.00 0.04 0.14 4.39 4.478 1.02 0.001 
0.497 0.66 0.05 2.68 0.07 0.23 3.95 4.286 1.09 0.007 
0.697 1.37 0.09 5.03 0.10 0.32 3.70 4.090 1.11 0.029 
0.898 2.51 0.15 8.03 0.13 0.41 3.52 3.934 1.12 0.085 
1.100 3.96 0.23 11.75 0.16 0.50 3.39 3.836 1.13 0.198 
1.299 5.68 0.32 16.31 0.19 0.60 3.29 3.818 1.16 0,397 
1.503 7.71 0.42 22.35 0.22 0.69 3.20 3.908 1.22 0.730 
1.652 9.38 0.51 27.03 0.24 0.76 3.15 3.912 1.24 1.083 
1.803 11.21 0.61 33.05 0.26 0.83 . 3.10 4.016 1.30 1.561 
1.951 13.16 0.72 41.47 0.29 0.89 3.06 4.306 1.41 2.170 
1.999 13.82 0.75 44.24 0.29 0.92 3.05 4.373 1.44 2.401 
2.102 15.29 0.83 50.61 0.31 0.96 3.02 4.525 1.50 2.96 
2.250 16.77 0.95 61.77 0.33 1.03 2.98 4.820 1.62 3.935 
2.400 19.11 1.08 82.22 0.35 1.10 2.95 5.639 1.91 5.151 
2.547 22.39 1.22 118.00 0.37 1.17 2.92 7.185 2.46 6.601 
2.701 25.13 1.37 179.40 0.40 1.24 2.89 9.714 3.36 8.432 
2.904 28.96 1.59 290.73 0.43 1.33 2.86 13.618 4.77 11.407 
3.101 32.92 1.81 432.58 0.45 1.42 2.83 17.770 6.29 14.997 


3.254 36.15 1.99 552.72 0.48 1.49 2.80 20.620 7.36 18.331 





Vem) = Speed of model in m/sec; 
vant) = Ryxamy + Rerow + Ryo; ‘ 
vey = Pay Vaaay”-Cexny(l +kx)-Skw in Newton; 
Cry = 0.075/(Log (Ry) - 2)? - 1800/Raxomys 
1+ky = 1.142 (based on (t/c), = 0.11); 
Soin = 1.036 m’; 
Ryamy = 20 cayVaemy’ Creal +Kp) Spam in Newton; 
Cerny = _:(0.075/(Log(Ryg) - 2); 
1+k, = 1.159 (based on (t/c), = 0.12); 
Sam) = 0.147 m’; 
vw) = 2.4 Vee’-Cewa(l tkw)-Swaa in Newton; 
rw) =: 0.075/(Log(Rywa) - 2)? - 1800/Rywey: 
1+k x = 1.077 (based on (t/c), = 0.063); 
Swim = 0.328 m’; 
intxm) = 0.000188, —4)Vaany’ in Newton, based on tum) = 
0.093 m; 
(tlc)x = 0.0918; 
(r/c), = 0.0908; 
Rum = Total resistance of model as measured, minus 
Ryximy> Rywaays Ryrewy and Riyrxay in Newton. 
(This can only be identified with the resistance of 
the canoe body at low values of the Froude num- 
ber since the appendages also posses a wave- 
making component); 
Fy = Froude number based on Leg (Lea) = 4.763 
m); 
Ry = Reynolds number based on Vay and Leg, with . 
Ym = 1.039 x 10*m*/sec (for 18.7°C); 
Crcam) = Friction coefficient of canoe body (x 10°); 
Crm) = Rywy/(20aVee-Scpo)> With pay, = 1000 kg/m? 
(x 10°); 
Scam =: 5.063 m’. 
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Table 7. Calculation of coefficient of total lift derived from 
results of side force measurements on 1:7 scale model of an 
International America’s Cup Class yacht shown in Figure 2. 


Vom) Bcos YapoyVey-Arr-Cose — (Cray) math 
TH(M) LTH(M)/exp 
14 10 0.04211 24.56 377.1 0.06513 0.06745 
14 10 0.02183 12.34 377.1 0.03272 0.03378 
1.4 10 0.06875 42.11 377.1 0.11170 0.11194 
1.6 10 0.05277 41.75 492.5 0.08477 0.08519 
1.6 10 0.03334 27.06 492.5 0.05494 0.05285 
1.6 10 0.01513 11.30 492.5 0.02294 0.02271 
1.8 10 0.01289 12.08 623.4 0.01938 0.01902 
1.8 10 0.02733 26.52 623.4 0.04251 0.04288 
1.8 10 0.04074 40.97 623.4 0.06572 0.06513 
2.0 10 0.03524 42.61 - 769.6 0.05537 0.05599 
2.0 10 0.02131 25.35 769.6 0.03294 0.03291 
2.0 10 0.01117 13.29 769.6 0.01727 0.01618 
1.4 20 0.09414 51.38 359.8 0.14280 0.14310 
1.4 20 0.07134 37.77 359.8 0.10500 0.10662 
1.4 20 0.11891 65.96 359.8 0.18330 0.18202 
1.6 20 0.09266 64.36 469.9 0.13700 0.14032 
1.6 20 0.06954 47.89 469.9 0.10190 0.10354 
1.6 20 0.04920 33.96 - 469.9 0.07227 0.07078 
1.8 20 0.03936 32.52 594.8 0.05467 0.05484 
1.8 20 0.05363 47.01 594.8 0.07903 0.07786 
1.8 20 0.06954 62.71 $94.8 0.10540 0.10336 
2.0 20 0.05871 63.16 734.3 0.08601 0.08591 
2.0 20 0.04494 49.79 734.3 0.06781 0.06380 
2.0 20 0.03083 33.47 734.3 0.04558 0.04103 
1.4 30 0.15420 69.50 331.6 0.20960 0.20885 
1.4 30 0.08147 35.19 331.6 0.10610 0.10633 
1.4 30 0.21630 92.06 331.6 0.27760 0.27895 
1.6 30 0.15890 92.83 433.1 0.21430 0.21162 
1.6 30 0.10720 62.88 433.1 0.14520 0.14339 
1.6 30 0.06484 34.43 433.1 0.07961 0.08083 
1.8 30 0.05336 34.02 548.2 0.06206 0.06314 
1.8 30 0.08812 61.64 548.2 0.11240 0.11520 
1.8 30 0.12360 90.36 548.2 0.16480 0.16470 
2.00 30 0.10130 90.07 676.8 0.13310 0.13329 
2.0 30 0.07316 61.77 676.8 0.09127 0.09281 
20 30 0.04444 33.26 676.8 0.04914 0.04937 


Vay = Speed of model in m/sec; 
= Heel angle in degrees; 
Bcos@ = Effective angle-of-attack in radians; 
SF = Total measured side force in Newton; 
Pm = Density of fresh water in towing tank (1000 kg/m’); 
Aura = Total lateral area of model (0.3907 m’); 
tram) = Coefficient of total hydrodynamic lift on model. 


The subscript "exp" refers to the measured value 
while the subscript "math" refers to the value ob- 
tained from the simple mathematical model given in 
this Appendix. 


Note that the coefficient of linear lift is given by the expres- 
sion 1.649(1 - 0.2321¢7).(@cos@), the non-linear lift coefficient 
by 0.4534(8cos¢)’ and the loss in lift due to flow separation by 
39.31.F,(@cos@)’. As can be seen, a coupling exists between 
Froude number and heel angle with respect to this last com- 
ponent. No significant physical reasoning can be based on this 
result since the experimental data covers too small a region 


with respect to Froude number, heel angle and angle-of-attack. 
On setting the left-hand side of this equation equal to zero and 
solving for Bcos@, the value of the zero lift angle of attack is 
found. By inspection, it can be seen that this value is 
dominated by a ¢? term. 


The effective value of the aspect ratio of the model as a 
whole follows from (see Paragraph 3.5.4): 


AR, = 82(8Cyyry/8a),/(422 - (8Cyyry/Aa),2) 
with (8C,.ru/8a), = 1.649(1 - 0.23214”) 
The results of this analysis are shown in Figure 11. 


Appendix 4 Example of Analysis of Model Test Results with 
respect to Resistance Due to Heel, Induced 
Resistance and Remaining Resistance due to 
(Lift) Attitude 


The results of the resistance measurements on the 1:7 scale 
model of the IACC yacht, shown in Figure 6, are here ana- 
lyzed according to the method outlined in Paragraph 3.5.5, 
3.5.6 and 3.5.7. The measurement values are first of all given 
and re-worked to yield resistance coefficients as given in Table 
8. 


On applying regression analysis to the measured data, using 
the mathematical model given in Paragraph 3.5.6, the fol- 
lowing result was obtained: 


(Cerna ~ Crrnoaye=o)man = 9-002319¢ + 0.3501¢?.Fy* 
+ (0.8141 + 5.470Fx2)(Croruan)man)/(TAR,) 
+ 74.9(Bcos¢ - 0.15)? 
where (CyrHaw)man = ~ 0.070929?+ 1.649(1-0.2321 47). (Bcos¢) 
The expression for AR, is given in Appendix 3. 


Note that in this case the resistance due to heel has a linear 
term in heel angle and a term dependent on Froude number to 
the 4th power. This latter term is clearly associated with the 
increase in wave resistance as the model heels. Note also that 
no terms appear in this mathematical model dependent on the 
3rd power of the angle-of-attack other than the C,,, term, 
which term was not obtained by regression analysis but by 
inspection, using the available data points displaying the exis- 
tence of this component of resistance. This Cp,, term is only 
utilized for angles-of-attack greater than 0.15 radians. 


The lack of another term in the angle-of-attack to the 3rd 
power indicates that the induced resistance associated with the 
non-linear lift, as identified in Appendix 3, is part of the 
induced resistance as depicted by the term dependent on the 
square of the (linear) lift coefficient. This, in turn, indicates 
that the non-linear lift occurs near the “leading edge” of the 
low aspect ratio components (the canoe body and the large 
bulb fitted to the keel), and is not due to separation from the 
lateral edges of these bodies, as discussed in Paragraph 
2.3.6.2. 


Finally, it should be noted that the coefficient of the induced 
resistance term, dependent on the square of the (linear) lift 
coefficient, can be directly compared to the efficiency factor 
1+o of the model in producing lift. In Paragraph 2.3.6.1 it is 
explained how this factor accounts for the effect of the non- 
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Table 8. Calculation of coefficient of total resistance associated 
with heel and leeway from measurements on 1:7 scale model 
of an International America’s Cup Class yacht (see Figure 6). 


Vem) Rmiay-Rrane=0 (Carnem7Crrncme=o)exp 


YapanV aan? Att 


Bcos¢ (Cerne -Crrnimye=0)math 
1.4 10 0.04211 0.62 382.9 0.00162 0.00201 
14 10 0.02183 0.31 382.9 0.00081 0.00086 
1.4 10 0.06875 1.69 382.9 0.00441 0.00470 
1.6 10 0.05277. 1.56 500.1 0.00312 0.00322 
1.6 10 0.03334 0.84 500.1 0.00168 . 0.00156 
1.6 10 0.01513 0.18 500.1 0.00036 0.00070 
1.8 10 0.01289 0.36 633.0 0.00057 0.00072 
1.8 10 0.02733 0.85 633.0 0.00134 0.00133 
1.8 10 0.04074 1.42 633.0 0.00224 0.00233 
2.0 10 0.03524 1.76 781.4 0.00225 0.00209 
2.0 10 0.02131 1.17 781.4 0.00150 0.00115 
2.0 10 0.01117 0.61 781.4 0.00078 0.00077 
1.4 20 0.09414 3.05 382.9 0.00797 0.00825 
1.4 20 0.07134 1.91 382.9 0.00499 0.00502 
1.4 20 0.11891 4.73 382.9 0.01235 0.01283 
1.6 20 0.09266 4.20 500.1 0.00840 0.00889 
1.6 20 0.06954 2.59 500.1 0.00518 0.00536 
1.6 20 0.04920 1.35 500.1 0.00270 0.00315 
1.8 20 0.03936 2.14 633.0 0.00338 0.00274 
1.8 20 0.05363 2.92 633.0 0.00461 0.00405 
1.8 20 0.06954 4.07 633.0 0.00643 0.00604 
2.0 20 0.05871 4.34 781.4 0.00555 0.00529 
2.0 20 0.04494 3.13 781.4 0.00401 0.00372 
2.0 20 0.03083 2.11 781.4 0.00270 0.00259 
1.4 30 0.15420 7.67 382.9 0.02003 0.01951 
1.4 30 0.08147 2.47 382.9 0.00645 0.00595 
1.4 30 0.21630 23.09 382.9 0.06030 0.06031 
1.6 30 0.15890 12.07 500.1 0.02414 0.02307 
1.6 30 0.10720 5.63 500.1 0.01126 0.01086 
1.6 30 0.06484 2.28 500.1 0.00456 0.00478 
1.8 30 0.05336 2.75 633.0 0.00434 0.00446 
1.8 30 0.08812 5.43 633.0 0.00858 0.00882 
1.8 30 0.12360 10.03 633.0 0.01585 0.01596 
2.0 30 0.10130 9.92 781.4 0.01270 0.01281 
2.0 30 0.07316 6.08 781.4 0.00778 0.00780 
2.0 30 0.04444 3.36 781.4 0.00430 0.00464 


Va) = Speed of model in m/sec; 
= Heel angle in degrees; 
Bcos¢ = Effective angle-of-attack in radians; 
Rina = Total measured resistance in Newton; 
THM)s=0 = Total measured resistance at zero heel, in New- 

ton; 

Pm) = Density of fresh water in towing tank (1000 
kg/m’); 

Arm) = Total lateral area of model (0.3907 m?); 

CerHm) = Coefficient of total hydrodynamic resistance of 
model; 

Crrame-o = Coefficient of total hydrodynamic resistance of 
model at zero heel. 

Note: The subscript "exp" refers to the experimental value 


while the subscript "math" refers to the value ob- 
tained from the mathematical model given in Appen- 
dix 4. 


Table 9. Calculation of the lift coefficient for a yawed, sur- 
face-piercing hydrofoil studied by Kuhn and Scragg (1993) for 
the full span. 


Fare AR, OC, /de CLactip) CLa2 Ciiaas) Crea) Cue=2) 


0.50 2.000 6.824 4.398 0.96 0.1170.586 0.292 0.145 
1.01 1.412 4.818 3.940 0.95 0.166 0.524 0.261 0.130 
: 1.122 3.828 3.587 0.94 0.209 0.476 0.237 0.118 
2.02 1.052 3.586 3.481 0.94 0.223 0.462 0.230 0.115 
1.026 3.502 3.441 0.94 0.228 0.457 0.227 0.113 





Speed of hydrofoil in knots; 

Froude number based on chord Jength c of hydro- 

foil (0.425 m); 

Fare = Factor for the calculation of effective aspect ratio 
when top of lifting surface is closer to the surface 
than c/4, equal to (1 + 0.422/F,3); 

- Effective aspect ratio equal to F,p,.AR, where AR 

is the geometric aspect ratio = 1.45/0.425= 3.412; 

Qmk,.ARJ/(2k, + cosA(ARZ/(cosA)* + 4)!” ); 

1+0.82(t/c)tanz,(0.117/(t/c) +3.2(t/c) +3.9(t/c)’) 

which, for vc = 0.09 and tanr, = 0.0963, efuals 

0.9178; 

Sweep angle of quarter-chord line (0 degrees); 

Linear-lift factor for the round tip is: 

Croitipy) = 1 - 0.135/AR,2?; 

Coefficient of the non-linear lift term for the round 

tip: Cyao = (Craagipy-TR + O-1)/AR,, with Chardin) = 

0.70. Since TR=1, it follows that c,,. = 0.8/AR,; 

= Lift coefficient for 8 degrees angle-of-attack is 

(Cre ttipBC/Ba)e-Lep far + Craz-@”)Cosd. Here the 

heel angle is zero and the factor accounting for the 

interaction with the canoe body equals 1.0, i.e. C, 

= Crain (OC /Oa)a + Cyy2-07. 
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elliptical lift distribution on the induced resistance. The value 
of this factor is given by: 


1+o = 0.8141 + 5.470F,?. 


For the lowest value of the Froude number here considered, 
the value of 1+o = 1.221. For the highest value of the Froude 
number this value is 1.644. This result indicates that the distri- 
bution of lift over the canoe body, keel fin, bulb and rudder as 
a whole is inefficient and that a design approach that only 
allows the high aspect ratio appendages to produce lift prob- 
ably constitutes a more efficient method from the viewpoint of 
induced resistance. This is one of the reasons behind the 
various concepts recently developed for various America’s Cup 
yachts, utilizing rotating appendages (Van Oossanen, 1992). 


The results of this analysis are also shown in Figure 12. 


Appendix 5 Calculation of Lift on a Yawed, Surface-Pier- 
cing Hydrofoil (Kuhn and Scragg, 1993) 


The subject hydrofoil was tested at a full span, submerged 
depth of 1.45 m and at a submerged depth of 0.76 m. This last 
case is referred to as the half span case. The chord length of 
0.425 m was constant over the height. The foil section adopted 
was a NACA 63A-009, which was also held constant over the 
height. The tip of the hydrofoil was rounded. 
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Table 10. Calculation of the lift coefficient for a yawed, sur- 
face-piercing hydrofoil studied by Kuhn and Scragg (1993) for 
the half span. 


VioFy Fare AR, OC,/8Q@ Chattipy Star Crees) Creasy Crtex) 


2 0.50 2.000 3.577 3.476 0.94 0.224 0.462 0.230 0.115 
4 1.01 1.412 2.525 2.880 0.93 0.317 0.379 0.188 0.094 
6 1.51 1.122 2.007 2.478 0.92 0.399 0.324 0.160 0.080 
8 2.02 1.052 1.880 2.367 0.91 0.426 0.310 0.153 0.076 
10 2.52 1.026 1.835 2.326 0.91 0.436 0.304 0.150 0.074 


Speed of hydrofoil in knots; 

Froude number based on chord .length c of hydro- 
foil (0.425 m); 

Fyne = Factor for the calculation of aspect ratio when top 
of lifting surface is closer to the surface than c/4, 
equal to (1 + 0.422/F,°); ; 


oi 
z 
tou 


AR, = Effective aspect ratio equal to F,,..AR, where AR 
is the geometric aspect ratio = 0.76/0.425 =1.788; 

OC Jaa = 2rks.ARJ/(2k, + cosA(AR,2/(cosA)* + 4)!” ); 

ks = ]+0.82(t/c)-tan7,(0.117t/c)+3.2(t/c) +3.9(t/c)*) 
which, for t/c = 0.09 and tanz, = 0.0963, equals 
0.9178; 

A = Sweep angle of quarter-chord line (0 degrees); 

CLacip) = Linear-lift factor for the round tip is: 
Cuaitip) = 1 - 0.135/AR2%; 

Ciao = Coefficient of the non-linear lift term for the round 


tip: Coo = (Croagipys TR + O.1/AR,, with Cyarcip) = 
0.70. Since TR=1 it follows that c,,. = 0.8/AR,; 

Cua- = Lift coefficient for 8 degrees angle-of-attack is 
(Craitip( OC /a)a.Lep far + Cro2-0?)cosd. Here the 
heel angle is zero and the factor accounting for the 
interaction with the canoe body equals 1.0, i.e. C, 
= Critip(OC Aaa + Cy. 


Tests were carried out with the hydrofoil piercing the water 
surface at 2, 4, 6, 8 and 10 knots, in the towing tank of the 
David Taylor Ship Research Center at Washington. The an- 
gles-of-attack studied were 2, 4 and 8 degrees. This latter 
angle-of-attack was only run at the lower speeds. 


This well documented study provides an excellent opportunity 
to test mathematical formulations such as discussed in this 
paper for the calculation of the three-dimensional lift of ap- 
pendages, for which the effect of the free surface on the lift 
characteristics needs to be addressed. Table 9 gives the results 
of the calculations carried out for the full span case and Table 
10 for the half-span case. 


Figure 16 shows the comparison between these results and the 
measurement values taken from the given reference. With the 
exception of 2 data points (for the 2 knot test runs at 4 and 8 
degrees angle-of-attack, for the half span and full span case 
respectively), the correlation is good. 


{Discussion and Closure follow] 
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Discussion 


Siu C. Fung, Member 


{The views expressed herein are the opinions of the discusser and not 
necessarily those of the Department of Defense or the Department of 
the Navy.] 


I wish to congratulate the author on a job well done. To 
me this paper could be considered as a mini-hydrodynamic 
textbook on sailing yacht designs rather than just an excellent 
paper. The only area 1 want to comment upon is the statistical 
approach on wave resistance which is mainly based on Gerrits- 
mas’ Delft Systematic Yacht Series. Since the author endorses 
the idea of using Gerritsmas’ DUT polynomials for wave resist- 
ance calculations, J wish to have his response on the following 
comments. 

1. A “minus” sign seems to be missing from the coefficient 
{a;) at Fy of 0.175; see Table 2. I believe the correct coefficient 
should read —31.914. 

2. The general rule of thumb in the relationship between 
the number of data points and the number of terms in a 
multiple regression equation is as follows: 


NP = NC*(NC + 3)/2 (1) 


where NP. is the number of data points required and NC the 
number of terms of the polynomial regression equation. 

The number of regression coefficients ranges from 5 to 9 
terms depending on the speed regime. That means the re- 
quired number of data points for speed range at Froude num- 
ber of 0.45. and below should be 54, and above that 20 data 
points are required. However, the data points used for statisti- 
cal analyses were 39 and 15 (down to 12 for Fy = 0.725 —- 
0.75), respectively. The applicability of a polynomial equation 
that is generated by excess amount of independent variables 
(or number of terms) than that required by equation (1) is 
always subject to question. By doing so, high payback values 
will be generated but the equation itself may fail to perform 
as a good predictor. 

3. Unfortunately, Gerritsma did not provide any guidance 
as to the range of values of the independent parameters for 
which the DUT polynomials are valid: The author postulated 
that “this can be approximately deduced from Gerritsmas’ 
Delft Series.” I would like to add the following to his 
statement: ‘ 

(a) Caution should be taken if any of the individual hull 
form parameters is larger or smaller than the standard devia- 
tion ‘from the mean value. 

(b) Regression equations for predictions should be avoided 
if the target ship possesses any hull form parameter which is 
close to the edge or outside the boundary of the database. 

4. In general, Gerritsma’s polynomials are able to provide 
a smooth speed-power curve, particularly when they are back- 
fitted to the original database in the low-speed regime. How- 
ever, occasionally discontinuities are exhibited in the high- 
speed regime. A case example can be seen in Fig. 19 (herewith) 
when the: equations are used to backfit Delft Series Model 
Number 22. Discontinuities exhibited in the predicted speed- 
power curves from speed-independent.regressions are very 
common and are usually caused by the following: 

{a) The statistical distribution of the hull form. characteris- 
tics which make up the resistance database is not homoge- 
neous throughout the speed range. The small database that 
is usually composed by systematic hull form series is a classic 
example. 

(b) The number of data observations varies with speed as 
does the influence of the hull form parameters on resistance. 
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Fig. 19 Gerritsmas’ DUT polynomials 


By using different sets of hull form parameters for different 
speed regimes in speed-independent multiple-regression will 
only magnify the probability of discontinuity occurrence. The 
speed-power curve predicted by Gerritsma’s DUT polynomi- 
als for the “Sample” yacht in Fig. 19 is a good example [the 
only differences between Model 22 and the “Sample” are 
pigs coefficient (C,) and longitudinal center of buoyancy 
(LCB)]. 

5. Longitudinal center of buoyancy seems to play a major 
factor on the discontinuity of the DUT polynomials predic- 
tions. In reality, LCB has strong ties to both section area 
shape and afterbody design, e.g., keel rise and transom. Unfor- 
tunately the Delft Systematic Yacht Series only consists 
of hull forms with dry (or emerged) transoms and were mainly 
focused on resistance at even keel. Otherwise, transom sub- 
mergence (e.g., transom area, transom width, buttock angle 
...) could be a set of useful parameters for statistical analysis. 
Because people do trim their yachts when surfing, so that the 
transom may increase the effective ship length, I wish the 
yacht community would develop a systematic investigation 
on yacht resistance at different initial trims, particularly in 
the high-speed range. 

Finally, my congratulations to the author on a well-written 
paper. I am honored to have the opportunity to comment on 
it. 
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Britton Chance, Jr., Member 


Dr. van Oossanen is to be congratulated for his contribu- 
tions to sailing yacht performance analysis and for the success 
of Australia II] with Joop Sloof and the late Ben Lexcen. 
Particularly noteworthy are the methodology of VPP, CFD, 
and experiment; this team put sailing yacht R&D on its present 
path. 

The author's practical guidance with respect to the tank 
testing of sailing yachts is most valuable. Given modern VPP 
requirements and capabilities, it seems clear that the matrix 
method of testing is preferable to the free-running method. 
Testing of appendages is problematical; it would appear that 
a combination of tank and wind tunnel testing is appropriate. 

The author provides boat and sail models useful in VPP 
analyses of performance with speed, heel, and leeway as inde- 
pendent variables. New elements are the inclusion of nonlin- 
ear effects on lift and drag important in maneuvering and 
other highly loaded conditions. I would recommend, as does 
Dr. Milgram in the companion paper, the substitution of side 
force for leeway as an independent variable: leeway is difficult 
to measure and Vmg is far less sensitive to leeway than to 
sideforce. 

The author seems to suggest that reduction in leeway will 
increase Vmg. This is not necessarily the case: the principal 
effect of reduction in leeway is to increase the gap between 
main and jib. Such may or may not be desirable aerodynami- 
cally. Also, when draft and general configuration are fixed, 
leeway reduction generally entails a wetted area/drag increase. 

Further, model and full-size experiments have suggested 
that some leeway is beneficial due to favorable hull-keel inter- 
action. This effect tends to be offset by increased hull form 
drag and increased sail-induced drag as leeway and heel in- 
crease and must be evaluated on a case-by-case basis. 

The author's remarks with respect to the design and effect 
on performance of Challenge Australia’s variable incidence 
keel and comment on the role of wind tunnel testing of ap- 
pendages would be appreciated. 


Olin J. Stephens, Member 


The author has accomplished a formidable task in compiling 
a virtually complete summary of the factors that influence the 
speed of sailing yachts, including the methods by which the 
factors can be evaluated. The usefulness of this paper is greatly 
enhanced by the inclusion of some “real life” data. His criti- 
cism of the lack of such information in the literature is well 
taken and one hopes that this example will be followed by 
others in this field. 

This discusser finds the scope of the paper much too broad 
for comment beyond that suggested by a few points of special 
interest to him. 

Such a point is the reference {item 17 of the second method 
of VPP modeling) to the predetermined maximum heel angle 
at which the sails are flattened or reefed. It is not clear to me 
how this angle is found, and it seems likely that a sophisticated 
optimizing process is needed to find the best point to shorten 
sail, as well as to decide whether it will be best to flatten, 
reducing the sail lift coefficient, or reef, reducing area and 
height. The balance between hull and appendage drag due 
to side force and heel as against sail drive as influenced by 
heel, area, and the lift and drag coefficients must be critical 
determinants of course and speed. 

Considering the comprehensive scope of this paper. there 
may be only one factor in performance that has not been 
treated. The time taken to tack and the distance lost in a tack 
is such a factor. The indication of one possible treatment of 
this subject can be found in a paper by Masuvama et al, 
entitled “Dynamic Performance of a Sailing Cruiser” which 
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was presented at the 11th Chesapeake Sailing Yacht Sympo- 
sium. This has been cited by the author in a different connec- 
tion. A simplified computer program, recently written by Da- 
vid Tabors, a Dartmouth student, yields similar results. 
Neither complete nor perfect, it confirms the feasibility of 
including, in a VPP, a module that might be useful in compar- 
ing racing performance. 

As a committed user of the towing tank in racing yacht 
design, at the Davidson Laboratory, Hydronautics, and a 
Lockheed facility in San Diego, 1 cannot altogether accept 
the author’s suggestion that between 1936 and 1983 the tank 
may have failed to pinpoint winning designs. It should be 
noted that the tank was widely used by competing designers 
all of whom could not win. Winning boats in the 6 and 12 
meter and 5.5 meter classes were routinely tested in model 
form as were ocean racers. That some were better than others 
should be no surprise, but there can be no doubt that there 
were steady improvements in performance and in the under- 
standing of the parameters that influence performance. The 
winning yachts had almost all gone through the tank. Short- 
comings on the part of both designers and tank personnel in 
the proper understanding of scale effects led to failures in the 
early seventies which led to renewed study and eventually to 
the better understanding we have today. The towing tank, in 
the sixties and seventies, was the foundation of the Pratt 
project at MIT. Some uniformed individuals were doubters. 
Those better informed attempted, with success, for a limited 
time, to counter these doubts. 

It is a fact that Dr. van Oossanen and his colleagues shook 
things up in 1983 because they, with the designer, were smart 
enough to look at fundamentals that those like myself had 
neglected, to our sorrow. It seems that their special consider- 
ation of drag due to side force and their use of VPP as an 
aid in creating a balanced design, well suited to expected 
conditions, Jed to a renewed appreciation of the many aspects 
of sailing craft hydrodynamics which aircraft studies had illu- 
minated. These have led to intense study of the sort exempli- 
fied in this excellent paper. 


David S. Greeley, Member 


I would like to congratulate the author for bringing together 
in one place such a wealth of information concerning the 
prediction of sailing yacht speed. The clear presentation of 
the material is welcome to those who work in the field and 
an invaluable reference for newcomers. 

My major comments concern the author’s approach to the 
calculation of lift and induced drag. The “build-up” approach 
to lift generation is satisfactory, and is common for computing 
the forces generated by aircraft, ships, and maneuvering un- 
derwater vehicles. The author has mentioned the most impor- 
tant interaction terms, such as the keel downwash seen by 
the rudder, the interaction between the canoe body and the 
keel, and the interaction between a bulb and the keel. 

The computation of induced drag is far more delicate, how- 
ever, and cannot be approached in the same “build-up” fash- 
ion. Classical aerodynamic texts (e.g., Glauert 1948) demon- 
strate that the induced drag of a biplane consists of the induced 
drag of each wing, plus a substantial interaction term. Such 
interaction terms are not mentioned in Section 4.4 of the 
paper. 

As mentioned by the author in Section 5.3, the CFD com- 
munity has found that a “wake integration” or “far-field” analy- 
sis of the computed flowfield is preferred for the accurate 
evaluation of induced drag for a yacht. This approach considers 
the kinetic energy left in the flow behind the vacht due to 
the generation of lift, rather than integrating the computed 
pressures over the underbody. The analysis is done using the 
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Trefftz plane, as reviewed by this discusser several years ago 
(Greeley 1989). Figure 20 herewith shows the vortex system 
shed behind a sailing vacht with a winged keel, and the 
intersection of the wake behind the yacht with the Trefftz 
plane. Of particular note is that in the Trefftz plane, the 
wake from the keel and the wake from the canoe body are 
indistinguishable, so that they must be considered as a unit. 
Thus, the author's attempt to separate the keel-induced drag 
and the canoe body-induced drag can be quite misleading. 
The same argument applies to the interaction between the 
keel and a bulb. 

Although the author's presentation of formulas for induced 
drag in terms of lift coefficients and aspect ratios is correct, 
1 find them potentially misleading for the unwary. Aspect ratio 
enters into induced drag only because one uses lift coefficients 
based on lateral areas. The actual induced drag of a wing or 
yacht depends on the lift (side force) generated (L), the semi- 
span of the wing (or draft of the yacht) (s), the forward speed 
(V), and the fluid density (p) as follows: 

Induced Drag = (1 + o) L?A27pV’s") 
where (1 + 0) accounts for the departure of the spanwise 
loading from an elliptical shape, as used by the author. As 
shown, the induced drag does not depend on the chord of 
the keel or the aspect ratio, although these latter quantities 
enter if lift or drag coefficients are used. For this reason, most 
people involved in U.S. America’s Cup campaigns prefer to 
deal with the “effective draft” of a yacht underbody rather 
than the effective aspect ratio. This is also consistent with the 
approach used in the IMS formulas, which account for the 
contraction of the keel wake behind the vacht due to the flow 
around a full canoe body. ] would suggest that the author 
might find dealing with (needed) refinements to the IMS 
technique a better approach to the estimation of induced drag. 

Finally, with regard to the author's discussion of induced 
drag not being proportional to the side force squared, Yates 
(1991) has shown that the quadratic relationship between 
induced drag and lift is quite fundamental, even if viscous 
flows. It is therefore not surprising that a regression analvsis 
of data (Appendix 3 of the paper) fails to extract a cubic term 
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in angle of attack for the induced drag, beyond the drag 
associated with hull attitude. 
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Frank DeBord Jr., Member 


I'd like to thank the author for making the effort to publish 
a paper which is a complete overview of this subject and 
provides sufficient detail for readers to gain insight into the 
finer points of predicting sailing yacht performance. Although 
individuals working in this field will probably never be able 
to agree on the detailed formulations given for calculating 
forces acting on a yacht, the material given will be of immedi- 
ate practical value for practicing yacht designers and should 
serve as a source of material that can be used by researchers 
to check new work. 

First, 1 have several general comments and then I will 
address some specific comments and questions on the material 
presented. The author states that prediction of @e speed of 
a sailing vacht is one of the most difficult problems in naval 
architecture. This point should be apparent when one reviews 
the number of effects discussed in the paper and the fact that 
most of the formulations given are still derived empirically. 
It is also demonstrated by the relatively large differences in 
performance exhibited by the vachts participating in the 1992 
America’s Cup challenger and defender elimination series, 
even though these were governed by a fairly restrictive rule 
and many of the competing syndicates expended a great deal 
of effort on design. 

The author states that tank testing, VPP analysis, and CFD 
must be used in conjunction with each other. I agree that 
this is the case and that most successful high-level design 
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programs use this approach. However, it is evident that the 
CFD codes are getting better and in some programs they 
have resulted in Jess reliance on tank testing. I would like to 
ask the author, as a yacht designer with a significant amount 
of towing tank experience, to comment on how he sees the 
relationship of these tools changing in the future. 

In terms of detailed comments and questions, the point 
that some widely used assumptions, such as side force being 
a linear function of leeway, are incorrect in certain cases is 
well taken. For this specific example, I wonder if the author 
has any comments on the variation of side force and dynamic 
roll moment with Froude number, possibly due to changes 
in trim. 

In Figs. 5 and 6 of the paper similar nonlinear drag versus 
angle of attack curves are shown. In one case the author labels 
the nonlinear effect as drag due to separation and in the other 
case it is labeled as drag due to keel wavemaking. | realize 
that the second case is for a rotating kee] configuration and 
that the canoe body has no angle of attack. However, anyone 
who has witnessed model tests of relatively light yachts with 
the canoe body and keel at large Jeeways has probably noticed 
a pronounced keel wave, even though the keel is fixed. Does 
the author have any recommendations on separating these 
two effects? 

Using current practice, as described by the author, perform- 
ance predictions are highly dependent on the method used 
to predict sail forces. These cannot only affect absolute per- 
formance but can also impact the relative performance of 
different yachts. Can the author comment on his level of 
confidence in the IMS sail coefficients as presented in the 
paper and can he make recommendations on a preferred 
approach to improving these? 

The author describes the two basic model testing tech- 
niques in use today. These include the “fixed” and “free” 
model techniques. For the fixed technique he mentions that 
heel angle is typically controlled by moving ballast weights. 
1 do not know the current practice at the Davidson Laboratory 
but historically the models were fixed in heel by the test 
apparatus. This was also the case at Arctec Offshore between 
1984 and 1992 and it is currently the method used at Canada’s 
National Research Council's Institute for Marine Dynamics. 
With respect to the “free” model technique, I would like the 
author’s comments on its applicability when formulations for 
the various components of drag are being developed. Also, 
how are results from this technique used in the author's VPP 
when the first step is to calculate hydrodynamic forces versus 
speed, heel, and leeway? 

I have several comments related to viscous drag calculation. 
First, the 1957 ITTC friction formulation is a correlation line 
rather than a flat plate friction line. As such it includes form 
drag, and it can produce negative “wavemaking” resistance 
at low speeds as well as form factors which are unrealistically 
low. Although this has little effect on predicted full-scale drag, 
it can be significant when comparing results with CFD results. 
A more correct method would be to use the Hughes formula- 
tion which is a flat plate friction formulation in the same form 
as the ITTC formulation. Also, when comparing tank results 
to CFD and trying to separate resistance components, why 
not use variable wetted area and wetted length as proposed 
by Teeters (1993)? This approach would solve the problem 
of separating viscous and wavemaking drag due to heel. These 
considerations, along with the author’s discussion of the Delft 
University parametric series results lead to the very important 
point that whenever wavemaking or residuary drag values are 
published, the method used to compute viscous drag must 
be explicitly stated. 
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Finally, with respect to seakeeping effects, can the author 
comment on the effects of a seaway on lift or drag due to lift? 
Also, does he have any insight into determining added drag 
in the tank in head seas versus the more realistic case of bow 
quartering seas? 

Again, I'd like to thank the author for providing us with 
this excellent overview of the performance prediction problem 
and the forces acting on a sailing yacht. I join him in his call 
for more frequent publication of concise, quantitative results 
in these areas. 


Gunnar C. F. Asker, Member 


Dr. van Oossanen’s comprehensive paper is an excellent 
treatise on the prediction of the speed of sailing yachts. He 
mentions three different methods available for the task. A 
fourth auxiliary empirical and simple method is stationary 
thrust measurement at dockside using now-available digital 
thrust meters capable of measuring thrusts from less than 1 
Tb to over 2000 Ibs. (Spring scales are available for measuring 
much higher thrusts.) This method is simple and especially 
useful for component evaluation in relation to each other. 

The relative performance between a soft sai] and new high- 
tech designs can be evaluated in model testing as well as in 
full-scale measurement. It is as useful in predicting perform- 
ance for true sailing yachts, as for motor sailers. For true 
windships the performance of conventional designs versus 
new developments can be determined. 

For motor sailers the thrust from the engine for various 
rpm’s is known and the percentage of the load from the 
use of sails can be determined for various wind forces and 
directions. Changes in angles of attack against the apparent 
wind, which in stationary dock testing equals true wind direc- 
tion, are readily obtainable. In wind-assist propulsion the aim 
is to develop a sail with high lift per unit of sail area. The 
influence due to changes in design and sail setting can be 
determined with thrust meters. 

I am working with aspirated asymmetric airfoil designs. 
The thrust influence from using the aspiration fan can be read 
with thrust meters. 

Because results are relatively easy to obtain, this simplistic 
method could be considered, especially in model scale, before 
finalizing designs with more sophisticated methods so well 
described in Dr. van Oossanen’s paper. 


Karl L. Kirkman, Member 


As a participant in design/technology for every America’s 
Cup since 1970, I have observed three eras: 

The era of Gift Technology when efforts pretty much had 
to be free to be accepted by a syndicate, and in which such 
work was mistrusted, leading to the loss of the Cup which 
led to 

Shopping Cart Technology in which syndicate manage- 
ments traveled the aisles of the technology stores and loaded 
up on items which caught their fancy, or were advocated 
or pushed by individual partisans, which raised costs and 
swamped management resources and manpower pools, 
forcing 

Educated Consumer Technology in which the potential im- 
pact of technologies is assessed by means of forecasting costs 
and benefits of alternatives, and investing accordingly. 

This evolution requires the ability to estimate the effects 
of technology on the race-winning probabilities which | will 
return to in more detail later. 

In preparation for the current America’s Cup technology 
efforts, I solicited the opinions of a number of knowledgable 
observers to hindcast their “predicted” pre-race probability 
of success of individual syndicates, and the budgets each 
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expended. A sort of consensus version of these results is shown 
in Fig. 21 with this discussion. ] have anchored the straight- 
line fit of the nominal effectiveness by the price of a minimal 
effort of 10 million, and almost sure success with a price of 
several hundred million. This presentation shows two points 
clearly: 

1. More money is better, but 

2. Effectiveness of utilization of resources is even more 
important for typical syndicates. 

A syndicate could gain as much effectiveness by careful 
focusing of expenditures as if they had an additional budget 
of perhaps 50-100 million dollars! 

Indeed, time and the availability of a pool of skilled partici- 
pants are terribly hard limits on any syndicate, and again 
careful management is important. 

Against this background, I have two points to make: 

1. The author notes that no definitive conclusions have 
been reached concerning the minimum requirements for 
model size and turbulence stimulation. 1 think that it more 
accurate to say that, while no hard published criteria enjoy 
widespread endorsement by an organization such as ITTC, 
we in the yacht-testing community have learned a great deal 
about the subject over the past 20 years, and that rational 
criteria can be developed once there is a consensus as to the 
purpose of the tests. 

To explain this apparent contradiction let me say that we 
have progressed with yacht prediction methods to the point 
that experiments cannot effectively continue to be used in 
old ways, that is, to validate particular design solutions. 

Accordingly, earlier criteria based upon experimental re- 
sults ranking the alternatives correctly have less weight be- 
cause we need to match absolute levels if we are to gain an 
understanding of the critical deficiencies of CFD models. 

Further, we must use experiments to support the validation 
of prediction tools. In this context, many of the old rules of 
thumb regarding experimental error are simply out the win- 
dow because we are confronted with producing experimental 
results for which the uncertainty levels must be small com- 
pared to the errors in prediction we are checking. I propose, 
for example, the working rule that the experimental uncer- 
tainty must be less than one fourth of the uncertainty level 
of validation sought. Under these circumstances, the experi- 
mental error must be small indeed. 

We have been developing such experiments in a number 
of areas including yacht testing, and given them the name 
Very Low Uncertainty (VLU) experiments to denote nontradi- 
tionally low levels of experimental uncertainty. 

When undertaking such testing, the need arises to invoke 
the following innovative concepts: 
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® Scientific, as opposed to industrial, standards for experi- 
mental planning; 

¢ Blocks of experiments involving common geometries in 
different types of test facilities; example tests of model-scale 
keels on Froude-scaled conventional models, Froude-scaled 
testing of full-scale keels, and wind tunnel testing of sub-scale 
models at appropriate Reynold’s numbers. 

¢ Nontraditional concern for determining components of 
resistance, and for improving extrapolation procedures. 

® Design of experiments to involve appendage models (that 
is, appendages designed to perform at model-scale conditions), 
not mode] appendages. 

® Interactive use of CFD with experiments to estimate bias 
effects and extrapolation techniques. 

¢ Development of new disciplines for validating prediction 
tools using systematic blind approaches, and aimed at charac- 
terizing the uncertainty level of a prediction tool result in 
terms of geometric and performance domain parameters. 

As an example of proposed standards for model size in VLU 
tests, let me suggest Fig. 22, which was developed for IACC 
class yacht testing; this is not entirely new and represen® 
application of suggestions made in my paper to the 1979 
Chesapeake Sailing Yacht Symposium on this subject. How- 
ever, in this case the selection of the cutoffs represents an 
interpretation of the requirements for VLU tests which might 
be more demanding than may be acceptable for other pur- 
poses. 

Further, let me point out the need for entirely new types 
of testing. As our ability to predict performance moves to 
unsteady phenomenon—seakeeping, tacking performance, 
unsteady aerodynamics of rigs, etc.—the need arises to de- 
velop experiments to support these tools with VLU experi- 
mental results. It was reported in the wake of the last America’s 
Cup competition that a number of syndicates were unable to 
obtain results of tests in rough water adequate to validate 
added resistance predictions. In anticipation of just this diffi- 
culty we expended significant resources at the University of 
Michigan to build a VLU capability which-supported signifi- 
cant advances in the computation of performance in waves. 


However, this was just a start and further developments are 


critical to continued progress so that we can model oblique 
seas, unsteady rig performance, etc. 

2. Turning to the “Final Remarks” section of the paper, I 
have had the opportunity in recent years to study how we 
can encourage innovation in marine design from which I 
would extend the author’s vision somewhat from that pre- 
sented: that the past was model testing, the present is VPP 
modeling, and the future is CFD. 

First, in designing competitive yachts we must think beyond 
the VPP level to the concept of probability of winning a race 
or series, and must engage in adaptive design—matching 
particular opponent’s characteristics. Thus, measures of sys- 
tem performance such as win probability as computed by a 
Race Modeling Program (RMP) are the future, not CFD. 

Having made this generalization, it is clear that speed pre- 
diction is a subset of the complete design process which must 
deal at early stages of design with: 

(a) optimization of particular designs, and 

(b) provision of robustness—acceptable off-design per- 
formance. 

Once one makes a commitment to this approach of simulat- 
ing the complete race course performance against opponents, 
whole new opportunities appear: using this simulation to eval- 
uate candidate technologies in the context of $/% win probabil- 
ity improvement, providing networks for design team mem- 
bers to interact with technology contributors at distributed 
locations (and to provide a tool to manage schedules and 
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deliverables), and particularly the ability to place intellectual 
effort into inventory are a few of the benefits. 
Thank you, Peter, for a fine addition to our literature. 


Author’s Closure 


The author would like to thank the discussers for their 
comments. Many interesting and important topics have been 
raised which have not been fully amplified in the body of the 
paper. The author particularly wishes to thank Olin Stephens 
for his contribution and kind remarks. 

In replying to the various questions and points raised by 
the discussers, the author has identified the subject matter 
concerned by way of providing a subtitle in addition to identi- 
fying the name of the discusser. 


DUT POLYNOMIALS FOR RESIDUAL RESISTANCE 


Dr. Fung raises some very important points in connection 
with the accuracy and applicability of the Delft polynomials 
for the prediction of the residual resistance of sailing yachts. 
As stated in the body of the paper, the Delft University of 
Technology (DUT) is currently reexamining the accuracy and 
suitability of these polynomials, repeating the tests with some 
of the earlier models when considered necessary. Some results 
of this work have already been presented to the IMS Interna- 
tional Technical Committee and the author is looking forward 
to the general release of these data. The author will now 
address in turn each of the comments made by Dr. Fung: 

l. The coefficient of a; in Table 2 for the Froude number 
value of 0.175 is indeed —31.914. 

2. It should be noted that the terms in the DUT polynomials 
were not obtained by multiple regression. The coefficients in 
both the low- and high-speed polynomials were based on a 
straightforward least-squares method after a study had been 
made of which terms in the independent parameters were 
most significant. If a multiple regression method had been 
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used, the polynomials for the different Froude numbers would 
not have yielded the same set of terms for each speed (i.e., 
not the same nine terms in all of the low-speed polynomials 
and not the same five terms in all of the high-speed polynomi- 
als). It is probably this factor, more than any other, which has 
led to doubts concerning the accuracy and validity of the 
polynomials, particularly since the value of the coefficients of 
the individual terms do not possess a smooth transition with 
increasing Froude number in all cases. 

3. The author agrees with Dr. Fung that when requiring 
accurate predictions the polynomials should not be used when 
one or more of the individual hull form parameters is greater 
or less than the average value of that parameter, plus or minus 
the standard deviation thereof, respectively. On that basis, 
the applicability requirements of the low-speed polynomials 
are as follows: 

—Ly:/By, between 3.042 and 3.986 

—Ly/Bya, between 2.609 and 3.419 

—By,/Tc between 2.41 and 10.504 

—cC, between 0.537 and 0.579 

—Lw/V cs between 4.435 and 7.053 

—LCB between — 1.254% and — 4.436% of the waterline length 

aft of the Ly,/2 location 

The applicability requirements of the high-speed polynomi- 
als are then as follows: 

—Lwy/Bw, between 3.577 and 4.423 

—Lw_/B,, between 3.019 and 3.775 

—Bw,/Tc between 7.04 and 14.188 

—C> between 0.538 and 0.560 

—Lw1/V cs" between 6.779 and 7.941 

—LCB between — 2.365% and —5.141% of the waterline length 

aft of the Lw,/2 location 

4. One disadvantage of using discrete polynomials for each 
Froude number is that an interpolation procedure needs to 
be adopted to obtain a continuous curve of resistance versus 
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speed. Another disadvantage is that, sometimes, the resistance 
values calculated for successive Froude numbers are too close 
together or too far apart, causing interpolation difficulties. 
The reasons for this to happen are correctly explained by Dr. 
Fung, namely, inaccuracies in the predictor equations due to 
an inhomogeneous distribution of the independent geometry 
parameters of the hull form. In the sample given by Dr. Fung, 
a low value of the prismatic coefficient (0.53) is combined 
with a high (negative) value of the longitudinal position of 
the center of buoyancy (—5.0% of waterline length), for the 
model with the smallest length-beam ratio and the lowest 
length-displacement ratio. In this region the prediction capa- 
bility of the predictor polynomials obviously breaks down. If 
the basic independent parameters of the DUT model series 
are recognized as being length-beam ratio (Lw,/Bw,), beam- 
draft ratio (Bw./T,), prismatic coefficient (Cp) and the longitudi- 
nal position of the center of buoyancy (LCB), it would require 
3* = 81 models to cover the chosen variation in each of these 
parameters, on varying each parameter systematically three 
times. In utilizing only 39 models, it will be clear that many 
of the so-called “corner” models, exhibiting combinations of 
the highest or smallest values of the independent geometry 
parameters, were not tested. 

It should also be noted that whereas the example chosen 
by Dr. Fung might have been expected to fall within the 
applicability region of the low-speed polynomials, this cer- 
tainly would not have been expected to fall within the applica- 
bility region of the high-speed polynomials. The chosen 
length-beam and length-displacement ratio values are lower 
than the minimum values of the 15 models considered in the 
development of the high-speed polynomials, for the speed 
region in excess of Froude number values of 0.45. 

5. The author agrees with Dr. Fung that the value of the 
longitudinal position of the center of buoyancy is often seen 
to influence the continuity of the resistance curves obtained 
from these polynomials significantly more than any of the 
other parameters. It should: be noted in this respect that the 
39 models display less systematic variation in this parameter 
than in any of the other independent parameters. Also, a 
change in LCB requires a change in the lines of the canoe 
body, whereas a change in length-beam or beam-draft ratio 
only can be achieved by maintaining the section or waterline 
shape. It follows that particular care should be taken in 
applying these polynomials for LCB values outside of the 
proposed applicability ranges given above and in applying 
these polynomials for yachts with section shapes which are 
significantly different to those of the three parent models 
used. 


VALUE OF WIND TUNNEL TESTING OF SAILING YACHT AP- 
PENDAGES 


In his contribution, Mr. Chance suggests that the testing 
of appendages in a wind tunnel is necessary as well. The 
results of towing tank tests reveal, however, that both the 
keel and nidder of a sailing yacht possess important lift and 
drag components dependent on the proximity of the lifting 
surface to the water surface. These lift and drag components 
cannot be modeled in a wind tunnel. For that reason, very 
little attention is given to wind tunnel testing in the body of 
the paper. Only the viscous resistance can be established in 
a wind tunnel, which in some cases might be important enough 
to want to carry out wind tunnel testing, for example, in the 
case of evaluating the viscous resistance of the large bulbs 
fitted to the keels of International America’s Cup class yachts. 
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USE OF SIDE FORCE IN LIEU OF ANGLE-OF-ATTACK 


The author does not agree to the usefulness of the substitu- 
tion of side force for leeway angle as suggested by Mr. Chance, 
which is nowadays often carried out by workers in the field. 
While he agrees that the measurement of leeway angle is 
difficult, particularly on the full-scale yacht, the use of leeway 
angle in VPP modeling will always lead to greater prediction 
accuracy. As explained in the body of the paper, the side force 
at some angle of heel is not zero when the leeway angle is 
zero. In addition, the side force is often not a linear function 
of leeway angle at all. Accordingly, when solving the velocity 
triangle to obtain the true and apparent wind speeds and 
angles, a more accurate result is obtained when the relation 
between side force and leeway is modeled as best as possible 
and included in the VPP, with the leeway angle treated as 
the independent variable. 


EFFECT OF LEEWAY ON OVERALL PERFORMANCE AND USE 
OF ROTATABLE APPENDAGES 


Mr. Chance raises the question of the effect of leeway on 
overall upwind and downwind performance. This topic has 
become the subject of considerable discussion ever since it 
has become popular to develop side force through the use of 
rotatable appendages, thereby restricting the leeway that the 
canoe body adopts. It has recently been the subject of many 
studies in various America’s Cup projects. 

In 1987, the 12-Meter Class yacht USA sported a forward 
rudder and only a relatively small strut to which the ballast 
bulb was attached. The development of side force was mainly 
brought about by the use of the forward and aft rudder. In 
the 1992 America’s Cup match about half of the contesting 
yachts utilized rotatable appendages for the production of side 
force. Figure 23 shows the different appendage configurations 
used on the 1992 crop of America’s Cup yachts. As in 1987, 
however, the best yachts proved to be those with a fixed keel 
fin (with trim tab). 

The reason for yachts that are fitted with more than one 
rotatable appendage (the rudder) generally not performing as 
well as those yachts utilizing leeway angle to control side 
force production is not yet entirely clear. Model tests have 
repeatedly revealed that the resistance associated with the 
production of side force, when the canoe body is at zero 
leeway, is appreciably less. This is entirely logical since the 
canoe body is known to produce a relatively sharp increase 
in resistance, even at moderate angles of leeway. This is indi- 
cated in Fig. 24 (van Oossanen 1992) in which the increase 
in resistance, as measured on a 1-to-7 scale model of an JACC 
yacht, is shown for the case where the rotatable keel fin is 
fixed at zero angle-of-attack relative to the centerline and then 
rotated. In this case, when the keel fin is rotated, the canoe 
body and the bulb are fixed at zero leeway. 

It is clear from Fig. 24 that the resistance increase associated 
with the production of side force is significantly less when 
the canoe body and large kee] bulb are held at zero leeway 
angle, particularly when the rudder is loaded up to the same 
side force as is the case when the complete hull is at the 
nonzero angle of leeway, i.e., when the rotatable keel is fixed 
on the centerline. 

The fact that this decrease in resistance is not seen to 
improve the overall performance on the race course can be 
attributed perhaps to the fact that it is much more difficult 
for the crew of a yacht fitted with more than one rotatable 
appendage to fine-tune the setting of these appendages for 
optimum performance. It should be realized that it takes 
between about 30 and 60 seconds for an IACC yacht, after a 
tack for example, to achieve the maximum speed in a given 
set of wind conditions. Any disturbance in the form of a helm 
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Fig. 23. Overview of types of appendages utilized on yachts 
participating in 1992 America’s Cup, particularly indicating the 
various types of rotatable appendages 


or sail shape adjustment will again lead to the yacht having 
to achieve a new equilibrium condition, requiring more time. 
The greater the inertia of the yacht the longer this will take. 
By having two appendages to continuously adjust instead of 
only the rudder, a much longer time is required for the crew 
to obtain the new optimum setting of the appendages for a 
new set of equilibrium conditions if this is done by “feel” or 
intuition only. It is necessary, therefore, in adopting more 
than one rotatable appendage, to also design and install some 
information system for the crew to use in setting these append- 
ages. In 1992, the author was successful in developing, build- 
ing, and fitting an accurate leeway-measurement device on 
Challenge Australia for the crew to use in setting the rotatable 
keel fin. Lack of time (basically due to lack of funding) pre- 
vented the crew from becoming accustomed to this method 
of sailing the yacht. Other methods were utilized on other 
yachts. On the whole, however, only the New Zealand yacht 
seemed to have been able to overcome this problem of dy- 
namic contro] during the 1992 America’s Cup. 


DETERMINATION OF OPTIMUM HEEL ANGLE FOR SAIL FLAT- 
TENING AND SAIL SHORTENING 


In answer to Mr. Stephen’s question concerning how, in 
the author's VPP, the optimum heel angle is determined at 
which the sails are flattened or sail area is shortened or both, 
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Fig. 24 Resistance due to heel and side force for a heel angle 


of 20 deg and a Froude number of 0.35, as measured on a 1-to- 

7 model of an IACC yacht. This result indicates that less induced 

resistance is incurred in the case when side force is generated 

through keel rotation, with canoe body at zero leeway, rather 

than through leeway angle with keel locked on centerline. All 
values are for full-scale yacht 


the author wants to point out the following. In every project 
the impact on performance of at least three values of an input 
heel angle at which the sails are flattened and then shortened 
is investigated. The optimum value is determined by investi- 
gating the results obtained. This is a relatively simple matter. 
In the program the sails are first flattened after which one of 
two methods of sail shortening can be adopted. The first of 
these shortens the foot of all the sails, leaving the hoist values 
the same. The second method covers the traditional slab reef- 
ing situation. 


METHOD FOR ASCERTAINING TACKING AND GYBING PER- 
FORMANCE 

The time and distance lost during tacking and gybing has 
indeed not been addressed in the paper. The author agrees 
that a mathematical model such as given by Masuyama et al 
(1993) can indeed be included in a VPP. For this to be really 
accurate, however, it is necessary to carry out further investi- 
gations to determine how some of the important coefficients 
in these maneuverability models depend on the geometry of 
canoe body, keel, and rudder. 


DISBELIEF IN TOWING TANK RESULTS IN U.S. BETWEEN 
1936 AND 1983 


The author’s statement that between 1936 and 1983, in 
the United States, the towing tank seemed to have failed to 
pinpoint the winning or best design in a number of instances, 
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was certainly not in reference to any of the well-known design 
projects carried out by Mr. Stephens. In fact Olin Stephens’s 
belief in model testing and his J-Class and 12-Meter Class 
designs such as Rainbow, Intrepid, Courageous, and Freedom, 
to name just a few, which were usually the result of diligently 
performed model testing programs, will forever constitute 
highlights in the history of yacht design. The author's remark 
specifically referred to the Mariner and some subsequent 
America’s Cup projects which seemed to have triggered a 
significant disbelief in the U.S. in the capability of the towing 
tank in general. In fact, when the author was working on the 
design and the testing of Australia 1] together with the late 
Ben Lexcen at NSMB during 198] and 1982, which fact did 
not go unnoticed by a number of U.S. clients visiting NSMB 
in Wageningen during that time, one of these, Mr. Geo B. 
Drake, naval architect and marine engineer from Cutchogue, 
N.Y., subsequently sent a letter to NSMB with a newspaper 
clipping which contained the following: 


... The Independence testing program is a departure 
from tradition. Rather than spend money on testing in 
tanks with models, Hood intends to test the real thing 
under real conditions. “Since 1967 (and Olin Stephens’s 
Intrepid), we seem to have gone backward in 12-Meter 
design,” Hood said. “Mariner, the Britton Chance design 
that Turner skippered in the 1974 elimination series, has 
become the classic case against testing tanks. In the tank, 
the Mariner model beat everything, something that was 
hardly duplicated at sea. I don’t think we know enough 
about tank testing,” Hood said. “And those who have 
done it have gotten into trouble. ...” 


Mr. Drake’s letter went on to state: 


“|... We all agree with Ted Hood—model testing is a 
waste of money. ...” 


This illustrates the sentiment in the U.S. in that period 
and, according to the author, is one of the main reasons for 
the U.S. losing the America’s Cup to Australia in 1983. This 
infamous letter was given to the author by the director of 
NSMB, Prof. J.D. van Manen, as an encouragement to turn 
the Australia II project into success. 


CALCULATION OF INDUCED DRAG 


Dr. Greeley pinpoints perhaps one of the weaknesses in 
the author’s VPP, i.e., the calculation of induced drag. The 
author adopts a modular calculation thereof, computing the 
induced drag associated with the lift on canoe body, keel, and 
rudder separately, neglecting any interaction between the 
individual components. Dr. Greeley states that since the far- 
field integration of kinetic energy yields a seemingly accurate 
assessment of total induced drag, it is more appropriate to 
consider the canoe body, keel and rudder as a single unit, 
rather than as separate bodies. To this the author agrees and, 
indeed, the results of the analysis of the resistance associated 
with the development of side force presented in Appendix 4 
of the paper, in which the induced resistance is seen to be 
basically a function of the square of the total lift on the model, 
confirms that this is the case. 

While the theoretical basis might be lacking, the author 
has nevertheless found that the summation of the induced 
resistance as calculated for the individual components will 
generally also yield the same level of total induced drag. If 
this is then a valid procedure, the expression of induced drag 
as the summation of the induced drag components of the 
canoe body, keel, rudder, and any other appendages, is of 
some importance since this will allow the study of the effect of 
changes in overall design of the canoe body and the individual 
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(lifting) appendages on induced resistance. Such is not possi- 
ble when expressing the total induced resistance as a single 
expression dependent only on maximum draft, or effective 
draft, utilizing a single (1 + o) term to account for the deviation 
of the total configuration from one with an elliptic spanwise 
loading. As can be seen from the analysis presented in Appen- 
dix 4 of the paper, this term accounts for almost a doubling 
of the induced drag at moderate Froude number values in 
that case (see also Fig. 13). When accounting for the induced 
drag of the individual components, the effect of the utilized 
(1 + o) factors is almost negligible. 


EXISTENCE OF A DRAG COMPONENT DUE TO LIFT, NOT 
PROPORTIONAL TO SQUARE OF LIFT; FURTHER DISCUSSION 
ON DEFINITION OF COMPONENTS OF RESISTANCE ASSOCI- 
ATED WITH LEEWAY AND HEEL 


Notwithstanding the recent work of Yates (1991) referred 
to by Dr. Greeley, the existence of a drag component due to 
lift, but not dependent on the square of the lift, is really 
beyond any doubt. In addition to the references given in the 
paper on this topic, mention should also be made of the 
summary on this subject given by Kitichemann (1978, pp. 160- 
168). This drag component is due to flow separation effects 
and is, according to many workers in the field, proportional 
to the cube of the angle-of-incidence of the flow. It is agreed 
that this component is often not discernible in the results of 
model tests and it is, indeed, not noticeable in the results 
presented in Appendix 4 of the paper. All so-called nonlifting 
bodies and all lifting bodies of small aspect ratio possess this 
component however. 

The author was recently involved in a project for which 
tests were also carried out with a model of a racing yacht 
with no keel or rudder fitted. Tests were subsequently carried 
out with this model with only the keel fitted to the canoe 
body, then with only the rudder fitted, and finally, with keel 
and rudder fitted. The forces and moments on canoe body, 
keel, and rudder were measured separately. In this way many 
of the interaction effects on side force and resistance, between 
canoe body and keel, between canoe body and rudder, and 
between keel and rudder, were determined. 

The side force on the canoe body was found to possess a 
dominant nonlinear component dependent on the square of 
the angle-of-attack Bcos@ as shown in Fig. 25. In accordance 
therewith, and in line with the discussion in paragraph 2.3.6.2 
of the paper, the resistance due to heel and side force revealed 
the canoe body to possess a component not dependent on the 
square of the lift coefficient. This is shown in Fig. 26, in which 
the total resistance associated with heel and side force is 
graphed against the square of the (total) lift coefficient. The 
subsequent analysis of these and similar measurements re- 
vealed that the coefficient of induced resistance could be 
accurately represented by an expression involving the sum of 
a term in the square of the coefficient of linear lift and a 
term in the product of the coefficient of linear lift with the 
coefficient of nonlinear lift. This last term is, in effect, a 
function of the cube of the angle-of-attack. 

All canoe bodies possess such characteristics. These charac- 
teristics are retained when the keel and rudder are added. In 
the latter case, however, it is sometimes difficult to pinpoint 
this behavior since the nonlinear lift component is relatively 
small compared to the lift of efficient lifting appendages, as 
can be seen from Fig. 2 in the body of the paper. 

Further correspondence with Dr. Greeley has revealed 
that he does not disagree with the existence of a resistance 
component dependent on the cube of the angle-of-attack but 
that he does not want to consider this as part of the induced 
drag. His definition of “induced drag” is restricted to that 
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Fig. 25 Results of side force measurements on a model of a 

Sailing yacht without keel and rudder fitted, at a Froude number 

of 0.35 and heel angles of 0, 15, and 30 deg, which indicate the 
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Fig. 26 Results of resistance measurements for same model 

for which side force is given in Fig. 25. Excess resistance 

associated with leeway and heel is plotted as a function of square 

of total side force. Here, total induced resistance (defined by the 

author as total resistance due to lift) is not proportional to square 
of total side force 
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component of resistance which is dependent on the square 
of the linear lift. He suggests that the resistance component 
dependent on the cube of the angle of incidence should be 
inchided in the “resistance due to leeway” category, which 
also includes the excess resistance associated with an increase 
in viscous and wave resistance because of the attitude of the 
hull. 

In the paper, the author has subdivided the total resistance 
associated with leeway and heel into induced resistance, de- 
fined as the total resistance directly associated with the level 
of the side force or lift, and the “resistance due to hull atti- 
tude,” which has no relation to the level of the side force or lift 
at all. As such, the induced resistance includes the resistance 
dependent on the square of the linear lift and the resistance 
dependent on the nonlinear lift to the power 1.5, i.e., depen- 
dent on the cube of the angle-of-attack, since the nonlinear 
lift is dependent on the square of the angle-of-attack. 

This subdivision of the additional resistance which is in- 
curred when a sailing hull adopts leeway and heel is perhaps 
more fundamental than the subdivision defined by Dr. Gree- 
ley, although the argument could be made that when using - 
potential flow CFD alongside model testing, it is possible to 
calculate the induced resistance associated with the square 
of the (linear) lift, which could be used to divide the additional 
resistance due to leeway and heel into a component dependent 
on Froude scaling (linear lift and induced drag associated 
with the square of the linear lift), well described by potential 
flow models, and the excess side force and resistance which 
then becomes a complicated mixture of components with both 
a Reynolds and a Froude origin. 

It is perhaps not important which of these definitions are 
adopted as long as it is recognized that the increase in resist- 
ance due to leeway and heel has components associated with 
the level of the side force and components that are not associ- 
ated with the side force at all. The following gives a schematic 
presentation of the two methods of subdividing the excess 
resistance associated with leeway and heel. 


Author's subdivision Discusser Greeley’s subdi- 


vision 





Resistance due to heel 
(at. zero lift); 


Induced resistance 
(resistance due to 
linear lift only); 


Resistance due to heel 
(at zero lift); 

Induced resistance 
(resistance due to 

total lift, divided 


into two components; 


Resistance due to leeway 
(all excess resistance, 
including that associated 
with nonlinear lift). 


Resistance due to 
attitude of hull 

(all excess resistance not 
associated with lift). 


Use oF CFD IN THE FUTURE 


Mr. DeBord raises the topic of how CFD might be used 
in the future in relationship to towing tank work, as tools in 
yacht design. In this respect the author wishes to state that 
while he is a firm believer in the towing tank for obtaining 
design guidance in important sailing yacht projects, he is at 
the same time convinced that the time will come when the 
towing tank will be used only to answer fundamental questions 
in relation to the physics of complicated hydrodynamic prob- 
lems or to carry out validation work in connection with the 
development of new numerical models. 

As stated in the body of the paper, CFD tools are at present 
generally not yet accurate enough and cannot be relied on 
solely in important sailing yacht projects. At the current rate 
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of development, however, it should soon be possible to en- 
tirely rely on CFD for the prediction of side force, wave 
resistance, and induced resistance. The reliable prediction of 
viscous resistance, for which accurate flow separation model- 
ing is necessary, will take longer to achieve. 


DEPENDENCE OF SIDE FORCE AND DYNAMIC ROLL MOMENT 
ON FROUDE NUMBER 


The analyses of towing tank tests usually reveal a significant 
effect of forward speed on side force and dynamic (restoring) 
roll moment. Particularly as the heel angle increases, when 
the proximity of the kee] and other lifting surfaces to the free 
surface is increased, the efficiency of the hull in producing 
side force is substantially reduced. In the author's VPP this 
is modeled by defining the effective aspect ratio of keel and 
rudder as a function of the depth of submergence, taking into 
account as best as possible the wave profile along the hull, 
the sinkage and trim of the hull and, in the case of the rudder, 
the bow-down trim due to the moment of the forward force 
on the sails. As a certain maximum depth of submergence is 
reached, the effective aspect ratio very quickly reduces to 
become approximately equal to the geometric aspect ratio. 
Part of this numerical model is given in Paragraph 4.3.6 of 
the paper. 


SEPARATION OF RESISTANCE INCREASE DUE TO FLOW SEPA- 
RATION AND FREE SURFACE EFFECTS 


The increase in resistance due to flow separation and free 
surface effects can be separated and independently accounted 
for by using a multiple regression analysis technique to “cross- 
fair’ the measured data for side force and resistance. This 
technique is described in Paragraphs 3.5.4 through 3.5.6.2 of 
the paper. In this way the effect of the independent variables, 
such as leeway angle (or keel rotation angle), heel angle, 
and Froude number can be accurately determined, yielding 
analytical expressions with which interpolation and extrapola- 
tion procedures can be performed. If the coefficients of the 
terms in the angle-of-attack are found to be dependent on 
Froude number, it may be assumed that the free surface has 
a major effect on the resistance associated with side force, 
which is usually the case when a surface wave appears at the 
location of the lifting appendages. If the independent variables 
are not utilized when cross-plotting or cross-fairing the mea- 
surement data, such an assessment cannot be made, impeding 
the accurate extrapolation of the model data to full scale. 

In finding that the increase in resistance associated with 
the development of side force is dependent on terms in the 
angle-of-attack only, without the coefficient thereof showing 
any major dependence on speed, this resistance increase is 
usually due to flow separation. Without performing further 
tests at various Reynolds numbers, no scaling relations for 
full-scale extrapolation can be deduced. On occasions the 
author has utilized a factor proportional to the coefficient of 
frictional resistance to derive an approximate extrapolation 
rule. Such practice without further supporting data, however, 
constitutes a dangerous procedure. 


RELIABILITY OF THE IMS SAIL COEFFICIENTS 


The author agrees with Mr. DeBord’s statement that the 
final performance prediction obtained with a VPP is signifi- 
cantly dependent on the sail coefficient values used. Even in 
the case of the author's VPP, in which the numerical model 
for side force and resistance is totally dependent on the hydro- 
dynamic properties of the hull, the final calculation of wind 
speeds and angles is still dependent on the lift and drag 
coefficients of the sail plan. A critical review of available sail 
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coefficient values or sail force models, in important projects, 
is therefore essential. 

The lift coefficient value of the various sails in the IMS sail 
model is basically a function of apparent wind angle only. No 
dependence, for example, on aspect ratio of the sails has 
been defined. This is a serious deficiency, particularly when 
comparing the performance of yachts with sail plans of differ- 
ent heights and different foot lengths. The aspect ratio only 
enters into the calculation of the induced drag of the sails. 

Another deficiency is that associated with the fact that all 
yachts are considered to make the same amount of leeway at 
a particular apparent wind angle. For upwind performance 
estimates involving different keel concepts, including those 
that allow the hull to sail at zero leeway angle for which the 
effective sheeting angles at a particular apparent wind angle 
are greater, the IMS sail coefficients yield inappropriate con- 
clusions. The author has spent considerable time in suc! 
projects developing expressions for sail lift and drag which 
allow for a more accurate estimation of some of these effects. 

Other shortcomings are associated with the fact that no 
account is made for differences in leech roach of the mainsail 
or leech hollow of jib or genoa, that is, in the spanwise load 
distribution associated with planform. The sail coefficient val- 
ues are based on a triangular planform as defined by I, J, P 
and E. A certain amount of mainsail-genoa overlap and main- 
sail leech roach can be accommodated by correcting the sail 
area factor. 

Of all the inaccuracies that are inherent in VPP’s today, 
the error involved in the calculation of the aerodynamic thrust 
and heel forces is probably the most significant, with the 
exception perhaps of the error involved in predicting perform- 
ance when VPP’s do not take the essential fourth equilibrium 
condition into account, namely, equilibrium of all forces in 
the horizontal plane. 


APPLICABILITY OF FREE-RUNNING MODEL TEST METHOD 
FOR DEVELOPMENT OF VPP MODELS . 


The basic disadvantage of the free-running test method in 
relation to the testing of sailing yacht models is that the effect 
of heel and leeway on side force and resistance cannot always 
be adequately separated. This is due to the fact that as leeway 
angle is increased, so is heel angle increased. It is not possible 
to test the model at, for example, a matrix of leeway angles 
comprising the values of 0, 3, 6, and 9 deg for heel angles of 
0, 10, 20, and 30 deg. When building a mathematical model 
for use in a VPP, this is a distinct disadvantage. It is possible 
to carry out free-running tests at three different positions of 
the vertical center of gravity of the model, thereby obtaining 
three distinct relationships for the heel angle as a function of 
leeway angle, with which some of the model-building prob- 
lems can be addressed when utilizing a multiple regression 
technique for the analysis of the results. The problem of 
determining the resistance due to heel at zero side force, at 
10, 20, and 30 deg of heel for example, cannot be addressed 
in the free-running test method. 

For this and other reasons the Maritime Research Institute 
Netherlands (MARIN) instigated a research project in 1988 
to simplify its free-running test method referred to in the 
body of the paper, so as to be able to combine the capability 
of significantly varying transverse stability characteristics of 
the model within the existing method by utilizing a servo- 
steered weight transversely across the deck of the model. This 
would allow, within one and the same test, to significantly 
alter the leeway-heel characteristic of any model. When the 
author left MARIN in 1989, however, MARIN had not yet 
succeeded in obtaining the desired results from a trial experi- 
ment. Since 1989 the author has used the captive test method 
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at the Delft University of Technology, up to now testing over 
30 models, from which results most of his VPP “building 
blocks” have been derived. In this laboratory, heel angle and 
transverse stability are varied by moving a weight transversely 
across the deck. 


SUITABILITY OF THE 1957 ITTC FRICTION FORMULATION 


Mr. DeBord correctly points out the error that can arise 
from using the 1957 ITTC friction coefficient formulation in 
certain cases. It is important, particularly when dealing with 
the results of tests with small models, to be aware of this. 

When using small models in combination with a three- 
dimensional extrapolation procedure, with a physically correct 
value of the form factor, special care must be taken when 
using the ITTC 1957 friction line when Reynolds numbers 
are less than about 8 X 10°. In this case the ITTC friction 
formula deviates from other friction formulas in that it inher- 
ently includes a so-called correlation allowance. The 1957 
ITTC formulation was specifically chosen to allow for a better 
correlation between the results of extrapolation of resistance 
data derived from small and large models when using a two- 
dimensional extrapolation procedure. When this friction line 
was adopted, insufficient knowledge was available on the ap- 
propriate value of the form factor to be applied and the two- 
dimensional extrapolation method, which employs an incre- 
ment in Cr (AC;), was still widely used. 

The problem of not obtaining a good correlation between 
the results obtained from small and large models stemmed 
from the problem of not using different levels of viscous 
resistance for different model sizes. Adopting the same AC; 
value for both small models and large models, as was often 
the case in two-dimensional extrapolation procedures (the 
value of +0.0004 was frequently adopted, irrespective of the 
Reynolds number), amounts to a lower level of viscous resist- 
ance for the small model than in the case of a large model. 
At a Reynolds number of 1 x 10°, for example, the addition 
of 0.0004 to the Schoenherr friction coefficient formula [often 
referred to as the American Towing Tank Conference (ATTC) 
friction line], frequently used prior to 1957, leads to a value 
of 0.004409 + 0.0004 = 0.004809, which is equivalent to 
using a form factor of 1.0907. For a larger model at, say, a 
Reynolds number of 5 x 10’, this same procedure leads to 
0.002289 + 0.0004 = 0.002689, equivalent to using a form 
factor of 1.1747. This discrepancy between the extrapolated 
results of tests with small and large models was incorrectly 
solved by adopting a friction line with an artificial increase 
in Cy at low Reynolds numbers, so as to obtain a higher level 
of the viscous resistance for small models. It is now clear that 
a constant form factor value and a physically correct C; value 
must be utilized to arrive at the properly scaled level of viscous 
resistance. This assumes that significant flow separation on 
the model is not present. 

The differences between the ITTC and ATTC friction lines 
are as given in Table 1] herewith. 

The physically correct level of flat plate frictional resistance 
at lower values of the Reynolds number is approximately given 
by the Hughes friction formulation (Hughes 1954). This has 
now been confirmed by various studies. After carrying out a 
significant range of flat plate tests and an exhaustive study of 
earlier flat plate tests, Hughes (1954) proposed the following 
formula for the flat plate friction coefficient: 


Cy = 0.066/log(Rx) — 2.03) 


This formula was modified by the 1957 ITTC to that used in 
ship hydrodynamics today. The Hughes friction line gives 
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Table 11 Differences between 1957 ITTC and Schoenherr 
(ATTC) friction formulation as a function of Reynolds 
number 
Reynolds number C, (ITTC) C, (ATTC) ¢,(ITTC) /C,(ATTC) 
1x 10 0.008333 0.007178 1.1607 
2 0.006882 0.006137 1.12214 
3 0.006203 0.005623 1.1031 
4 0.005780 0.005294 1.0918 
5 0.005482 0.005057 1.0840 
6 0.005254 0.004875 1.0777 
7. 0.005073 0.004727 1.0732 
8 0.004923 0.004605 1.0691 
9 0.004797 0.004500 1.0660 
1 x 10€ 0.004688 0.004409 1.0633 
2 0.004054 0.003872 1.0470 
3 0.003741 0.003600 1.0392 
4 0.003541 0.003423 1.0345 
5 0.003397 0.003294 1.0313 
6 0.003285 0.003193 1.0288 
7 0.003195 0.003112 1.0267 
8 0.003120 0.003044 1.0250 
g 0.003056 0.002985 1.0238 
i x 10 0.003000 0.002934 1.,0225 
2 0.002669 0.002628 1.0156 
3 0.002500 0.002470 1.0121 
4 0.002390 0.002365 1.0106 
is) 0.002309 0.002289 1.0087 
6 0.002246 0.002229 1.0076 
7 0.002195 0.002180 1.0069 
8 0.002162 0.002138 1.0112 
g 0.002115 0.002103 1.0057 
1 x 10 0.002083 0.002072 1.0053 


slightly lower values than the ATTC friction line. Some hydro- 
mechanics laboratories, particularly those testing smal] mod- 
els, still use the Schoenherr line, as do many aerodynamic 
laboratories. Sometimes a variation of the Schoenherr line is 
used, yielding still lower C; values, such as the Hama (1954) 
friction formulation 


1C,)” = 3.46log(Rx) — 5.6 


Granville (1956) gives a useful summary of the various friction 
formulations proposed prior to 1956 and their validity. 

It follows that when adopting form factor values from partic- 
ular experimental studies, it is important to also adopt the 
same friction coefficient formulation. A form factor of 1.1 
found from low-speed tests, using the ITTC friction coefficient 
formulation, in the Reynolds number range around 1 x 10°, 
corresponds to a form factor of about 1.17 when using the 
ATTC friction coefficient formulation. 


USE OF VARIABLE WETTED AREA AND LENGTH IN MODEL 
TESTING 


The technique of determining wetted area and wetted 
length after each run in the test tank from underwater photo- 
graphs or video is not mentioned in the body of the paper. 
The author agrees with Mr. DeBord that this technique, in 
then being able to base the calculation of frictional resistance 
on the actual wetted area, for whatever Froude number, lee- 
way or heel angle under consideration, wil] allow for a more 
accurate separation of the viscous and wave resistance due 
to heel. It will also allow for a more accurate determination 
of the wave resistance as a function of speed, when upright. 
When needing an accurate value of the viscous resistance 
to subtract from the total resistance measured, or whenever 
carrying out fundamental research projects in which the fric- 
tional resistance and form values must be determined more 
accurately than can otherwise be carried out, this technique 
should be adopted. 

It should be noted however that in many model test facilities 
the requirement of making underwater photographs or video 


Predicting the Speed of Sailing Yachts 


will Jead to a significant cost increase which, in cases involving 
projects not requiring this accuracy, cannot be justified. 


EFFECT OF WAVES ON SIDE FORCE AND RESISTANCE DUE 
TO LIFT 


Very few studies have been carried out on the effect of 
waves on side force and induced resistance. Most studies 
involving waves have addressed the added resistance problem 
only, and on this topic a significant amount of data are now 
available. 

Generally, the results of model tests indicate that the effect 
of waves on side force is small. The results published by 
Kapsenberg (1991) indicate that for a 12-Meter yacht both an 
increase and a decrease in side force can occur of up to 
about 10%, although the average value for the considered 
wavelengths revealed a decrease in side force. This is contrary 
to the result given by Pedrick (1974), who found, for another 
model and in other wave conditions, that waves generally 
increased the side force. 

Changes in the level of the induced resistance due to waves, 
are heavily obscured by the relatively large added resistance 
component which comes into play when carrying out tests in 
waves. There is evidence to suggest, however, that the induced 
drag in waves is. completely in line with the change in the 
value of the square of the side forces, which relationship is 
found from tests in calm water. This supposition was adopted 
by Pedrick (1974) who found that on correcting the induced 
resistance in this way for the effect of waves, the remaining 
added resistance was independent of leeway angle, which 
could therefore be considered to be added resistance due to 
waves alone. 


DETERMINATION OF ADDED RESISTANCE IN BOW-QUARTER- 
ING SEAS 


Mr. DeBord also raises the question of how the level of 
the added resistance in head seas relates to the level of added 
resistance in bow-quartering seas and how this affects the 
prediction of sailing yacht performance when carrying out 
head sea tests only. It should be noted in this respect that 
only a few test facilities in the world allow for testing a model 
in bow-quartering seas, so this question is important to those 
carrying out studies of the added resistance of sailing yachts. 
Note also that when wind and sea directions coincide, a yacht 
never sails in head seas. 

Reference should be made here to work by Spens et al 
(1967), Pedrick (1974), Kapsenberg (1991), and others, who 
carried out tests in bow-quartering seas rather than in head 
seas. From this work, and from recent work involving calcula- 
tions of the added resistance in waves by Gerritsma et al 
(1993) and Sclavounos et al (1993), it is apparent that the head 
seas case constitutes the wave direction in which the added 
resistance is greatest. As the wave direction changes from 
the head seas to beam seas, the corresponding change in 
magnitude of the added resistance is found to be almost pro- 
portional to a cosine function. 

It follows that in carrying out tests in head seas, the results 
must be reworked to obtain the relative lower added resistance 
value for the appropriate bow-quartering seas condition. 
When it is assumed that wind and wave directions coincide, 
it is important to have carried out enough study of the perform- 
ance to know the true wind angle at which the best speeds 
to windward are obtained, to enable the added resistance to 
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be determined at these angles and added to the still water 
resistance. This exercise will result in a somewhat greater 
true wind angle for optimum performance. 

It is also important to take into account any wave spreading 
that might occur in the geographical area for which the per- 
formance of the yacht is to be optimized. Wave spreading also 
substantially reduces the magnitude of the added resistance, as 
demonstrated by Gerritsma and Keuning (1988). 


STATIONARY THRUST MEASUREMENT AT DOCKSIDE 


The author wishes to thank Mr. Asker for his contribution. 
The method of measuring the thrust in a line holding a yacht 
or vessel at a mooring is indeed a useful technique. Although 
strictly not covered by the subject matter of the paper, this 
technique can be used to gather data on the thrust produced 
by a particular set of sails or types of sails and, as such, is 
certainly worthy of mention. 


MINIMUM MODE SIZE IN RELATION TO VERY LOW UNCER- 
TAINTY EXPERIMENTS 


The author finally wants to thank Mr. Kirkman for his 
important contribution which covers much more ground than 
the speed prediction problem alone. His outline of the types 
of “Very Low Uncertainty” model tests that are now required, 
and the associated large size of models, is not presented from 
the viewpoint of speed prediction for any one design as a 
standard procedure, but rather from the point of view of 
validating specific CFD tools, investigating new scaling laws 
for Reynolds-dependent lift and resistance components, etc., 
necessitating nonconventional experiments. The author 
agrees that the types of tests mentioned, constituting the 
fundamental R&D in sailing yacht research today and tomor- 
row, will from now on become more important, particularly 
as the competition for the America’s Cup becomes more fierce 
and demanding. 


Is SAILING YACHT MODEL TESTING PART OF THE PAST? 


Mr. Kirkman proposes that model testing of the type cov- 
ered in the paper, purely from the viewpoint of speed predic- 
tion, is part of the past. While the author agrees that the 
demand for such model testing will become less and less as 
CFD becomes more accurate, it cannot yet be dispensed with 
at present. It is the author’s opinion that it will probably take 
another decade or so for all aspects of VPP models and CFD 
to become accurate enough so as not to have to carry out any 
model tests in important projects except for those of a very 
special nature as discussed above. 
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